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ABSTRACT 


An  experiment  was  performed  to  test  the  effectiveness  of  a  multivariate 
method  of  analysis  for  distinguishing  earthquakes  from  explosions.  The  data 
base  for  the  experiment  consisted  of  digital  recordings  made  at  twenty-seven 
stations,  including  four  large  arrays,  of  one  hundred  thirty-three  Eurasian 
events.  Spectral  magnitudes  were  measured  in  three  frequency  bands  of  the 
six  phases  P,  Lg,  long-period  P,  long-period  S,  LR,  and  LQ.  Complexities  of 
the  P— waves  were  measured  in  three  time  windows,  and  corner  frequencies  and 
low-frequency  spectral  levels  were  computed  for  the  P-wave  spectra.  When 
signals  could  not  be  detected,  spectral  magnitudes  were  measured  of  noise 
samples,  and  these  were  taken  to  be  detection  thresholds.  Signal  magnitudes 
and  detection  thresholds  from  all  the  stations  were  used  to  find  maximum- 
likelihood  estimates  of  the  spectral  magnitudes  for  each  event.  A  "training 
ggj-"  q£  presumed  earthquakes  and  explosions  was  used  to  find  discrimination 
variables  by  effectively  normalizing  the  network  spectral  magnitudes  to  elim¬ 
inate  bias  with  respect  to  m^. 

On  account  of  gaps  in  the  data,  many  different  combinations  of  different 
numbers  of  variables  were  used  in  the  classification  of  the  one  hundred  thirty 
three  events.  The  variables  which  were  applicable  to  the  most  events  were 
short-period  P-wave  spectral  ratios.  The  most  effective  discriminants  were 
spectral  ratios  of  Love  waves  to  P  waves. 

The  multivariate  discrimination  resulted  in  certain  misclassif ications 
which  were  attributable  to  insufficient  data  with  acceptable  signal- to-noise 
ratios.  The  misclassification  of  certain  deep  earthquakes  as  explosions  may 
be  the  result  of  systematic  errors  in  computing  the  magnitudes  of  deep  events, 
or  it  may  reflect  actual  differences  between  deep  earthquakes  and  shallow  ones. 
Both  deep  earthquakes  and  multiple  explosions  are  prone  to  misclassification  on 
account  of  anomalous  complexity  values.  The  overlap  between  the  earthquake 
and  explosion  populations  might  have  been  reduced  had  magnitude  corrections 
and  relative  weights  been  applied  to  the  data  from  individual  stations  before 
network  estimates  of  the  magnitudes  were  computed. 

Discriminants  involving  only  short-period  data  tend  to  misclassify  ex¬ 
plosions,  but  this  situation  can  be  rectified  (at  the  expense  of  increasing 
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the  false  alarm  rate)  by  assigning  more  weight  to  the  explosions  in  the  train¬ 
ing  set.  Investigations  were  made  of  the  influence  exerted  upon  discrimination 
by  both  the  number  and  the  geographical  distribution  of  events  in  the  training 
set.  Because  the  earthquakes  and  the  explosions  occurred  for  the  most  part 
in  regions  distant  from  each  other,  caution  should  be  taken  in  applying  the 
results  of  this  experiment  to  the  discrimination  of  unknown  events,  especially 
those  with  epicenters  in  regions  not  sampled  by  this  data  base. 
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INTRODUCTION 


Although  many  studies  of  seismic  discrimination  have  been  carried  out  in 
the  past,  they  have  tended  for  the  most  part  to  concentrate  on  some  single 
technique,  such  as  differences.  The  purpose  of  the  experiment  discribed 

in  this  report  is  to  develop  a  systematic  technique  of  discrimination  based  upon 
many  different  potential  discriminants.  The  effectiveness  of  the  multivariate 
discrimination  method  is  then  to  be  tested  by  applying  it  to  a  large  body  of 
data  collected  for  events  which  encompass  a  wide  range  of  magnitudes  and  which 
occurred  in  diverse  tectonic  settings.  The  utility  of  the  individual  dis¬ 
criminants  will  then  be  examined  in  terms  of  the  number  of  events  of  the  data 
set  to  which  they  may  be  applied  and  their  power  to  distinguish  between  earth¬ 
quakes  and  explosions. 

The  individual  discriminants  employed  in  this  experiment  have  been  ex¬ 
amined  in  detail  in  previous  reports  (von  Seggern  and  Sobel,  1977;  Sobel  et 
al,  1977a,  1977b,  1977c;  Sobel  and  von  Seggern,  1976,  1978)  as  they  apply  to 


von  Seggern,  D.  H. ,  and  P.  A.  Sobel  (1977).  Study  of  selected  Kamchatka 
earthquakes  in  a  seismic  discrimination  contect,  SDAC  Report  No. 

XR-76-10,  Teledyne,  Geotech,  Alexandria,  Virginia. 

Sobel,  P.  A.,  D.  H.  von  Seggern,  E.  I.  Sweetser,  and  D.  W.  Rivers  (1977a). 
Study  of  selected  events  in  the  Pamirs  in  a  seismic  discrimination 
context,  SDAC  Report  No.  TR-77-3,  Teledyne  Geotech,  Alexandria,  Vir¬ 
ginia. 

Sobel,  P.  A.,  D.  H.  von  Seggern,  E.  I.  Sweetser,  and  D.  W.  Rivers  (1977b). 

Study  of  selected  events  in  the  Baikal  Rift  Zone  in  a  seismic  dis-  _ 

crimination  context,  SDAC  Report  No.  77-5,  Teledyne  Geotech,  Alexandria, 

Virginia. 

Sobel,  P.  A.,  D.  H.  von  Seggern,  E.  I.  Sweetser,  and  D.  W.  Rivers  (1977c). 

Study  of  selected  events  in  the  Caucasus  in  a  seismic  discrimination 

context,  SDAC  Report  No.  TR-77-6,  Teledyne  Geotech,  Alexandria,  Virginia. 

Sobel,  P.  A.,  and  D.  H.  von  Seggern  (1976).  Study  of  selected  events  in  the 
Tien  Shan  region  in  a  seismic  discrimination  context,  SDAC  Report  No. 
TR-76-9,  Teledyne  Geotech,  Alexandria,  Virginia. 

Sobel,  P.  A.,  and  D.  H.  von  Seggern  (1978).  Analysis  of  selected  seismic 
events  from  Asia  in  a  seismic  discrimination  context,  SDAC  Report  No. 
TR-78-5 ,  Teledyne  Geotech,  Alexandria,  Virginia. 
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selected  regions  within  Eurasia.  Although  this  experiment  was  in  many  ways 
a  logical  extension  of  those  earlier  reports,  it  was  necessary  to  a  cer¬ 
tain  extent  to  sacrifice  the  constraint  of  a  purely  regional  study  in  order  to 
analyze  all  the  chosen  events,  which  occured  in  many  different  parts  of  the 
Eurasian  "Area  of  Interest"  (AI) .  The  multivariate  analysis  used  for  this 
experiment  was  applied  previously  to  earthquakes  and  explosions  in  the  South¬ 
western  United  States  by  von  Seggern  and  Rivers  (1979).  That  study  analyzed 
a  smaller  data  base,  consisting  only  of  events  whose  classification  was  known 
a  priori;  in  this  experiment  a  known  "training  set"  is  used  to  develop  dis¬ 
criminants  which  may  then  be  applied  to  the  classification  of  events  regarded 
as  unknowns. 


von  Seggern,  D.  H. ,  and  D.  W.  Rivers  (1979).  Seismic  discrimination  of 
earthquakes  and  explosions  with  application  to  the  Southwest  United 
States,  SDAC  Report  No.  TR-77-10,  Teledyne  Geotech,  Alexandria,  Virginia. 
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DATA  BASE 


Records  were  analyzed  for  133  Eurasian  events,  most  of  which  occurred 
m  the  last  months  of  1977.  These  events  are  listed  In  Tah.e  I,  alony  with 
their  epicenters  red  origin  times,  which  were  taken  from  two  separate  sts 
furnished  to  us.  The  two  sets  of  epicenters  were  frepuentlv  at  variance, 
with  the  result  that  several  events  were  listen  bv  the  NFIE  as  occurr  ng 
outside  the  Area  of  Interest  which  the  data  base  had  been  chosen  to  samp  e 
(see  Eig.  1).  The  epicenters  must  then  be  regarded  as  too  prone  to  error 
to  permit  discrimination  to  be  carried  out  on  the  basis  of  geography.  Ha 
such  a  geographical  discriminant  been  Judged  reliable,  several  events  cou 
have  been  eliminated  from  further  consideration  because  they  were  locate 
either  outside  the  political  boundaries  of  the  Soviet  Union  and  China  or  in 
deep  water  on  either  side  of  the  Kuril  Islands.  Some  possible  implications 
of  these  epicentral  differences  will  be  discussed  in  the  section  on  magnitude 
calculation.  A  table  of  depths  for  the  events  will  be  presented  later,  when 
the  effect  of  using  such  information  for  discrimination  is  evaluated. 

The  previously  mentioned  Eurasian  studies  of  Sobol  and  von  Seggern  showed 

the  importance  of  tectonic  setting  upon  seismic  discrimination.  For  insta^’ 

deep  earthquakes  occurring  on  a  descending  plate  behind  an  islan  arc  ar 

efficient  at  generating  surface  waves  than  are  the  shallow  earthquakes  which 

occur  in  front  of  the  island  arc  at  the  trench.  Teleseismic  P-waves  from  a 

region  characterized  by  a  low-velocity  zone  in  the  upper  mantle  are  more 

strongly  attenuated  at  high  frequencies  than  are  signals  from  regions  without 

such  a  zone.  Strike-slip  and  dip-slip  motions  along  faults  di  cr  in 

SUCtl  3  r  rQ„  These  and  other  differences 

efficiency  with  which  they  generate  surface  waves. 

between  earthquakes  in  different  regions  tend  to  make  more  i  ten 
bleros  of  distinguishing  earthquakes  from  explosions.  The 

fore  separated  into  groups  occurring  in  ten  geographical  regions  of  o  a 

indicated  in  the  last  column  of  Table  I.  These  regions  are  defined 

The  data  used  in  the  experiment  were  recorded  at  a  network  of  27  stations, 
which  are  listed  in  Table  III.  Only  digital  data  were  analyzed  For  eac , 
event  and  each  station,  the  records  for  as  many  as  six  phases  (short-perio 
and  Eg;  long-period  P,  S.  ER,  and  LQ)  were  made  available  for  analysis  In 
practice,  fewer  records  were  available  because  the  Eg  window  was  retained 
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ABLE  I  :  EVFN75  INCLUDED  IN  m|'  E  DT  SCI  T  H  T  S  A -T  ON  EX^  EF  T  M  EN 


EVENT 

1 

3 

4 
6 

7 

8 
9 

10 

14 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 


4  1 

45 

46 

47 

48 

49 

50 
53 

55 

56 

57 

58 

59 

60 
6  1 
62 

63 

64 

65 

66 

67 

68 

69 

70 

72 

73 

74 

75 

76 

77 

78 

79 

80 
8  1 


DATE 

26 

JUL 

77 

1 

NOV 

77 

1 

MOV 

77 

4 

NOV 

77 

4 

NOV 

77 

5 

NOV 

77 

5 

NOV 

77 

6 

NOV 

77 

30 

JUL 

77 

10 

a  nr. 

77 

17 

AUG 

77 

20 

AUG 

77 

1 

SEP 

77 

5 

S  E15 

77 

10 

SEP 

77 

30 

SEP 

77 

10 

NOV 

77 

10 

NOV 

77 

12 

NOV 

77 

12 

NOV 

77 

13 

NOV 

77 

15 

NOV 

77 

17 

NOV 

77 

18 

NOV 

77 

18 

NOV 

77 

18 

NOV 

77 

9 

OCT 

77 

18 

NOV 

77 

18 

NOV 

77 

16 

OCT 

77 

18 

NOV 

77 

16 

OCT 

77 

18 

NOV 

77 

13 

OCT 

77 

19 

NOV 

77 

20 

NOV 

77 

16 

OCT 

77 

19 

OCT 

77 

19 

OCT 

77 

20 

OCT 

77 

29 

OCT 

77 

26 

OCT 

77 

26 

OCT 

77 

26 

OCT 

77 

27 

OCT 

77 

28 

OCT 

77 

29 

OCT 

77 

29 

OCT 

77 

29 

OCT 

77 

30 

OCT 

77 

31 

OCT 

77 

20 

NOV 

77 

20 

NOV 

77 

20 

NOV 

77 

20 

NOV 

77 

21 

NOV 

77 

22 

NOV 

77 

22 

NOV 

77 

22 

NOV 

77 

22 

NOV 

77 

23 

NOV 

77 

26 

NOV 

77 

26 

NOV 

77 

27 

NOV 

77 

27 

NOV 

77 

28 

NOV 

77 

30 

NOV 

77 

OFIGIN  mINF 


16 

59 

57.6 

3 

54 

26. 0 

17 

56 

42.7 

10 

5  1 

40.0 

23 

54 

44. 7 

2 

9 

38.0 

4 

6 

58.5 

13 

3  1 

4  1.1 

1 

56 

58.0 

21 

59 

58.7 

4 

26 

57.7 

21 

59 

58.7 

2 

59 

57.5 

3 

2 

57.8 

16 

0 

3.3 

6 

59 

55.6 

4 

57 

46.0 

9 

22 

68.0 

5 

9 

16.0 

12 

27 

5.  3 

21 

2 

29.  3 

20 

20 

46.6 

4 

23 

54.0 

5 

20 

11.3 

5 

33 

19.7 

1  1 

27 

27.6 

1  1 

0 

0.  3 

15 

10 

4  1.8 

17 

23 

24.4 

20 

3 

35.0 

21 

55 

39.7 

15 

2 

48.8 

23 

12 

49. S 

20 

38 

52.  1 

11 

51 

14.2 

1 

4  1 

23.0 

21 

5 

37.3 

5 

1 

57.0 

21 

20 

41.6 

8 

18 

7.0 

3 

6 

57.7 

5 

38 

52.2 

7 

1  1 

31.3 

13 

14 

30.9 

7 

20 

28.9 

21 

15 

1.9 

4 

14 

56.5 

6 

26 

40.6 

10 

33 

56.7 

21 

38 

38.6 

9 

40 

38.5 

1 1 

1 

22.0 

18 

55 

28.0 

20 

57 

33.1 

23 

40 

35.9 

19 

43 

33.7 

0 

8 

3.7 

6 

56 

14.2 

1  1 

33 

45.0 

19 

16 

12.0 

10 

28 

7.0 

15 

44 

41.0 

22 

46 

52.2 

2 

9 

7.0 

3 

57 

0.0 

9 

2 

26.0 

4 

6 

57.5 

LATITUDE  LONGITUDE  F  EG  TIN 


69.  532  N 
55.436  N 
16.581  N 
23.000  N 
30.700  N 
42.900  N 

36.  6  35  N 
36.313  N 
49.777  N 
50.  923  N 
49.8  14  N 
64.223  V 
73.  376  N 
50.092  N 
57.294  N 
47.800  N 

37.  100  N 
33.000  N 
38.000  N 
35.757  N 
26.471  N 
38.172  N 
28.000  N 
32.693  N 
32.  640  M 
32.655  N 
73.  626  V 
32.719  N 

3  2.6  10  N 

4  8.  4  00  N 
60.081  N 
36.  466  ’I 
32.  698  M 
37.301  N 
36.450  N 
30.600  N 
49.4  19  N 
36.44  0  N 
49.  349  N 
56.426  M 
49.  84  1  N 
49.000  N 
46.400  N 
51.500  N 
53.  500  N 
39.569  N 
47.  00,0  N 
40.  283  N 
47.  137  N 
45.528  N 
55.  96  2  N 
56.800  K 
39.900  N 
37.419  N 
32.  395  N 
36.66  3  N 
36. 257  N 
36.478  M 
43.  000  N 
40.000  N 
34.000  M 
37.000  N 
39.465  N 
28.000  N 
50.000  N 
43.200  N 
49.957  N 


90.583  E 
130.545  E 
68.  745  E 
101.000  F 
8 1.279  E 

45. 300  F 
69. 128  E 
71.121  E 
78. 163  E 

110.761  E 
78.151  E 
99.577  E 
54.581  E 
78.96  1  E 
106.240  E 
48.145  P 

71.800  E 
89.000  E 
91.000  E 
71.331  E 
93.  1  1  4  E 

74.158  E 
90.000  E 
88.388  E 
88.431  E 
88.469  E 
53.  158  E 
88.316  E 
88.315  E 

152.900  F. 
143.373  E 
71.034  E 
88.422  E 
78.058  E 
7  1.279  E 
93.300  E 
155.  431  E 
71.273  E 
155.632  E 
164.133  E 
78.174  E 

155.800  E 

153.500  E 
153.400  E 
160.000  E 

73.406  E 

152.300  E 
63.455  E 

153.158  E 
145.988  E 
162.786  E 

108.500  E 
73.900  E 
71.824  E 
87.765  E 
71.206  E 
70.729  E 
71.236  E 
89.000  E 
75.000  E 
83.000  E 
71.000  E 

117.938  E 
90.000  E 
79.000  E 
47.600  E 
78.931  E 


14 

3 

10 

6 
7 
1  3 
10 
10 

9 

4 

9 

3 
14 

9 

4 

1  2 
10 

7 

6 

10 

7 

10 

7 

7 

7 

7 

14 

7 

7 

2 

3 

10 

7 

10 

10 

7 
2 

10 

2 

1 

9 
2 
2 
2 
1 

8 
2 

1  1 

2 

3 
1 

4 

8 

1  0 

7 

10 
10 
10 

8 
8 
7 

10 

5 
7 
9 

1  3 
9 
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TABLE  I  (CONTINUED) 


EVENT 

DATE 

0  R  TG IN 

TIME 

143 

2 

DEC 

77 

12 

57 

10.7 

144 

2 

DEC 

1  77 

16 

15 

44.9 

145 

3 

DEC 

77 

17 

6 

21.9 

146 

4 

DEC 

77 

4 

1  3 

47.0 

147 

4 

DEC 

77 

1  1 

39 

2.8 

148 

5 

DEC 

77 

23 

37 

37,  1 

1  49 

6 

DEC 

77 

10 

52 

53.5 

1  50 

151 

7 

7 

DEC 

DFC 

77 

77 

2 

16 

3 

1  9 

37.0 
3  3.9 

1  52 

153 

8 

8 

DEC 

DEC 

77 

77 

•> 

5 

2 

45 

54.0 

20.0 

1  54 

8 

DEC 

77 

1  3 

57 

4.4 

1  55 

8 

DFC 

77 

23 

37 

22.5 

1  56 

9 

DEC 

77 

4 

23 

36 . 0 

1  57 

10 

DEC 

77 

21 

58 

51.3 

158 

12 

DEC 

77 

1  1 

6 

42.0 

159 

13 

DEC 

77 

6 

58 

89.0 

160 

13 

DEC 

77 

1  1 

34 

20.0 

161 

162 

163 

1* 

15 

15 

DEC 

DEC 

DEC 

77 

77 

77 

5 

15 

15 

15 

7 

23 

42.  3 
51.8 
30.7 

164 

16 

DEC 

77 

7 

1 1 

4  1.6 

165 

16 

DEC 

~n 

9 

8 

5  0.7 

16s 

16 

DEC 

77 

10 

15 

27.5 

167 

16 

DEC 

77 

17 

55 

14.2 

16  8 

16 

DEC 

77 

23 

17 

17.7 

169 

19 

DEC 

11 

6 

57 

33.3 

170 

18 

DEC 

11 

16 

a  7 

17.  1 

171 

18 

DEC 

11 

19 

9 

18.7 

177 

18 

n?r 

11 

20 

43 

5.9 

173 

19 

nrc 

11 

1° 

1  2 

25.  0 

175 

20 

npr 

7  7 

-T 

27 

38.9 

176 

70 

DEC 

77 

20 

52 

1 0. 0 

177 

21 

DEC 

77 

0 

30 

46. 3 

179 

17° 

1  90 

21 

21 

21 

DFC 

DEC 

DFC 

77 

ln 

77 

16 

24 

20 

39 
17 

40 

3  3.0 
13.6 
5.0 

197 

22 

ppr 

77 

14 

5 

45.1 

193 

22 

DEC 

11 

19 

34 

5.0 

184 

1  P  5 

23 

23 

DEC 

DEC 

77 

77 

7 

a 

31 

0 

44.2 

54.1 

186 

24 

DFC 

77 

3 

27 

57.0 

187 

25 

DFC 

77 

8 

33 

37.0 

18R 

79 

DFC 

"T  *» 

17 

38 

4  2.0 

1  8° 

26 

DEC 

7  7 

4 

2 

57.-’ 

190 

26 

DEC 

77 

5 

15 

21.3 

191 

26 

DEC 

77 

23 

4 

34.0 

197 

27 

DFC 

77 

1 

10 

11.0 

19  3 

7T 

DFC 

11 

12 

3  1 

0.0 

194 

2  8 

DEC 

11 

1* 

14 

4  6.0 

1 0  5 

31 

DEC 

11 

3 

24 

3  8.6 

264 

20 

FA? 

76 

4 

3 

39.  3 

2  *  * 

70 

v  nr 

li 

3 

56 

57.7 

266 

19 

F  A- 

78 

3 

47 

0.0 

2*7 

1 1 

jnv 

7  9 

57 

0.0 

268 

28 

JOL 

78 

2 

47 

0.0 

7  r  0 

Q 

*vr. 

79 

1« 

0 

0.0 

777 

24 

MIG 

7  8 

18 

0 

0.0 

271 

29 

A  nr 

78 

2 

0 . 0 

272 

5 

FED 

7  8 

n 

22 

0.0 

2^3 

21 

r  Fr 

78 

1* 

0 

0.0 

274 

25 

r? 

7  6 

4 

5  9 

57.0 

275 

5 

re? 

iE 

4 

27 

0.0 

276 

13 

DFC 

7  5 

4 

5  6 

r  7  .  3 

277 

19 

!'  H  v 

if 

7 

56 

57.9 

279 

23  , 

TUT 

76 

2 

1 0 

50.0 

LATITUDE 

52.  929  !! 
44.  554  N 
41.932  N 

56.200  n 
48.250  N 
55.669  N 
41.425  K 
41.000  N 
35.65  1  N 
52.  900  11 
41.000  N 
50.  425  1! 
36.263  ll 
54.  400  11 
51.379  II 
51.400  N 
35.  400  11 

42.  300  11 
36.54  5  11 
4  3.  24  1  1! 

43.  643  u 
4  3.  233  1! 
51.630  N 
33.  305  1! 
36.991  N 
43.  059  1! 
55.278  II 
39.  865  11 
51.049  H 
39.733  M 
39.764  ri 
39.794  11 

55.200  M 
41.052  II 
52.999  h 
36. 196  n 
52. POn  u 
52. 956  n 
53.100  11 

44.6  10  II 
39.599  II 

51.200  II 
5 On  OOP  11 

40.900  11 
4  9.88  1  N. 
39.917  u 
33.300  n 
28.000  11 

54.7  00  ff 
56.0  00  11 

3  9.191  fi 
50.054  n 

4  9.790  II 
50.000  1! 
50.000  II 
50.000  N 
64.000  II 
66.000  N 
50.000  II 
43.000  II 
66.000  ri 
47.500  H 
*0.  000  N 
49.  799  11 
49.  856  11 
u?.  900  |! 


LONGITUDE 

R  EGK 

159.714 

E 

1 

146.459 

E 

2 

131.124 

E 

5 

163.100 

E 

1 

146.590 

E 

2 

162.21 1 

E 

1 

69.730 

E 

8 

72.000 

E 

8 

94.818 

E 

7 

8Q.700 

E 

9 

72.000 

E 

8 

149.808 

E 

2 

70.574 

E 

10 

160.600 

E 

1 

156.590 

E 

1 

157.500 

E 

1 

88.400 

E 

7 

133.200 

E 

5 

71.150 

E 

10 

45.168 

E 

13 

45. 372 

E 

13 

146.763 

E 

2 

159.460 

E 

1 

97. 560 

E 

7 

59.760 

E 

1 1 

47.082 

7 

13 

160.561 

E 

1 

77.334 

E 

P 

157.831 

E 

1 

77.696 

E 

8 

77.714 

E 

9 

69. 331 

E 

8 

159.200 

E 

1 

47.915 

E 

13 

159.904 

E 

1 

68.667 

E 

10 

159.500 

E 

1 

159.902 

E 

1 

163.400 

E 

1 

32.901 

E 

13 

77.449 

E 

9 

156.900 

E 

1 

91.000 

E 

q 

69.700 

E 

8 

73.  141 

E 

9 

7 1.967 

E 

7 

90. 900 

E 

7 

00.000 

E 

7 

16 1.500 

E 

1 

162.000 

41.100 
73.337 
78.149 
78.000 
79.000 
78.000  „ 

125.000  E 
112.000  E 
78.000  E 
89.000  E 
96.000  E 
47.500  V 
79.000  E 
78. 196  E 
78.007  E 

78.100  E 


1 

6 

9 

9 

9 

9 

9 

3 

4 
9 
8 

14 

12 

9 

9 

9 

9 
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(U)  Figure  la.  Overall  view. 

Location  of  events  used  in  the  discrimination  experiment. 
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TABLE  II 


(U)  Regions  within  the  Eurasian  Area  of  Interest 


imber 

Name 

Events 

1 

Kamchatka 

50 

58 

64 

143 

146 

148 

156 

157 

158 

165 

169 

171 

176 

178 

180 

182 

183 

186 

193 

194 

2 

Kuril  Islands 

36 

47 

49 

55 

56 

57 

60 

62 

144 

147 

154 

164 

3 

Siberia 

3 

18 

37 

63 

269 

4 

Lake  Biakal 

16 

21 

65 

270 

5 

Eastern  China 

77 

145 

160 

6 

Central  China 

6 

25 

195 

7 

Tibet 

7 

24 

27 

29 

30 

31 

32 

34 

35 

39 

46 

68 

75 

78 

151 

159 

166 

190 

191 

192 

8 

Tien  Shan 

59 

66 

73 

74 

i49 

150 

153 

170 

172 

173 

.175 

185 

188 

272 

9 

East  Kazakh 

14 

17 

20 

53 

79 

81 

152 

187 

189 

264 

265 

266 

267 

268 

271 

275 

276 

277 

278 

10 

Pamirs-Hindu  Kush 

4 

9 

10 

23 

26 

28 

38 

41 

45 

48 

67 

69 

70 

72 

76 

155 

161 

179 

11 

Turkmen 

61 

167 

12 

West  Kazakh 

22 

274 

13 

Caucasus 

8 

80 

162 

163 

168 

177 

184 

14 

Northwest  Russia 

1 

19 

33 

273 
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TABLE  III 


Stations  Used  in  the  Discrimination  Experiment 


Alaskan  Stations  (6):  AT-AK 

BF-AK 

CN-AK 

NJ-AK 

TN-AK 

UC-AK 

Special  Data  Collection  System  [SDCS]  (2):  HN-ME 

RK-ON 

Seismic  Research  Observatories  ISRO]  (10) :  ANMO 

ANTO 

BOCO 

CHTO 

GUMO 

MAIO 

NWAO 

SHIO 

SNZO 

TATO 

Alternate  Seismic  Research  Observatories  [ASRO]  (5):  CTAO 

KAAO 

KONO 

MAJO 

ZOBO 

Arrays  (4) :  ILPA 

KSRS 
LASA 
NORSAR 
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primarily  at  KSRS  only,  and  because  most  of  the  SRO  and  ASRO  records  were  of 
long-period  data  only.  The  short-period  SRO  records  consisted  mainly  of  de¬ 
tectable  signals,  whereas  the  long-period  records  were  retained  of  the  time 
window  in  which  a  given  signal  was  predicted  to  arrive,  even  if  no  signal 
could  be  seen  above  the  noise.  The  sampling  rate  for  long-period  data  was 
one  sample  per  second,  and  for  short-period  data,  except  for  NORSAR  and  LASA, 
the  rate  was  twenty  samples  per  second.  The  NORSAR  and  LASA  short-period 
data,  recorded  at  the  rate  of  ten  samples  per  second,  were  interpolated  to 
twice  that  rate  when  the  signals  were  processed. 

Long-period  data  at  NORSAR,  ILPA,  and  KSRS  were  recorded  at  seven  indi¬ 
vidual  sensors,  and  LASA  data  at  nine  sensors;  these  were  beamformed  using  the 
back  azimuth  for  the  event  being  processed  and  an  assumed  constant  value  for 
the  velocity  of  the  wave.  The  velocities  used  for  beamfonrrfng  are  given  in 
Table  IV.  In  everv  case  the  short-period  data  channel  with  the  largest 
signal-to-noise  ratio  was  chosen.  Short-period  data  for  KSRS  and  LASA  were 
taken  exclusively  from  the  first  sensor  in  the  record,  providing  that  channel 
was  not  malfunctioning. 

The  number  of  data  windows  processed  for  each  event  is  shown  in  Table  V. 

In  this  table  the  columns  headed  "LP"  and  "LS"  refer  to  long-period  P  and  S, 
respectively.  Many  more  seismograms  were  actually  available  than  were  processed, 
but  it  was  decided  that  certain  types  of  records  should  be  deleted  from  the 
data  set.  Principal  among  these  deleted  records  were  those  of  body  waves  for 
which  the  epicenter- to-station  angular  separation  A  exceeded  100°.  Deletion  of 
these  records  insured  that  no  waveforms  were  distorted  due  to  proximity  of  the 
propagation  path  to  the  earth's  core.  No  S-waves  were  processed  for  which  A 
<15°,  since  no  reliable  B-factor  could  be  assigned  to  them  and  since  the  long- 
period  S-wave  is  difficult  to  separate  from  surface  waves  at  short  distances. 

On  account  of  the  previously  mentioned  uncertainty  in  the  epicenters  of  many 
of  the  events,  magnitudes  determined  at  regional  distances  were  considered 
unreliable,  and  thus  signals  from  nearby  events  were  not  processed  unless 
their  amplitudes  could  be  determined  to  be  consistent  with  teleseismic  measu¬ 
rements.  Almost  all  the  records  deleted  for  this  reason  were  actually  of 
noise  rather  than  of  signal.  (The  reason  for  processing  noise  records  will  be 
explained  later).  A  particular  example  of  disagreement  between  regional  and 
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TABLE  IV 


(U)  Velocities  Used 

for  Beamforming 

Phase 

Distances 

Velocity  (km/sec) 

P 

A<20° 

111.2  / 

(14.80-0.1256  • 

A) 

P 

A>20° 

111.2  / 

(11.14-0.0712  • 

A) 

S 

A<20° 

111.2  / 

(26.98-0.2325  • 

A) 

S 

A>20° 

111.2  / 

(19.48-0.1118  • 

A) 

LR 

all 

3.6 

LQ 

all 

3.6 

TABLE  V 

NUMOEI  OP  DATA  HI NDOH5  \  N  A  I,  Y  7  E  D 


FVFNT 

p 

LG 

I  P 

LS 

I  R 

I.Q 

EVENT 

? 

LG 

I.  ° 

LE. 

1.P 

1 1 

1 

1  4 

0 

8 

6 

1  1 

1  1 

144 

4 

0 

10 

9 

1  1 

1  1 

3 

14 

1 

1  1 

1  1 

1  3 

1  4 

145 

1  1 

1 

12 

1  1 

1  2 

1  2 

4 

5 

2 

9 

8 

9 

10 

1  4b 

10 

1 

1  1 

1  2 

1  1 

1  1 

6 

; 

1 

1  1 

10 

15 

1  6 

147 

12 

1 

1  1 

1  1 

12 

1  2 

7 

) 

0 

7 

8 

12 

1  3 

148 

6 

1 

1  1 

6 

7 

7 

8 

1 

0 

8 

7 

1  1 

12 

149 

1  ) 

2 

9 

7 

i  ;> 

1  0 

4 

1  2 

0 

1  3 

10 

15 

15 

1  50 

1  3 

1 

1  2 

9 

1  i 

1  2 

1  0 

1  0 

2 

1  3 

1  1 

12 

1  3 

151 

1  1 

) 

1  1 

1  1 

1 1 

1  2 

1  4 

1 1 

1 

1  1 

0 

13 

1  3 

152 

12 

2 

1  1 

1  1 

1  3 

1  3 

1  6 

7 

1 

9 

8 

8 

8 

1  5  3 

9 

1 

1  1 

10 

1  1 

1  2 

1  7 

12 

1 

1  1 

8 

12 

12 

154 

1  1 

1 

1  1 

12 

1  ) 

1  3 

i  a 

6 

1 

7 

7 

7 

7 

1  55 

9 

2 

9 

5 

1  0 

0 

1  9 

n 

1 

8 

8 

1  2 

12 

156 

R 

1 

12 

10 

12 

1  2 

70 

4 

1 

9 

9 

9 

9 

157 

10 

1 

1  1 

12 

1  3 

1  3 

7  1 

12 

1 

9 

8 

a 

1  1 

158 

8 

1 

1  3 

14 

15 

1  3 

22 

6 

1 

4 

2 

a 

10 

159 

10 

2 

12 

1  1 

1  2 

1  2 

2  3 

9 

1 

1  0 

9 

10 

1  0 

160 

1  1 

1 

15 

13 

1  5 

1  4 

24 

1  1 

1 

1  1 

9 

10 

12 

16  1 

1  3 

2 

10 

8 

1  1 

12 

2  5 

1  1 

1 

12 

12 

14 

14 

162 

1  1 

2 

1  0 

a 

1  1 

1  1 

26 

1  3 

1 

1  3 

12 

14 

16 

163 

12 

1 

10 

8 

9 

12 

27 

1  1 

2 

9 

6 

15 

17 

164 

16 

1 

14 

12 

14 

1  3 

28 

5 

4 

1  0 

a 

12 

12 

165 

13 

1 

13 

12 

12 

1  1 

29 

2 

2 

10 

7 

13 

1  4 

166 

14 

1 

10 

1  1 

12 

12 

30 

1  3 

2 

12 

12 

15 

13 

167 

10 

1 

q 

6 

10 

1  0 

3  1 

15 

4 

14 

12 

14 

14 

168 

7 

0 

10 

8 

12 

1  1 

32 

12 

1 

12 

9 

13 

16 

169 

14 

1 

13 

1  2 

14 

9 

33 

12 

1 

9 

a 

12 

12 

170 

14 

2 

1  1 

9 

12 

10 

34 

6 

2 

9 

a 

10 

1  1 

17  1 

14 

1 

14 

12 

14 

1  4 

35 

14 

2 

1  1 

14 

15 

15 

172 

1  1 

2 

13 

1  1 

15 

1  3 

36 

4 

0 

10 

10 

12 

1  3 

173 

10 

3 

1  1 

8 

13 

12 

37 

12 

1 

1  1 

12 

14 

1  3 

175 

10 

1 

1  1 

6 

12 

12 

38 

1 1 

1 

1  1 

8 

8 

12 

176 

6 

1 

14 

12 

1  3 

14 

39 

14 

2 

14 

1  1 

16 

15 

177 

10 

2 

12 

7 

10 

1  3 

4  1 

5 

1 

8 

7 

10 

1  3 

178 

6 

1 

10 

12 

1  3 

1  3 

45 

8 

1 

13 

9 

13 

14 

179 

6 

1 

1  1 

8 

10 

1  1 

46 

1 1 

1 

12 

a 

14 

14 

180 

6 

1 

12 

1  1 

14 

14 

47 

14 

0 

1  1 

13 

15 

13 

182 

12 

1 

1  3 

13 

12 

14 

48 

8 

1 

10 

10 

10 

12 

183 

12 

1 

14 

1  3 

14 

1  3 

49 

12 

1 

9 

12 

12 

10 

184 

10 

1 

10 

a 

10 

9 

50 

14 

1 

1  1 

'12 

12 

1  1 

185 

1  1 

2 

10 

8 

12 

1  2 

53 

16 

1 

12 

10 

10 

9 

186 

4 

1 

9 

9 

1  1 

1  1 

55 

10 

1 

10 

10 

12 

10 

187 

9 

2 

9 

9 

10 

1  1 

5b 

2 

0 

10 

10 

1  1 

1  1 

188 

1  1 

2 

a 

7 

12 

12 

57 

10 

1 

12 

13 

1  3 

13 

189 

1  1 

1 

1 1 

9 

13 

1  3 

58 

10 

0 

1  3 

14 

12 

1  1 

190 

12 

3 

10 

7 

1  3 

1  3 

59 

7 

3 

1  1 

9 

13 

13 

19  1 

12 

1 

9 

8 

1  1 

1  1 

60 

7 

0 

14 

1  1 

12 

1  3 

1  92 

8 

0 

8 

8 

10 

10 

61 

12 

2 

10 

10 

13 

14 

193 

9 

2 

1  1 

10 

9 

10 

62 

13 

0 

1  1 

1  1 

12 

12 

194 

9 

2 

12 

10 

1  1 

10 

63 

12 

1 

14 

12 

1  4 

13 

195 

12 

2 

1  1 

1  1 

1  4 

15 

64 

12 

1 

1  3 

14 

13 

13 

264 

3 

1 

6 

6 

5 

6 

65 

12 

1 

15 

15 

16 

17 

265 

5 

2 

6 

5 

7 

8 

66 

9 

1 

1  1 

10 

14 

15 

266 

13 

2 

9 

10 

1  1 

10 

b  7 

1 1 

3 

8 

8 

9 

1 1 

267 

8 

1 

10 

9 

1  1 

1  1 

68 

14 

2 

1  1 

10 

13 

15 

268 

10 

1 

9 

9 

12 

12 

69 

1 

2 

7 

6 

9 

9 

269 

14 

1 

1  1 

12 

12 

1  3 

70 

10 

2 

12 

7 

15 

14 

270 

6 

1 

1  1 

1  1 

1  1 

12 

72 

10 

1 

10 

10 

12 

13 

271 

15 

2 

10 

10 

14 

12 

7  3 

8 

2 

1  3 

12 

14 

16 

272 

7 

0 

10 

10 

12 

12 

74 

3 

2 

8 

6 

9 

9 

273 

12 

1 

10 

12 

14 

14 

75 

9 

2 

1  1 

9 

1  1 

12 

274 

1 

0 

0 

0 

0 

0 

76 

1 1 

2 

12 

10 

12 

13 

275 

0 

0 

0 

0 

0 

0 

77 

12 

0 

1  3 

10 

13 

13 

276 

2 

0 

0 

0 

1 

2 

78 

8 

1 

7 

7 

11 

13 

277 

4 

1 

6 

4 

9 

7 

79 

6 

2 

1  1 

12 

7 

12 

12 

13 

8 

12 

13 

14 

10 

1  3 

12 

278 

4 

1 

6 

5 

7 

6 

80 

8  1 

143 

9 

12 

1  3 

0 

2 

1 

8 

1  1 

1  1 

TOTAL 

1275 

165 

1382 

1247 

1539 

1566 
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teleseisraic  magnitudes  is  that  of  the  high-frequency  component  of  the  short- 
period  P-waves.  Because  the  signals  from  regional  events  travel  only  through 
the  crust  and  not  through  the  mantle,  these  high  frequency  components  are  less 
attenuated  than  are  those  of  teleseismic  signals.  This  lack  of  attenuation 
introduces  a  distance-dependent  magnitude  bias  which  is  a  function  of  both 
geographical  region  and  frequency,  and  which  therefore  cannot  be  removed  by 
addition  of  the  standard  B-factor  for  P-waves.  This  problem  was  especially 
severe  at  KAAO  (for  which  many  of  the  AI  events  are  regional),  since  the  spectra 
of  all  events  close  to  that  station  were  observed  to  be  nearly  flat  due  to  the 

absence  of  attenuation,  and  hence  they  all  appeared  very  strongly  to  be  ex- 

\ 

plosion-like. 

Another  reason  that  certain  seismograms  were  deleted  from  the  data  set 
was  that  it  was  believed  that  the  signal  may  actually  have  arrived  in  a  time 
window  other  than  that  which  was  spanned  by  the  seismogram.  In  some  cases, 
the  signal  which  appeared  within  the  seismogram  was  believed  to  be  associated 
with  a  different  event.  Other  records,  which  included  visually  acceptable 
signals,  were  deleted  because  they  contained  too  short  a  noise  sample  before 
the  signal  arrival  to  permit  accurate  measurement  of  the  signal-to-noise  ratio 
[SNR].  Several  seismograms  were  deleted  from  the  data  base  because  the  ampli¬ 
tude  of  both  the  signal  and  noise  were  grossly  in  error,  usually  by  two  or 
more  orders  of  magnitude,  apparently  on  account  of  a  calibration  error.  With 
certain  rare  exceptions,  this  problem  was  confined  to  HN-ME  and  RK-ON,  the  two 
stations  for  which  calibration  data  had  to  be  entered  into  the  computer  by  hand 
in  order  to  create  the  digital  seismogram  tape.  Finally,  the  data  set  con¬ 
tained  nine  short-period  P  windows  from  GUMO,  and  they  were  all  deleted  because 
the  signal  and  noise  levels  were  both  quite  large  in  every  case.  A  short- 
period  P  wave  magnitude  bias  appeared  to  exist  at  NORSAR,  since  the  magnitude 
of  almost  every  event,  especially  at  high  frequencies,  was  noticeably  larger 
there  than  at  almost  all  the  other  stations  in  the  network;  nevertheless,  these 
signals  were  all  retained.  A  table  listing  the  relative  contribution  of  each 
station  to  the  data  base  will  be  presented  after  a  description  of  the  discri¬ 
mination  parameters  which  were  measured.  Table  V  shows  that  after  the  unusable 
seismograms  were  deleted  from  the  data  base,  no  seismograms  for  event  275  re¬ 
mained.  This  event  was  therefore  deleted  from  further  consideration. 
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DISCRIMINATION  PARAMETERS 


The  six  phases  (short-period  P  and  Eg,  long-period  P,  S,  LR,  and  LQ) 
which  comprised  the  data  for  each  station  and  each  event  were  processed  to 
compute  twenty-seven  parameters,  twenty-five  of  which  were  to  be  used  for 
discrimination.  These  twenty-seven  parameters,  listed  in  Table  VI,  are  now 
described  in  detail.  Whenever  no  signal  could  be  detected  in  a  given  data 
window,  the  noise  in  the  window  was  measured  as  if  it  had  been  a  signal, 
and  then  it  was  flagged  with  a  minus  sign  to  denpte  noise.  A  total  of  thirty 
measurements  were  actually  made  oik  the  signals,  but  the  conversion  from 
spectral  amplitudes  to  magnitudes  (to  be  described  later)  combined  the  three 
amplitudes  of  SH  and  of  SV  into  three  magnitudes  for  S. 

The  algorithm  for  the  computation  of  the  P-wave  spectral  parameters 
(numbers  4-9  in  Table  VI)  is  now  presented,  and  then  the  time-domain  analysis 
for  the  complexity  and  frequency— domain  analysis  for  the  spectral  amplitudoj 
will  be  described. 

Fitting  the  Observed  Spectrum  to  a  Source  Model 

Define  the  modulus  of  the  observed  ground  displacement  spectral  density 
as 


|AQ(f)|  =  |  fl  x(t)e_l2lTftdt| 


(1) 


Details  on  the  actual  computation  of  |A^(f)|  are  given  in  the  next  section. 
Hereafter,  the  bars  denoting  modulus  will  be  dropped,  and  the  amplitude  spec¬ 
trum  is  implied  throughout.  A  correction  for  the  instrument  response  1(f) 
gives 


A'  (f)  =  A  (f)/I(f)  (2) 

o  o 

This  spectrum  relates  to  the  source  spectrum  as 

A’(f)  =  A  (f)‘G(A)'e_TTft:*,E(f)  (3) 

o  s 

where  t*  is  the  attenuation  coefficient,  G(A)  is  the  divergence  term,  and  E(t) 
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TABLE  VI 


Parameters  Measured  for  Each  Station  and  Each  Event  (Whenever  Possible) 
Number 


Symbol 

P, 


V"2) 

fc(-2> 


t*(-2) 


V~3) 

fc(-3) 


Explanation 
low-  -i 

middle- 
high- 


frequency  P 


0.469  -  0.938  Hz 
1.094  -  1.876 
2.031  -  3.438 


low-frequency  spectral  level 
corner  frequency 

attenuation  coefficient 


assuming  P-wave  spectrum 

_ 2 

decreases  as  f 


9 

t*(-3)  . 

3 

10 

comp^ 

5-10  sec 

11 

comp  2 

5-15  sec 

12 

comp^ 

5-30  sec 

13 

Lgi 

low- 

14 

Lg2 

middle- 

15 

Lg3 

high- 

16 

LP1 

low- 

17 

LP2 

middle- 

18 

LP3 

high- 

19 

LSX 

low- 

20 

LS, 

middle- 

21 

LS3 

high- 

22 

LR1 

low- 

23 

lr2 

middle- 

24 

LR3 

high- 

25 

LQX 

low- 

26 

lq2 

middle- 

27 

LQ3 

high- 

same  as  4-6,  assuming  P-wave  spectrum  decreases  as  f 


coda  interval  complexity 


-3 


frequency  Lg 


^  period  P 


0.039  -  0.055  Hz 
0.063  -  0.086 

0.094  -  0.188 


>  frequency  long-period  S 


frequency  Rayleigh  wave 
(vertical  component) 


0.016  -  0.031  Hz 
0.039  -  0.053 

0.063  -  0.086 


-28- 


is  a  random  error.  Taking  the  logarithm  of  (3)  gives 

ln[A'(f)]  -  In  [A  (f )  *G(A)  ]-irft*  +  ln[E(f)]  (4) 

o  ® 

It  is  assumed  that  the  last  term  is  a  normally  distributed  variable  with  zero 
mean.  A  simple  source  spectral  shape  is  introduced  as 

f  . 

V£)  ■  °o  <rrr)  (5) 

c 

where  ft  is  the  long-period  spectral  level  at  the  source.  By  absorbing  the 

constant  G(A)  term  into  ft  ,  the  long-period  spectral  level  ft  '  at  the  station 

o  o 

rather  than  at  the  source  can  be  used  when  (5)  is  substituted  in  (4)  to  give 


f 

In [A^ (f )  ]  =  lnfij  +  b-ln(— £-)— Tift*  +  ln[E(f)] 


(6) 


c 


The  observations  A1 (f ) 

o 

fi0’  b’  fc’  3nd  pluS 
let 


are  thus  a  non-linear  combination  of  the  unknown  parameters 
a  random  error  term.  To  simplify  the  following  expression, 


Y  =  In  [  A '  (f)  ] 
i  o 

i 

f  . 

(7) 

F  =  lnftj  +  b  *  In  1— 

r-rVHv* 

(8) 

L  0/ 

PL  =  lnftj 

(9a) 

p2  =  f  c 

(9b) 

V» 

(9c) 

p.  =  -irt* 

(9d) 

The  iterative  least-squares  method  for  fitting  non-linear  functions  will 
be  used.  Starting  with  an  initial  estimate  p°  =  (p°,  p^,  P^»  P^)  the  un¬ 
known  parameters,  the  function  is  expanded  in  a  Taylor  series  about  these 
points: 


dF. 

l 


0  ,  0 

ln( 

'  P2 

-  n\  +  „°f 

P1  +  P3 

kfi  + 

o  1  +  p, t 
P  2/  4 

+  (Px  - 

pl) 

dFi  + 

dpl 

<P2  "  P2> 

+  (P3  - 

0* 

P3 

dFi  + 

dp3 

.  0 
<P4-  P4>  - 

dp, 

dF. 

l 

.  0 
bP,. 


(10) 
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The  derivatives  are  given  by 


(Ha) 


dFi,  P3fi 

,  0  0,.  .  0, 

dp2  P2(fi  +  ?2) 


?1  =  In/ 

,°  V 


fi  +  P2 


(lib) 


(He) 


(lid) 


Substituting  these  in  (10)  gives 


Fi  =  F1  +  (pl  “  Pl}  +  (p2  ‘  p2}  P3fl 

p2^Fi  +  P^) 


0 

+  (Po  ~  Po)  In  f  ?2  \  + 

Fi  +  p2 


(p4  -  P4>fi 


Let  be  the  residual  between  the  data  and  the  initial  model 


Ri  =  Yi  '  F1 


Since  plus  a  random  error  from  the  normal  density  N(0,o2) 

Ri  =  Fi  "  F1  +  N(0>°2) 
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Using  (12),  the  residual  becomes 


p  ,  ,  ,  0,  p3fi 

Ri  =  ^P1  "  pl)  +  ^p2  '  p2  _ 


p2(fi  +  p2> 


+  (P3  “  P3)  In 


+  (P4  -  P4)fi+N(o»a2) 


(15) 


Equation  (13)  is  termed  the  "equation  of  condition."  If  _n  points  in  the  spectrum 
are  available,  the  equations  of  condition  can  be  expressed  in  matrix  form  as 


Vi 


W  R 
n  n 


W, 


W 


W, 


0. 

P3fl 


p3(fi  +  P°) 


P3fn 


P?(f  + 

1  2  n 


0, 

P2 


W1  .  In 


W  .  In 


(rh) 

n  v2 


Wlfl 


W  f 
n  n 


APi 


APr 


Ap. 


Ap  / 


(16) 


where  arbitrary  weights  VL  have  been  entered.  Letting  the  residual  vector  be  R 
and  the  coefficient  matrix  be  B,  the  solution  vector  Ap  =  (Ap^,  Ap^,  Ap^,  Ap^) 
is  given  by 


Ap  =  (BTB)_1BTR  (17) 

The  solutions  Ap  are  added  to  p°  to  obtain  new  estimates  of  the  parameters, 
and  the  process  (13),  (16),  and  (17)  is  iterated  until  some  convergence  criterion 
is  satisfied. 

As  in  all  non-linear  fitting  problems,  proper  convergence  depends  on  a 
satisfactory  initial  estimate  of  the  parameters.  These  initial  estimates  must 
be  close  to  the  true  parameters  in  order  for  the  linearization  via  the  Taylor 
expansion  in  (10)  to  be  valid.  Preliminary  work  with  the  spectral  fitting 
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indicated  that  good  initial  estimates  were  not  possible  in  most  cases.  There¬ 
fore,  the  parameter  vector  £  was  reduced  by  assuming  a  range  of  values  for  P3 
and  P4  and  solving  the  non-linear  problem  for  the  remaining  parameters,  which 
were  long-period  spectral  level  and  corner  frequency,  with  each  set  of  (P3> 
p  )  values.  The  exponent  b(-P3)  for  decay  of  the  source  spectrum  was  set  at 
either  2  or  3.  These  are  the  exponents  found  in  most  source  model  derivations 
for  explosions  or  earthquakes  (von  Seggern  and  Rivers,  1979).  The  f  model 
is  favored  for  explosions  and  f"3  for  earthquakes  by  these  authors.  The  range 
of  t*  in  P4  was  assumed  to  be  between  0.0  and  1.0,  with  0.1  increments.  For 
each  combination  of  t*  and  b  then,  the  set  of  equations  of  condition,  reduced 


from  (16) ,  is 


W1R1 


W, 


0, 

P-f 


n. 


pJ(fH 


p2) 


(18) 


W  R 
n  n 


W 


n 


W 


p°Cf 


n 


Po> 


For  each  value  of  b,  the  assumed  t*  which  results  in  the  best  fit  of  the 
data  to  the  model,  as  determined  by  the  mean-square  of  the  residuals,  is  taken 
to  be  the  correct  estimate.  Note  that  under  the  two  assumptions  b  -  2  and  b  =  3 
the  t*  estimates  may  differ  greatly.  For  any  particular  spectrum,  convergence 

was  not  reached  for  most  of  the  trial  t*  values. 


An  initial  estimate  of  f£  =  1.0  Hz  was  found  to  be  satisfactory  in  most 
test  cases  and  was  used  throughout  in  computing  the  results  of  this  report  un 
less  otherwise  noted.  The  initial  value  of  tt'Q  was  taken  from  the  average  of 
the  spectral  points  in  A'q  (f)  between  0.5  and  1.0  Hz  for  a  given  signal. 
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Weighting  was  done  according  to  the  inverse  of  frequency,  that  is,  i\. 

The  number  of  points  n  of  the  spectrum  which  were  fit  by  (18)  was  var¬ 
iable,  being  determined  by  the  width  of  the  frequency  band  lying  between  the 
low- and  high-frequency  cutoffs  at  which  the  signal  spectrum  fell  below  the 
level  of  the  noise  spectrum.  These  cutoff  frequencies  are  shown  in  Figure  2. 

In  order  to  facilitate  the  convergence,  these  cutoffs,  and  sometimes  initial 
estimates  of  f  and  n'Q,  were  specified  for  certain  signals  rather  than  allowing 
them  to  be  calculated  by  the  computer.  The  convergence  algorithm  was  always 
subjected  to  the  further  constraint  that  0.2  Hz  <  f£  <  5.0  Hz.  Plots  were  made 
of  all  spectra  so  that  cases  requiring  specification  of  cutoffs  and/or  initial 
values  could  be  identified  visually. 

No  spectral  fit  parameters  were  assigned  to  those  spectra  for  which  the 
number  of  useful  points  n  was  less  than  six  or  for  which  no  convergence  was 

attained. 

Conplexity  of  the  Short-Period  P  Seismogram 

For  this  computation  a  window  about  the  short-period  P  wave,  sampled  at  20 
s/sec,  is  defined  from  10  sec  before  the  P  start  time  to  35  sec  after  it.  A 
bandpass  filter  with  corner  frequencies  at  0.5  and  5.0  Hz  and  rolloffs  of  12 
db/octave  is  applied  in  the  time  domain.  The  filter  output  points  are  squared 
and  smoothed  over  50  points  with  a  rectangular  sliding  window  to  produce  a 
smooth  envelope  trace,  and  the  square  roots  of  the  resulting  envelope  points 
are  then  taken.  Integration  is  performed  on  the  final  output  x  (t)  through 


with  integrated  values  defined  at  ^  =  200,  300,  400,  500  and  900  points  (10, 
15,  20,  25,  and  45  sec)  into  the  window.  Define  the  integrated  P-wave  signal 

as 


S  =  (A2  -  Ax)  -  1/2  A1 

The  second  term  here  represents  a  correction  for  the  noise  contribution  to  the 
signal  amplitude.  Now  define  three  coda  integrals  as 
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C1  =  (A3  -  A2)  -  1/2  Ax 

C2  =  (A4  '  A2)/2  '  1/2  A1 

C3  =  (A5  -  A2)/6  -  1/2  Ax 

The  ratio  of  coda  to  signal,  or  complexity,  is 

C  =  C./S  i  =  1,  3 

P.  i 

for  three  distinct  definitions  of  coda:  5-10  sec,  5-15  sec,  and  5-30  sec  after 
the  P  arrival  time.  If  the  complexity  was  calculated  to  be  negative,  it  was 
interpreted  as  measuring  noise  rather  than  signal  coda,  and  it  was  set  equal  to 
zero  so  that  it  would  be  ignored  by  further  processing. 

Spectral  Amplitudes  for  Short-Period  P 

For  the  short-period  spectra,  the  folding  frequency  is  10  Hz  since  the 
data  are  sampled  at  20  s/sec.  The  combination  of  instrument  response  and 
attenuation  in  the  earth  will  reduce  aliasing  effects  to  insignificant  levels 
in  nearly  all  cases. 

If  a  signal  is  present,  the  spectrum  is  computed  according  to  (1)  using 
an  FFT  algorithm  on  6.4  sec  (128  points)  of  data  immediately  following  the 
starting  point  of  the  signal.  A  noise  spectrum  is  computed  on  6.4  sec  of 
data  immediately  prior  to  this  starting  point.  The  mean  of  the  noise  from  the 
start  of  the  available  data  on  the  subset  tape  to  the  signal  start,  a  window 
of  varying  length,  is  used  to  compute  the  mean  of  the  trace.  This  mean  is 
subtracted  from  both  the  noise  and  the  signal  window,  and  a  half-cosine  bell 
function  is  applied  to  5%  of  the  data  window  at  both  ends  before  Fourier  trans¬ 
formation. 

The  resulting  65-point  spectrum  is  smoothed  once  with  Hanning  weights, 
and  the  instrument  response  is  removed  according  to  (2)  in  order  to  obtain  the 
ground  displacement  spectral  density  in  nm-sec.  The  average  spectral  density 
is  computed  for  three  bands: 
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EVENT  21  NJAK 


EVENT  21  NJAK 


Low  frequency  cutoff:  0  45  Hz 
High  frequency  cutoff  6  29  Hi 
fl  o  =  44.215  nm  »ec 
tc  =  1  457  Hi 


Selected  au toma tici II y 


fi  „  =  974.314  nm  sec 
tc  =  0  272  Hi 


LEGEND  Dotted  line  on  bottom  =  noise  spectrum 

Lower  solid  line  =  observed  signal  spectrum  corrected  for  instrument  respoiise  but  not  ettenuition 
Upper  solid  line  =  observed  signsl  spectrum  corrected  for  instrument  response  end  for  ettenuetion 
Dashed  line  -  theoretical  fit  to  observed  date 

Figure  2.  Example  of  the  P-wave  spectral  fit. 


0.46875  -  0.9375  Hz  (pts.  4-7) 

1.09375  -  1.8759  Hz  (pts.  8-13) 

2.03125  -  3.4375  Hz  (pts.  14-23) 

The  averages  in  the  three  signal  spectrum  bands  are  compared  to  those  in 

the  three  noise  spectrum  bands;  and,  for  any  band  where  the  noise  exceeds  the 
signal,  the  value  computed  for  the  signal  window  is  given  a  negative  sign  to 
indicate  a  noise  level  estimate.  If  the  trace  is  judged  to  have  no  visible 
signal  and  is  so  flagged  for  processing,  then  the  program  computes  only  the 
spectrum  for  6.4  sec  following  the  predicted  arrival  time  and  automatically 
appends  negative  signs  to  the  values  for  all  three  spectral  bands. 

Spectral  Amplitudes  for  Lg 

A  25.6  sec  (512  point)  window  for  Lg  is  processed  in  all  cases.  The  start 
time  of  the  window  is  picked  such  that  the  energy  arriving  with  a  velocity  of 
about  3.5  km/sec  is  at  the  start.  Mean  removal  and  tapering  is  done  for  the  P 
wave.  The  resulting  257-point  spectrum  is  reduced  to  65  points  by  successive 
applications  of  the  relation 

A'  (£i>  -  -r  A(f2i  -  2> +  -r  A« 2i  -  i>  +  ~r A  (C^  <19) 

Instrument  response  is  then  removed,  and  three  spectral  averages  are  computed 
over  bands  Identical  to  those  used  for  the  P  waves.  Flagging  of  noise  values 
is  done  in  the  same  manner  as  for  the  P  wave. 

Spectral  Amplitudes  for  Long-Period  P  and  S 

Since  the  long-period  data  is  sampled  at  1  s/sec,  the  folding  frequency 
is  .5  Hz.  Due  to  the  long-period  instrument  responses,  which  all  fall  off 
from  roughly  .05  Hz  with  slopes  of  at  least  18  db/octave,  aliasing  is  not 
expected  to  be  significant. 

Signal  windows  for  both  P  and  S  are  the  64  sec  (64  points)  following  the 
signal  start  time  while  noise  windows  are  the  64  sec  preceding  this  time. 

Mean  removal  and  tapering  is  done  in  the  same  manner  as  for  the  short-period 
P  wave.  The  resulting  33-point  spectrum  is  expanded  to  65  points  by  use  of 
the  relations 
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A'(fi>  -  A(f(±  +  d/2)  1-1.  3,  5,  ... 

a'(£,)  *  -rA(£i/2)  +-rM£i/2  + 1>  1  ■ 2- 4. 6-  ••• 


(20a) 

(20b) 


Instrument  response  is  removed,  and  three  bands  in  the  long  period  body 
wave  spectra  are  defined  as 

.0390625  -  .0546875  Hz  (pts  6-8) 

.0625000  -  .0859375  Hz  (pts  9-12) 

.0937500  -  .1875000  Hz  (pts  13-25) 

Average  values  of  the  points  in  each  band  are  computed.  Flagging  of  noise 
values  is  done  in  the  same  manner  as  for  short-period  P  waves. 

For  the  shear  wave,  the  horizontal  seismometer  components  are  rotated  to 
give  motion  perpendicular  (SH)  and  parallel  (SV)  to  the  projection  of  the  ray 
path  on  the  earth's  surface.  The  SV  component  was  measured  on  the  radial 

channel  only. 


Spectral  Amplitudes  for  Rayleigh  and  Love  Waves 

Variable  windows  are  set  for  the  Rayleigh  and  Love  waves  in  the  rotated 
long-period  traces  by  inspection.  In  general,  it  was  attempted  to  window  the 
entire  dispersed  waveforms  between  periods  of  50  to  10  seconds.  Due  to  back¬ 
ground  noise,  such  a  window  could  not  be  well  defined  in  many  cases,  and  the 
tendency  was  to  overestimate  the  window  length  required.  A  window  of  noise 
on  the  seismogram,  equal  in  length  to  the  LR  or  LQ  window  and  immediately 
preceding  it,  was  used  for  the  noise  spectral  estimate.  The  trace  mean  was  re¬ 
moved  as  previously  described  for  the  other  phases,  and  again  a  5%  cosine  taper 

was  applied. 

The  computed  spectra  for  LR  and  LQ  were  reduced  to  65  points  by  equation 
(19),  repeated  as  needed;  and  the  instrument  response  was  removed.  The  three 

spectral  bands  for  LR  and  LQ  were 

.0156250  -  .0312500  Hz  (pts  3-5) 

.0390625  -  .0546875  Hz  (pts  6-8) 

.0625000  -  .0859375  Hz  (pts  9-12) 

Again,  the  averages  of  the  spectral  amplitudes  over  the  points  in  each 
band  were  used  for  the  three  spectral  estimates.  Noise  estimates  for  when  no 
signal  was  visible  or  when  the  noise  spectral  estimate  exceeded  the  signal 
spectral  estimates  were  again  indicated  by  negative  signs. 
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DATA  CONTRIBUTION  OF  EACH  STATION 


Table  VII  shows,  for  each  of  the  twenty-seven  discrimination  parameters 
described  in  the  previous  section,  the  total  number  of  observations  that  were 
made  of  that  parameter.  The  listing  is  separated  in  order  to  show  the  number 
of  instances  in  which  the  measurements  were  of  signals  and  the  number  of  in¬ 
stances  in  which  they  were  of  noise.  It  can  be  seen  that  most  of  the  signal 
measurements  were  of  the  short-period  parameters.  The  contribution  of  each 
station  to  the  data  base,  both  of  signal  and  of  noise  measurements,  is  shown  in 

Table  VIII. 

By  summing  the  number  of  signal  and  noise  observations  of  short-period  P 
at  each  station,  as  given  by  Table  VIII,  it  is  possible  to  determine  not  only 
the  total  contribution  of  each  station  to  the  data  base  for  that  parameter  but 
also  the  relative  effectiveness  of  each  station  in  detecting  the  signals.  This 
procedure  is  illustrated  in  Table  IX.  TL  jse  stations  with  a  poor  detection 
rate,  such  as  HN-ME  and  RK-ON,  may  be  readily  identified.  It  is  also  apparent 
that  the  detection  rate  is  frequency-dependent;  in  particular,  P  waves  at 
NORSAR  have  a  slightly  greater  probability  of  being  detected  in  the  high-fre¬ 
quency  band  than  in  the  middle-frequency  b^nd,  but  at  most  other  stations  the 
middle-frequency  band  is  the  more  sensitive  one.  The  high  SNR  at  the  SRO  and 
ASRO  stations  leads  to  large  detection  probabilities  for  that  network,  but  it 
must  be  noted  that  the  total  number  of  signal  and  noise  observations  was,  for 
every  station,  less  than  one  half  of  the  133  events. 

Comparing  the  number  of  observations  (signal  or  noise)  of  P  with  the  per¬ 
centage  of  those  which  were  actually  signals,  as  shown  in  Table  IX,  suggests 
that  at  many  stations  seismograms  which  contained  no  detectable  P  signal  were 
systematically  excluded  from  the  data  base.  That  the  gaps  in  the  data,  at  least 
at  ANMO,  were  in  fact  weak,  unobserved  signals  is  shown  in  Figure  3.  This 
systematic  exclusion  of  weak  signals  would  result  in  a  bias  of  the  data  towards 
larger  magnitude  estimates  for  every  event.  As  will  subsequently  be  explained, 
measurements  of  the  background  noise  level  at  those  stations  which  did  not 


-38- 


lb 


TABLE  VII 


Total  Numbe  of  Observations  of  Each  Discrimination  Parameter 
Refer  to  Table  VI  for  Explanation  of  the  27  Parameters 


Parameter 

Signal 

Noise 

P 

556 

722 

1 

P2 

810 

468 

p 

684 

594 

3 

M-2) 

654 

0 

fc(-2) 

654 

0 

t*(-2) 

654 

0 

M-3) 

654 

0 

f  c(-3) 

654 

0 

t * (—3) 

654 

0 

comp^ 

466 

0 

comp  ^ 

428 

0 

comp 

312 

0 

J 

Lgl 

31 

136 

Lg2 

25 

142 

Lg. 

22 

145 

-3 

LP 

114 

1270 

1 

LP 

93 

1291 

L 

LP 

71 

1313 

J 

LS-. 

174 

1075 

1 

LS_ 

153 

1096 

2 

LS- 

127 

1122 

j 

LR 

291 

1249 

1 

LR„ 

438 

1102 

L 

LR 

320 

1220 

314 

1253 

J- 

lq7 

394 

1173 

L. 

lq3 

290 

1277 
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TABLE  VIII  (CUNTINBED) 
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Z0B0 


many 


detect  a  given  signal  were  used  in  the  maximum- likelihood  estimator  of  Ringdal 
(1976)  in  order  to  minimize  the  bias  which  would  affect  the  estimation  of  the 
event  magnitude  of  only  the  signals  were  used.  As  a  result  of  the  failure  in 
so  many  instances  to  include  in  the  data  base  those  seismograms  which  contained 
only  noise,  however,  the  full  power  of  the  maximum-likelihood  estimator  could 
not  be  used.  Had  the  situation  which  is  reflected  in  Figure  3  been  anticipated, 
the  resulting  bias  in  the  estimation  of  the  magnitudes  of  the  small  events 
could  have  been  at  least  partially  alleviated  by  the  substitution  of  average 
noise  levels  typical  of  each  station  in  place  of  the  missing  instantaneous 
measurements  of  the  noise  levels  at  the  predicted  signal  arrival  times.  Al¬ 
though  von  Seggern  and  Rivers  (1978)  showed  that  average  noise  levels  can  in 
fact  be  used  effectively  in  maximum -likelihood  magnitude  estimations,  this 
procedure  was  not  employed  in  this  experiment. 

In  contrast  to  the  situation  for  short-period  P,  long-period  seismograms 
such  as  those  for  LR  were  included  in  the  SRO  data  base  even  if  no  signal  was 
detected.  Table  X  shows  the  number  of  LR  seismograms  reported  from  each  station 
and  the  percentage  of  seismograms  which  contained  detectable  signals.  It  can 
be  seen  that  five  of  the  SRO  stations  reported  LR  seismograms  for  more  than  one 
hundred  events,  although  signals  were  detected  in  only  about  one  third  of  the 
cases.  It  was  thus  possible  to  utilize  more  fully  the  power  of  the  maximum- 
likelihood  magnitude  estimator  for  LR  (and  other  long-period  data)  than  for 
short-period  P.  The  LR  detection  probability  is  shown  in  Table  X  to  be  fre¬ 
quency-dependent,  LR^  being  detected  less  frequently  than  LR^  and  LR^.  It  is 
also  shown  that  certain  stations,  such  as  ILPA,  are  more  efficient  at  detecting 
LR^  relative  to  LR^  than  are  others,  such  as  KSRS.  Finally,  Table  X  shows  that 


Ringdal,  F.  (1976).  Maximum-likelihood  estimate  of  event  magnitude,  Bull. 
Seism.  Soc.  Am.,  66 ,  789. 

von  Seggern,  D.  H. ,  and  D.  W.  Rivers  (1978).  Comments  on  the  use  of  truncated 

distribution  theory  for  improved  magnitude  estimation,  Bull.  Seism. _ Soc . 

Am.,  68,  1543. 
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the  LR  data  base  was  heavily  influenced  by  the  inclusion  of  a  few  certain 
stations.  In  particular,  it  should  be  noted  that  was  detected  for  more  than 
forty  events  at  only  four  stations:  ILPA,  KSRS,  KAAO,  and  MAIO.  The  detection 
capability  of  these  stations  was  due  to  the  beamforming  capability  of  the  two 
arrays  and  to  the  proximity  of  three  of  the  stations  to  tectonic  regions  9  and 
10,  in  which  over  one  fourth  of  the  133  events  occurred.  Had  these  four  stations 
not  been  included  in  the  data  base,  the  number  of  detections  of  LR,  and  hence 
the  value  of  the  M  :m^  discriminant,  would  have  been  significantly  diminished. 
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CALCULATION  OF  MAGNITUDES 


The  amplitude  parameters  which  were  described  previously  were  measured 
for  each  signal  or  noise  window  and  then  were  converted  to  magnitudes  by  means 
of  formulas  which  will  now  be  described.  It  should  be  noted  that  the  amplitudes 
were  measured  in  the  frequency  domain  rather  than  the  time  domain  (units  of  nm- 
sec  rather  than  nm) ,  so  the  magnitudes  are  not  the  same  as  those  which  are  con¬ 
ventionally  calculated  from  the  visual  measurement  of  seismograms.  This  dis¬ 
tinction  should  create  no  difficulty,  since  it  is  the  difference  between  magni¬ 
tudes  for  different  parameters  and  different  events,  rather  than  the  magnitudes 
themselves,  which  is  important  for  discrimination.  The  same  formulas  and  B- 
factors  are  applicable  to  both  the  time-domain  and  frequency-domain  magnitudes, 
since  the  difference  between  the  two  scales  should  be  constant  for  each  para¬ 
meter. 

It  has  already  been  mentioned  that  the  small  distances  between  several 
events  and  some  of  the  stations  used  in  this  report  presented  certain  diffi¬ 
culties  for  the  calculation  of  magnitudes,  and  this  problem  will  now  be  described. 

The  observed  amplitudes  of  the  surface  waves  LR  and  LQ  were  converted  to 
magnitudes  by  correcting  for  attenuation  and  for  geometrical  spreading  of  the 
wavefront  (a  correction  which  employs  the  stationary  phase  approximation)  by 
means  of  the  formula  (cf.  Sato,  1967): 

Magnitude  =  log10  (amplitude)  +  log1Q  ^eYr J r  sin  (r/Re) 

where  r  is  the  epicenter-to-receiver  distance  in  km,  Re  is  the  radius  of  the 
earth,  and  Y  is  a  frequency-dependent  attenuation  coefficient  taken  from 
Mitchell  et  al  (1976)  and  listed  in  Table  XI.  The  second  term  on  the  right- 
hand  side  of  equation  (21)  is  plotted  in  Figure  4  as  a  function  of  the 

Sato,  R. ,  (1967).  Attenuation  of  seismic  waves,  J.  Phys.  Earth,  15 ,  32. 

Mitchell,  B.  J.,  L.  W.  B.  Leite,  Y.  K.  Yu,  and  R.  B.  Herrmann  (1976).  Atten¬ 
uation  of  Love  and  Rayleigh  waves  across  the  Pacific  at  periods  between 
15  and  110  seconds,  Bull.  Seism.  Soc.  Am.  66 ,  1189. 
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TABLE  XI 


Constants  Used  in  the  Calculation  of  Spectral  Magnitudes 


1)  LR  and  LQ  Attenuation  Coefficients 


frequency  (Hz) 
Y  (km--*-) 


band  1 

0.016-0.031 

0.0001 


band  2 

0.039-0.055 

0.00015 


band  3 

0.062-0.086 

0.0002 


2)  Takeoff  Angles  as  a  Function  nf  Distance 


V  \  '  ^  +  A3 

•  A3  +  A2  *  A 

2  +  A,  A  +  A  , 
1  o’ 

where 

A 

\ 

A3 

A2 

A1 

A 

0 

0.0°<A<  9.0° 

0.0 

0.0 

0.0 

-0.055556 

47.9 

9.0°<A<11.5° 

0.0 

0.0 

0.0 

-0.24 

49.56 

11.5°<A<15.5° 

0.0 

0.0 

0.0 

-0.475 

52.2625 

15.5°£A<20.0° 

0.0 

0.0 

0.0 

-1.644444 

70.388889 

20.0o£A<30.0° 

-8.9961xl0_4 

0.088244 

-3.1096 

45.407 

-188.74 

30.0°<A<85.0° 

-1.3478xl0-6 

3.0859xl0~4 

-0.025787 

0.025787 

23.882 

85.0°<A<100.0° 

0.0 

-7.8222xl0“4 

0.22419 

-21.429 

697.66 

-50- 


epicenter-to-receiver  angular  separation  A.  It  is  clear  from  the  figure  tnat 
small  changes  in  the  assumed  value  of  A  are  unimportant  for  teleseismic  events, 
but  at  distances  of  less  than  about  5°  they  are  sufficient  to  cause  substantial 
fluctuations  in  the  calculated  magntiudes.  It  was  on  account  of  this  sensi¬ 
tivity  to  errors  in  the  epicenter  location  that  LR  and  LQ  were  in  a  few  cases 
ignored  at  small  distances  if  the  calculated  magnitudes  conflicted  with  those 
which  were  measured  at  teleseimic  stations  or,  much  more  commonly,  if  only  noise 
was  observed  at  all  stations,  in  which  case  no  comparison  of  signal  magnitudes 
could  be  made. 

For  Lg,  magnitudes  were  calculated  by  means  of  the  following  formula: 

2 

Magnitude  =  log10  (amplitude)  +  log.^  r  •  (22) 

The  distance  dependence  of  this  empirical  formula,  plotted  in  Figure  5, 
is  stronger  than  that  of  equation  (21) ,  so  epicenter  mislocations  are  problem¬ 
atic  at  those  regional  distances  at  which  Lg  is  ordinarily  observed.  The  data 
base  of  Lg  detections  was  insufficient  to  test  the  validity  of  equation  (22), 
and  thus  it  is  not  known  how  much  error  was  introduced  into  the  magnitude  cal¬ 
culation  by  the  use  of  this  empirical  formula.  It  should  be  noted  that  Romney 

-3 

(1959)  deduced  an  amplitude  decay  proportional  to  r  for  Lg  propagating  across 
the  United  States. 

Body-wave  magnitudes  were  computed  from  the  formula: 

Magnitude  =  log^  (amplitude)  +  B(A)  (23) 

where  B  is  a  term  which  corrects  for  geometrical  spreading  of  the  wavefront  and 
for  anelastic  attenuation.  This  "B-factor"  is  of  course  dependent  upon  the 
depth  of  the  hypocenter,  but  since  all  133  events  were  to  be  regarded  as  pos¬ 
sible  explosions  only  B-factors  appropriate  for  surface  events  were  used.  ^ or 
the  short-period  P-wave  magnitudes  and  P3>  B-factors  were  taken  from 

the  study  of  Veith  and  Clawson  (1972),  which  used  conventional  measurements 


Romney,  C.  (1959),  Amplitudes  of  seismic  body  waves  from  underground  nuclear 
explosions,  J .  Geophys .  Res. ,  64 ,  1489. 

Veith,  K. ,  and  G.  Clawson  (1972).  Magnitudes  from  short-period  P-wave  data, 
Bull.  Seism.  Soc.  Am.  62,  435. 
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ORRECTION  ADDED  TO  LG  -  ALL  3  FREQUENCY  BANDS 


of  ra^  and  hence  frequencies  in  the  bands  P1  and  For  the  long-period  P- 

wave  magnitudes  LP^  LP2>  and  LP^  B-factors  were  taken  from  the  study  of 
Gutenberg  and  Richter  (1956),  which  used  P-waves  with  frequencies  of  about 
0.25  hz  (Gutenberg,  1945).  Although  this  frequency  is  higher  than  that  of 
even  the  high-frequency  band  LP^,  less  frequency-dependent  error  is  introduced 
into  the  magnitude  calculations  by  using  the  Gutenberg-Richter  B-factors  for 
the  long-period  bands  than  was  introduced  by  using  the  Veith-Clawson  B-factors 
for  all  three  short-period  bands.  That  this  is  so  may  be  demonstrated  by 
considering  the  fact  that  because  geometrical  spreading  of  the  wavefront  is 
independent  of  frequency,  the  frequency  dependence  of  the  shape  of  the  B(A) 
curve  is  determined  entirely  by  anelastic  attenuation.  [A  frequency-dependent 
geometrical  spreading  effect,  namely  P-wave  diffraction  by  the  earth  s  core, 
is  unimportant  since  no  measurements  were  made  on  body  waves  for  which  A  ex¬ 
ceeded  100°].  An  assessment  of  the  effect  of  attenuation  may  be  made  by  evalu¬ 
ating  the  quantity  exp  ( — tt f t * )  for  different  frequencies  f.  This  is  done  in 
Table  XII  for  an  assumed  value  of  0.5  sec  for  the  attenuation  coefficient  t*. 

It  can  be  seen  from  the  table  that  the  frequency  difference  between  the  long- 
period  bands  and  the  values  used  by  Gutenberg  and  Richter  is  unimportant.  It 
might  appear  that  B-factors  measured  at  frequencies  near  one  hz  are  inadequate 
for  the  high  frequency  band  P3,  but  it  will  be  shown  subsequently  that  the 
shape  of  the  B(A)  curve  (as  opposed  to  its  absolute  level,  which  may  be  con¬ 
sidered  arbitrary)  is  controlled  mainly  by  geometrical  spreading  rather  than 
attenuation,  so  it  is  not  strongly  frequency-dependent. 

The  B (A)  curve  is  shown  in  Figure  6  for  the  long-period  P  magnitudes  and 
in  Figure  7  for  the  short-period  P  magnitudes.  As  was  the  case  for  surface 
waves,  the  P-wave  magnitudes  are  insensitive  to  small  errors  in  epicenter 
location  if  the  source-to-receiver  separation  is  teleseismic,  but  small  location 
errors  can  lead  to  substantial  errors  in  the  magnitude  calculation  for  smaller 
distances.  This  problem  is  compounded  by  uncertainties  in  the  B(A)  curves 


Gutenberg,  B.  (1945).  Amplitudes  of  P,  PP  and  S  and  magnitudes  of  shallow 
earthquakes,  Bull.  Seism.  Soc.  Am.,  35 ,  57. 

Gutenberg,  B.,  and  C.  F.  Richter  (1956).  Magnitude  and  energy  of  earthquakes, 
Annali  Geofisica,  9_,  1. 
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TABLE  XII 


Relative  Effect  of  Anelastic  Attenuation  on  P-Wave 
Amplitudes  at  Different  Frequencies 


Frequency  Band 

lowest 

Frequency  (hz) 

highest 

Frequency  (hz) 

mean  F  (hz) 

-nft*|  .  „  c 
e  |t*=0.5sec 

LPi 

0.039 

0.055 

0.047 

0.929 

1 

LP9 

0.063 

0.086 

0.075 

0.890 

0.094 

0.188 

0.141 

0.801 

Gutenberg-Richter  B(A) 

(typically) 

0.250 

0.675 

Ft 

0.469 

0.938 

0.704 

0.331 

1 

P9 

1.094 

1.876 

1.485 

0.331 

2 

P3 

2.031 

3.438 

2.735 

0.014 

Veith-Clawson  B (A) 

(typically) 

1.000 

0.208 

-55 
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themselves  at  distances  of  less  than  about  25°.  Booth  et  al  (1974)  have 
measured  the  distance  dependence  of  P-wave  amplitudes  at  frequencies  of  1.0 
hz  and  0.063  hz,  and  their  results  show  certain  differences  from  the  B(A) 
curves  used  in  this  report  for  the  distance  range  0°£A<25°.  For  the  1.0  hz 
case  these  differences  are  less  than  about  0.3  magnitude  units,  but  for  the 
0.063  hz  case  the  difference  is  a  whole  magnitude  at  A  =  10°.  Discrepancies 
between  these  different  determinations  of  B(A)  for  a  given  frequency  are  pro¬ 
bably  due  to  lateral  heterogeneity  in  earth  structure,  and  so  they  are  depen¬ 
dent  upon  the  geographical  region  in  which  the  (short)  source-to-receiver  path 
occurs.  As  a  further  complication,  triplication  of  the  travel-time  curve 
causes  several  P  phases  to  arrive  close  together  for  small  A,  so  at  short 
distances  the  60  sec-long  window  which  was  used  for  long-period  P  contains 
several  different  phases  which  interfere  constructively  or  destructively  to 
produce  the  spectral  amplitude  which  is  measured  within  that  window.  A  final 
source  of  error  in  the  magnitude  calculation  is  introduced  by  the  assumption 
of  zero  depth  for  all  events;  as  can  be  seen  from  Figure  9  of  Veith  and 
Clawson  (1972),  the  error  caused  by  this  assumption  is  a  slowly  increasing 
function  of  depth  for  the  distance  range  30°<A<90°,  out  at  distances  of  less 
than  about  10?  the  error  increases  quite  rapidly  with  depth  below  the  base 
of  the  crust.  In  view  of  all  these  different  sources  of  possible  error,  it 
is  not  surprising  that  the  P-wave  magnitudes  measured  at  near-regional  and 
regional  distances  were  sometimes  in  conflict  with  measurements  made  at  tele- 
seismic  distances.  In  those  instances  in  which  such  disagreements  were  found, 
the  observations  made  at  the  close-in  stations  were  deleted.  Most  of  the 
deleted  close-in  measurements  of  P,  however,  were  actually  measurements  of  only 
noise.  Such  measurements  were  deleted  since  the  teleseisraic  observations  were 
also  of  noise,  and  no  comparison  of  signal  magnitudes  could  be  made. 


Booth,  D.  C.,  P.  D.  Marshall,  and  J.  B.  Young  (1974).  Long  and  short  period 
P-wave  amplitudes  from  earthquakes  in  the  range  0°-ll4°,  Geophys.  J.  R. 
Astr.  Soc. ,  39 ,  523. 
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An  example  of  the  arrival  of  more  than  one  P  phase,  in  this  case  for 
short-period  rather  than  long-period  P,  is  shown  in  Figure  8.  The  figure  shows 
the  P-wave  signature  from  event  14  as  recorded  at  stations  MAIO  and  KAAO,  for 
which  the  distances  from  the  epicenter  were  19.1°  and  16.6°,  respectively. 

It  might  appear  that  the  seismograms  show  an  earthquake  with  a  small  precursor, 
but  in  fact  event  14  is  a  presumed  explosion.  The  first  arrival  at  KAAO  is  due 
to  a  ray  which  penetrates  the  mantle  to  a  depth  of  192  km,  and  the  large  motion 
following  it  is  due  to  a  ray  which  penetrates  to  a  depth  of  428  km.  In  cases 
such  as  this,  the  6.4  sec  P-wave  windows  were  adjusted  to  begin  at  the  start  of 
the  second  arrival,  and  the  6.4  sec  noise  window  was  moved  back  so  that  it  did 
not  include  any  of  the  first  P  phase.  It  is  interesting  to  note  that  at  these 
two  SRO  stations  all  single  explosions  at  Semipalatinsk  look  like  multiple 
events. 

It  has  already  been  noted  that  no  records  of  long-period  S  were  analyzed 
for  source-to-receiver  separations  of  less  than  15°.  The  B-f actors  which  were 
used  for  S  were  taken  from  Gutenberg  and  Richter  (1956).  As  is  shown  in  Figure 
9,  the  calculated  magnitude  for  S  is  strongly  sensitive  to  location  errors  at 
distances  of  less  than  20°.  A  further  uncertainty  in  the  magnitude  is  intro¬ 
duced  by  use  of  the  formula 

Magnitude  =  y  log^  (SH^  +  SV^)  +  B(A)  (24) 

when  the  amplitudes  for  both  SH  and  SV  exceeded  the  threshold  SNR  of  1.5  (or 
when  they  were  both  less  than  this  threshold,  in  which  case  the  magnitude  is 
a  noise  level) ,  and  use  of  the  formula 

Magnitude  =  log^  (SH)  +  B(A)  (25) 

when  the  SH  amplitude  exceeded  the  threshold  and  the  SV  amplitude  did  not,  or 
use  of  the  formula 

Magnitude  =  log^  (SV)  +  B(A)  (26) 

in  the  opposite  case. 

A  rigorous  calculation  of  the  magnitude  of  the  low-frequency  spectral 
level  Q,0  can  be  made  by  means  of  the  formula 
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Figure  8.  P-wave  records  at  MAIO  and  KAAO  from  event  14. 


(27) 


Magnitude  =  log1Q  (amplitude)  +  |  log1Q  ^ 

where  the  second  term  on  the  right-hand  side  is  a  correction  for  the  geome¬ 
trical  spreading  of  the  wavefront  of  the  ray  which  traverses  the  source-to- 
receiver  angular  separation  &  and  then  emerges  at  angle  iQ  with  respect  to  the 
vertical  (Ben-Menahem  et  al,  1965) .  Implicit  in  this  formula  is  the  assump¬ 
tion  that  the  source  is  at  zero  depth,  and  hence  no  distinction  is  made 
between  angle  of  incidence  at  the  source  and  angle  of  emergence  at  the  receiver. 
The  emergence  angle  may  be  computed  in  terms  of  a  fourth— order  polynomial  in 
A,  as  is  shown  in  Table  XI.  The  results  of  this  polynomial  computation  are 
illustrated  herein  as  Figure  10.  The  term  ^  which  occurs  in  equation  (27) 
may  be  obtained  simply  by  taking  the  slope  of  the  curve  in  Figure  10,  and  the 
results  of  this  procedure  are  shown  in  Figure  11.  The  geometrical  spreading 
correction  may  then  be  computed,  and  it  is  shown  in  Figure  12.  A  comparison 
of  Figures  7  and  12  shows  that  the  shape  of  the  B(A)  curve  is  determined  pri¬ 
marily  by  geometrical  spreading,  as  was  asserted  previously.  The  curves  in 
these  two  figures  differ  significantly  from  each  other  (apart  from  an  ir¬ 
relevant  additive  constant)  only  at  distances  of  less  than  20°  or  greater 
than  95°.  In  these  same  distance  ranges  the  geometrical  spreading  correction 
must  be  regarded  as  being  poorly  determined,  and  it  was  decided  that  for  this 
report  equation  (27)  would  be  approximated  by: 

Magnitude  =  log1Q  (amplitude)  +  B(A)  (28) 

The  use  of  equation  (28)  might  seem  to  lead  to  large  errors  for  measurements 
of  made  at  near-regional  distances,  but  in  fact  no  theoretical  fit  was  made 
to  most  of  the  spectra  which  were  observed  at  such  small  distances.  The  reason 
that  the  theoretical  source  spectra  could  not  be  computed  is  that,  as  has  al¬ 
ready  been  noted,  the  P-wave  signals  which  traversed  such  short  and  mostly 
crustal  paths  were  characterized  by  quite  large  amplitudes  at  high  frequencies. 


Ben-Menahem,  A.,  S.  W.  Smith,  and  T.  L.  Teng  (1965).  A  procedure  for  source 

studies  from  spectrums  of  long-period  seismic  body  waves,  Bull. _ Seism. 

Soc.  Am.,  55,  203. 


-62- 


-65- 


GEOMETRICRL  SPRERDING  CORRECTION  TERN  RDDED  TO  LONG-PERIOD  ,_EVEl 


The  observed  spectra  thus  appeared  nearly  flat,  and  so  no  suitable  corner 
frequency  could  be  found  which  yielded  a  satisfactory  source  spectrum.  Com¬ 
puting  the  source  spectrum  was  made  even  more  difficult  by  the  imposition  of 
the  constraint  t*  >_  0.1,  so  the  case  of  vanishing  attenuation  could  not  be 
treated  satisfactorily.  As  a  result,  few  spectral  fits  were  made  at  near- 
regional  distances,  and  the  failure  of  equation  28  to  apply  at  such  distances  is 
unimportant.  A  larger  source  of  error  in  the  magnitude  calculation  at  small 
distances  is  the  assumption  of  zero  depth  for  all  events.  For  small  distances 
and  even  only  moderate  depths,  there  is  a  large  difference  between  the  take¬ 
off  angle  i^  and  the  angle  of  emergence  i0 ,  sc  equation  (27)  fails  to  describe 
adequately  the  geometrical  spreading  of  the  wavefront.  Even  if  all  the  events 
had  occurred  at  the  surface,  the  correct  values  of  i.  and  which  should  be 
used  for  any  particular  given  event  are  rather  uncertain  due  to  local  variations 
in  the  earth  structure  and  hence  in  the  raypath.  A  further  complication  in 
the  deter  mination  of  for  small  A  is  the  existence  of  triplications  in  the 

travel-time  curve  which,  as  was  shown  in  Figure  8,  leads  to  narrowly  separated 
arrivals  of  rays  having  different  paths  and  hence  different  amplitudes  (the 
lesser  of  which  may  not  be  seen  at  some  stations).  It  was  on  account  of  this 
uncertainty  that  the  approximation  given  by  equation  (28)  was  used  in  place 
of  equation  (27).  In  theory,  a  different  geometrical  spreading  correction 

could  be  determined  for  every  given  source  region-to-station  path,  using  values 
d! 

of  i0  and  — (calculated  from  travel-time  residuals)  which  taka  into  account 
the  effects  of  local  structure.  In  practice,  it  might  lead  to  less  error  in 
determining  simply  to  use  only  those  stations  for  which  A>30°,  although 
for  several  events  only  spectra  measured  at  shorter  distances  are  sufficiently 
noise-free  to  permit  the  calculation  of  the  theoretical  fit. 

"Magnitudes"  were  created  for  the  corner  frequencies  by  taking  the  log¬ 
arithm,  since  the  earthquake-explosion  discriminant  (to  be  illustrated  later) 

has  the  form  log  (ft  )  versus  log  (f  ).  The  values  of  the  attenuation  co- 

c 

efficient  c* ,  having  been  used  to  calculate  the  P-wave  spectral  fits,  are  not 
to  be  used  further,  so  no  processing  is  performed  upon  discrimination  para¬ 
meters  numbers  6  and  9. 

"Magnitudes"  were  computed  for  the  complexity  variables  simply  by  taking 
the  logarithms  of  the  observed  values.  During  this  procedure  all  complexity 
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measurements  made  in  certain  distance  ranges  were  deleted  from  the  data  set. 
First,  no  observations  were  retained  for  A<30°  since  in  this  range  the  P-wave 
coda  is  liable  to  contain  the  arrivals  of  other  P  phases  on  account  of  multi- 
pathing.  Second,  precaution  was  taken  to  avoid  the  inclusion  of  PcP  in  the 
toda  by  retaining  no  observations  of  the  short  complexity  window  comp^  for  A> 
78°,  of  comp ^  for  A>74°,  and  of  comp^  for  A>66°. 

It  has  been  mentioned  that  data  records  containing  only  noise  and  no 
signal  were  flagged  by  means  of  a  minus  sign  but  were  otherwise  treated  as  if 
they  had  been  signals.  The  data  base  for  each  event  thus  consisted  of  as  many 
as  twenty-seven  parameters  measured  at  as  many  as  twenty-seven  stations,  and 
eighteen  of  these  parameters  could  be  measurements  of  either  signal  or  noise. 
(For  noisy  P  waves,  the  spectral  fit  parameters  and  complexities  are  set  equal 
to  zero).  The  measurements  of  noise  levels  were  converted  to  magnitudes  in 
the  same  manner  as  if  they  had  been  signals,  and  they  were  then  assigned  a 
special  flag  so  that  they  could  be  distinguished  as  detection  thresholds  rather 
than  as  actual  detections. 

For  each  event,  the  magnitudes  at  the  individual  stations  were  next  com¬ 
bined  to  give  a  network  estimate  of  the  magnitude  for  each  parameter.  For 
the  low-frequency  spectral  levels,  corner  frequencies,  and  complexities  this 
network  estimate  consisted  merely  of  the  arithmetic  mean.  Network  averages  of 
t*  are,  of  course,  meaningless,  so  those  two  parameters  are  set  equal  to  zero. 
For  the  spectral  magnitudes  the  network  estimates  were  computed  using  the 
maximum-likelihood  technique  of  Ringdal  (1976).  This  technique  introduces 
the  noise  levels  into  the  computation  as  upper  bounds  for  the  magnitudes  of 
undetected  signals,  thereby  eliminating  the  magnitude  bias  which  would  result 
from  the  inclusion  of  only  those  signals  with  a  large  SNR.  Because  so  many 
observations  were  of  noise  rather  than  of  signals,  this  magnitude  bias  would 
have  been  severe  had  not  the  noise  measurements  been  made  and  the  Ringdal 
method  used.  There  were  many  cases  in  which  only  noise  levels  and  no  signals 
were  measured  at  all  the  stations  for  which  there  were  data  records.  In  these 
cases  we  have  taken  Ringdal's  formula  for  the  probability  of  non-detection  of 
an  event  of  a  given  magnitude  by  a  station  with  a  given  noise  level  and  have 
used  it  to  calculate  the  likelihood  that  an  event  of  that  magnitude  would  not 
have  been  detected  by  any  of  the  stations  in  the  network,  given  the  noise 
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levels.  By  incrementing  the  assumed  magnitude,  we  then  take  as  an  upper  bound 
for  the  event  magnitude  that  value  for  which  the  likelihood  of  no  detection 
equals  ninety  per  cent.  These  upper  bounds  are  flagged  as  such  so  that  they 
will  be  treated  separately  from  the  magnitudes  which  were  calculated  on  the 
basis  of  at  least  one  actual  detection.  In  order  to  carry  out  this  calculation 
of  the  upper  bounds  it  is  necessary  to  assume  a  value  of  the  standard  deviation 
which  characterizes  random  measurements  of  the  magnitude  of  the  parameter  in 
question.  We  have  assumed  in  all  cases  a  value  of  0.40,  since  this  value  is 
typical  of  time-domain  magnitude  standard  deviations  and  is  at  least  consistent 
with  the  values  found  for  the  frequency-domain  magnitude  in  those  cases  which 
included  detections. 

An  event-by-event  listing  of  the  network  estimates  of  the  magnitudes  for 
each  parameter  is  reproduced  as  Appendix  I.  It  should  be  noted  that  the 
maximum-likelihood  estimates  were  calculated  subject  to  the  constraint  that 
the  standard  deviation  of  the  most  likely  Gaussian  distribution  be  in  the  range 

0.15  <  o  <0.63  (29) 

The  number  of  signal  and  noise  observations  used  to  deduce  these  values  are 
shown  in  Appendix  II. 


DISCRIMINATION  VARIABLES 


After  the  deletion  of  the  attenuation  coefficients,  parameters  numbers 
b  and  9,  there  remain  twenty-five  parameters  for  which  the  network-averaged 
magnitudes  for  each  event  are  to  be  used  for  discrimination.  Discrimination 
which  is  based  upon  these  magnitudes  as  raw  data,  however,  can  suffer  from  a 
magnitude  bias  and  can  therefore  yield  misleading  results.  One  source  of  this 
bias  is  the  failure  of  certain  parameters  such  as  Mg  to  scale  with  equally 
for  both  the  earthquake  and  the  explosion  populations.  Another  source  of  bias 
is  the  difference  in  the  magnitude  distributions  for  the  two  populations,  a 
difference  which  is  illustrated  in  Figure  13.  The  criteria  for  assignation  of 
the  events  to  one  or  the  other  of  the  two  populations  will  be  explained  sub¬ 
sequently.  Asymmetry  between  the  two  distributions  can  lead  to  spurious  "dis¬ 
criminants  .  As  a  hypothetical  example,  consider  classical  M  -  m  discrimi- 

s  b 

nation  as  carried  out  upon  a  set  of  events  consisting  of  small  explosions  and 
large  earthquakes  (Figure  14).  Depending  upon  the  exact  shape  of  the  distri¬ 
bution  of  the  events  in  the  Mg  -  m^  plane,  a  straightforward  statistical  ana¬ 
lysis  might  determine  that  the  two  populations  of  events  were  best  separated 
by  a  discriminant  line  having  a  nearly  horizontal  slope.  The  criterion  upon 
which  future  events  would  be  classified  is  that  events  with  large  M  are  always 
earthquakes  and  events  with  small  Mg  are  always  explosions.  Although  this  hypo¬ 
thetical  example  is  an  extreme  one,  it  points  out  that  a  magnitude  bias  can  be 
built  into  a  discriminant  function  if  raw  data  are  used. 

The  remedy  for  this  bias  is  to  perform  discrimination  based  not  on  the 
individual  variables  Mg  and  m^  but  rather  on  the  single  variable  (M  -  m^) . 

In  this  manner  one  implicitly  imposes  the  constraint  that  the  Mg  versus 
discriminant  line  has  unit  slope.  If  this  line  has  some  (unknown)  intercept 
c  on  the  -  axis,  all  events  for  which  >c  lie  above  the  line  M  = 

c>  hence  they  are  earthquakes.  Similarly,  those  events  for  which 

Ms  ”  “b  <c  are  explosions.  The  statistical  analysis  then  consists  of  finding 
the  Ms  -  intercept  c  rather  than  both  the  intercept  and  the  slope.  Of  course 
one  may  specify  a  priori  the  intercept  also  (as  in  the  case  of  the  classical 
relation  Mg  =  +  1.5),  in  which  case  the  variable  (M&  -  n^)  may  be  used  as 

a  discriminant  without  reference  to  any  statistical  analysis  of  a  given  pop¬ 
ulation  of  earthquakes  and  explosions. 
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Figure  13.  Distribution  of  middle-frequency  P-wave  magnitudes 
for  explosions  and  earthquakes. 


Even  with  a  large  body  of  data  encompassing  a  wide  range  of  magnitudes 
for  both  earthquakes  and  explosions,  an  unrestricted  analysis  of  raw  data 
might  be  expected  to  exhibit  a  magnitude  bias.  Such  a  bias  would  tend  to 
occur  at  least  for  the  pair  of  variables  Mg  and  i^since  Mg  for  large  earth¬ 
quakes  tends  to  lie  above  the  discriminant  line  by  a  larger  amount  than  Mg  for 
large  explosions  tends  to  lie  below  it.  An  analysis  of  the  raw  data  might 
thus  attempt  to  classify  events  (particularly  if  Mg  was  frequently  not  detected 
for  small  events)  on  the  basis  of  a  discriminant  line  of  slope  greater  than 
unity.  We  therefore  believe  that  the  constraint  of  unit  slope  should  be  im¬ 
posed  upon  the  Mg  versus  discriminant  by  analyzing  data  in  terms  of  the 

single  variable  (M  -  m^) . 

We  have  extended  the  analogy  of  visual  (Mg  -  n^)  to  our  data  base  by 
"normalizing"  all  spectral  magnitudes  to  those  of  the  high-frequency  P-wave 
magnitudes  P3  by  subtracting  from  all  magnitudes  the  quantity  a»? y  where  a 
is  a  coefficient  which  is  to  be  determined.  The  short-and  long-period  P-wave 
magnitudes  P2>  P3,  LP^  LP2  and  U>3  are  thereby  converted  from  magnitudes  to 
spectral  ratios.  The  worth  of  these  spectral  ratios  to  the  problem  of  discri¬ 
mination  lies  in  the  relative  enhancement  of  high  frequencies  for  explosions. 
For  the  high-frequency  band  P3  this  enhancement  is  due  to  the  higher  corner 
frequencies  of  explosions  than  of  earthquakes  and  to  the  (admittedly  arguable) 
slower  high-frequency  decay  rate  of  explosions  than  that  of  earthquakes,  f 
versus  f“3.  With  signals  originating  in  so  many  diverse  tectonic  regions, 
however,  these  differences  may  be  overwhelmed  by  differences  in  t*.  The  lower 
frequency  bands  are  expected  to  contain  less  energy,  relative  to  that  in  P3, 
for  explosions  than  for  (deeper)  earthquakes  on  account  of  cancellation  of  the 
signal  by  pP.  Subtracting  from  the  S,  LR,  and  LQ  magnitudes  converts 

them  to  measures  of  the  relative  shear  wave-to-compressional  wave  excitation 
ratio  for  the  event.  This  quantity  is  expected  to  be  higher  for  earthquakes 
than  for  explosions,  but  we  note  that  this  will  be  affected  strongly  by  the 
focal  mechanism  and  depth  of  the  earthquakes. 

In  order  to  maintain  a  strict  analogy  with  classical  Mg  -  differences, 
all  magnitudes  should  have  been  normalized  to  the  middle-frequency  P-wave  mag¬ 
nitude  P  rather  than  to  Py  since  visual  measurements  are  usually  made  of 
signals  of  about  1  sec  period.  The  higher  frequency  band  was  chosen  instead 
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in  order  to  utilize  the  largest  possible  frequency  contrast  between  the  short- 
period  and  the  long-period  data.  As  a  result  of  this  choice,  the  spectral 
analog  of  the  versus  m^  discriminant  will  have  a  slope  which  is  steeper  than 
the  unit  slope  which  has  thus  far  been  assumed.  Normalization  of  all  magnitudes 
to  generates  two  short-period  discrimination  variables,  namely  the  P-wave 
spectral  ratios  of  P^  to  P^  and  of  ¥  ^  to  P^.  In  order  to  more  fully  charac¬ 
terize  the  shape  of  the  P-wave  spectrum,  which  is  assumed  to  be  a  valuable 
discriminant,  an  additional  spectral  ratio  was  used,  namely  that  of  P^  to  P,. 
That  the  three  short-period  discrimination  variables  are  not  redundant  can  be 
seen  by  considering  the  expressions  for  the  spectral  ratios: 


V 

V 

V 


1 

2 

3 


«13  P3 
“23  p3 
“12  P2 


(30) 


It  can  easily  be  shown  that  the  three  discrimination  variables  V^,  and 
are  linearly  independent  if  043  /  0132  *  “2  3* 


In  order  to  determine  suitable  values  of  the  coefficients  a,  we  have  plotted 
in  Figure  15  P^  as  a  function  of  P^  for  both  the  earthquake  and  explosion  pop¬ 
ulations.  Those  events  which  were  classified  a_  priori  on  the  basis  of  their 
epicenters  and/or  origin  times  as  being  explosions  were: 


Semipalatinsk  Area  -  14,  17,  20,  79,  81,  189,  266,  267,  268, 

276,  277,  278 
Novaya  Zemlya  -  19,  33 

West  Kazakh  (Volga  River/Caspian  Sea)  -  22,274 
Siberia  and  Lake  Baikal  -  1,  16,  18,  21,  269,  270,  273 


In  addition,  three  other  events  were  identified  by  epicenter  and  origin 
time  data  as  being  explosions,  but  inspection  of  their  P-wave  signals  revealed 
them  to  be  multiple  events.  These  presumed  shot  arrays  are  events  53,  265, 
and  271,  and  they  are  denoted  by  a  special  symbol  in  Figure  15.  All  other 
events  were  classified  a  priori  as  earthquakes.  It  has  already  been  noted  that 
events  occurring  in  different  regions  can  appear  dissimilar  on  account  of 
propagation  path  effects  even  if  their  source  mechanisms  are  identical.  There 
thus  exists  the  possibility  that  the  P^  versus  P^  "discriminant"  illustrated 
in  Figure  15  does  not  so  much  separate  earthquakes  from  explosions  as  it 
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separates  events  in  "explosion  regions"  (defined  above)  from  events  in  all 
other  regions  listed  in  Table  II.  This  possibility  will  be  investigated  further 
in  a  later  section  of  this  report.  The  same  separation  of  events  into  earth¬ 
quake  and  explosion  populations  was  used  in  Figure  16  to  determine  the  coeffici¬ 
ents  a  which  were  to  be  used  to  calculate  discrimination  variables  based  on  the 
normalization  to  of  the  other  fourteen  spectral  magnitudes. 

Many  data  points  in  Figures  15  and  16  have  arrows  attached  to  them  pointing 
either  down  or  to  the  left;  such  points  are  the  previously  mentioned  90%  con¬ 
fidence  limit  upper  bounds  to  the  true  magnitudes,  calculated  from  measurements 
of  only  noise  and  not  signal.  The  true  magnitude  of  each  point  is  assumed  to 
lie  in  the  direction  of  the  arrow.  Not  appearing  in  the  figures  are  those  points 
for  which  the  values  of  both  the  abscissa  and  ordinate  magnitude  variables  were 
approximated  by  upper  bounds,  since  these  points  cannot  contribute  any  useful 
information  to  the  discrimination  between  the  two  populations  of  events.  The 
dashed  lines  in  Figures  15  and  16  were  visually  fit  to  be  either  coincident 
with  or  parallel  to  other  lines  (not  shown)  which  as  nearly  as  possible  separate 
the  explosions  from  the  earthquakes.  The  two  populations  overlap  for  every  pair 
of  magnitude  variables  shown  in  the  figures,  so  in  no  case  was  perfect  sepa¬ 
ration  possible.  The  reason  that  the  lines  plotted  on  the  figures  are  in  some 
cases  parallel  to,  rather  than  coincident  with,  the  "best"  discriminant  lines 
is  that  it  is  only  the  slopes  of  these  lines,  and  not  their  intercepts,  which 
affect  their  discrimination  capability.  That  this  is  so  may  be  seen  from 
equation  (30),  which  may  be  interpreted  geometrically  in  terms  of  Figures  15 
and  16  as  meaning  that  the  value  of  the  discrimination  variable  for  a  given 
point  on  one  of  those  plots  is  given  by  the  vertical  distance  from  the  given 
point  to  the  dashed  line.  If  the  intercept  of  the  line  is  changed  but  not  its 
slope,  a  constant  will  be  added  to  the  value  of  the  discrimination  variable  for 
every  point  on  the  plot.  This  additive  constant  does  not  affect  the  separation 
between  the  earthquake  and  explosion  populations,  and  hence  it  may  be  chosen 
arbitrarily.  In  Figure  15  and  16  this  constant  has  been  chosen  in  such  a  manner 
as  to  utilize  as  fully  as  possible  those  magnitudes  which  were  approximated 
by  upper  bounds.  Specifically,  the  intercepts  of  the  lines  were  visually  ad¬ 
justed  to  maximize  the  number  of  explosion  upper  bounds  and  minimize  the  number 
of  earthquake  upper  bounds  lying  on  the  "explosion"  side  of  the  discriminant 
line.  The  reason  that  such  adjustment  is  advantageous  is  that  all  upper  bounds 
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Fipure  lba.  Low-frequency  P-wave  magnitude  P,. 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P„. 
Earthquakes  are  denoted  by  crosses,  single  explosions  by 
circles,  and  multiple  explosions  by  triangles. 
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Figure  16c.  Low-frequency  long-period  P-wave  magnitude  LP^ 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  Py 
Earthquakes  are  denoted  by  crosses,  single  explosions  y 
circles,  and  multiple  explosions  by  triangles. 
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MIDDLE-FREQ,  LONG-PERIOD  P  VS,  HIGH-FREO,  P 


Figure  L6d.  Miudle-f requency  long-period  P-wave  magnitude  LP^- 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P  . 
Earthquakes  are  denoted  by  crosses,  single  explosions  By 
circles,  and  multiple  explosions  by  triangles. 
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Figure  16e.  High-frequency  long-period  P-wave  magnitude  LP  . 

Spectral  magnitudes  of  presumed  earthquakes  ana  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P  , 
Earthquakes  are  denoted  by  crosses,  single  explosions  5y 
circles,  and  multiple  explosions  by  triangles. 
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Figure  16f.  Low-frequency  long-period  S-wave  magnitude  LS^. 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P  . 

Earthquakes  are  denoted  by  crosses,  single  explosions  By 
circles,  and  multiple  explosions  by  triangles. 
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Figure  16g.  Middle-frequency  long-period  S-wave  magnitude  LS^, 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P„. 
Earthquakes  are  denoted  by  crosses,  single  explosions  by 
circles,  and  multiple  explosions  by  triangles. 
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Figure  16h.  High-frequency  long-period  S-wave  magnitude  LS  . 

Spectral  magnitudes  of  presumed  earthquakes  ana  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P  • 

Earthquakes  are  denoted  by  crosses,  single  explosions  By 
circles,  and  multiple  explosions  by  triangles. 
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Figure  16i. 
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.  Low-frequency  Rayleigh-wave  magnitude  LR^ . 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P_. 
Earthquakes  are  denoted  by  crosses,  single  explosions  By 
circles,  and  multiple  explosions  by  triangles. 
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16 j •  Middle-frequency  Rayleigh-wave  magnitude  LR^. 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P  . 
Earthquakes  are  denoted  by  crosses,  single  explosions  by 
circles,  and  multiple  explosions  by  triangles. 
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Figure  16k.  High-frequency  Rayleigh-wave  magnitude  LR^. 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P  . 
Earthquakes  are  denoted  by  crosses,  single  explosions  By 
circles,  and  multiple  explosions  by  triangles. 
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Figure  161.  Low-frequency  Love-wave  magnitude  LQ  . 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-f requency  P-wave  magnitude  P  . 
Earthquakes  are  denoted  by  crosses,  single  explosions  by 
circles,  and  multiple  explosions  by  triangles. 
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Figure  16m.  Middle-frequency  Love-wave  magnitude  LQ^. 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P 
Earthquakes  are  denoted  by  crosses,  single  explosions  By 
circles,  and  multiple  explosions  by  triangles. 
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Low-frequency  Lg-wave  magnitude  Lg  . 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P  . 
Earthquakes  are  denoted  by  crosses,  single  explosions  by 
circles,  and  multiple  explosions  by  triangles. 
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Figure  16p.  Middle-frequency  Lg-wave  magnitude  Lg^. 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P  . 
Earthquakes  are  denoted  by  crosses,  single  explosions  by 
circles,  and  multiple  explosions  by  triangles. 
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HIGH-FRE3.  LG  VS.  HIGH-FRE3.  P 

Figure  16q.  High-frequency  Lg-wave  magnitude  Lg  . 

Spectral  magnitudes  of  presumed  earthquakes  and  explosions 
as  a  function  of  the  high-frequency  P-wave  magnitude  P,. 
Earthquakes  are  denoted  by  crosses,  single  explosions  By 
circles,  and  multiple  explosions  by  triangles. 
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lying  below  the  discriminant  line  may  be  presumed  to  lie  within  the  explosion 
population,  even  though  the  true  positions  of  these  points  are  unknown.  Phy¬ 
sically,  this  geometrical  argument  is  equivalent  to  identifying  as  explosions 
all  those  events  for  which  the  spectral  ratio  defined  by 

spectral  ratio  =  amplitude  variable^/ (amplitude  variable^)01  (31) 

is  less  than  some  specified  threshold  which  explosions  are  presumed  to  exceed. 
Adjusting  the  intercept  of  the  discriminant  dine  to  reflect  this  threshold  and 
then  retaining  as  useful  data  those  upper  bounds  which  fall  below  the  line  and 
rejecting  those  which  lie  above  it  is  an  application  of  the  principle  of  neg¬ 
ative  discrimination.  Moving  the  discriminant  lines  parallel  to  themselves  in 
Figures  15  and  16  would  have  the  effect  merely  of  changing  the  number  of  upper 
bounds  which  are  retained  as  useful  data  points.  If  for  some  reason  it  was 
t>el  “Ved  that  negative  discrimination  should  not  be  applied,  the  intercepts  of 
the  discriminant  lines  could  be  chosen  to  be  so  low  that  only  actual  measure¬ 
ments,  and  no  upper  bounds,  are  retained  in  the  data  base. 

Table  XIII  lists  the  discriminant  lines  which  are  shown  in  Figures  15  ana 
16.  The  lines  were  fit  visually,  and  their  exact  values  were  chosen  to  be  con¬ 
sistent  with  those  of  a  companion  study  employing  a  somewhat  larger  data  base 
(Rivers  et  al,  1979).  It  can  be  seen  that  the  slope  of  the  short-period  P- 
wave  spectral  ratio  increases  with  increasing  separation  of  the  ratioed  fre¬ 
quency  bands  and  then  stays  constant  as  the  separation  is  increased  from  P  /P 

13 

to  LP^/Py  This  same  constant  slope  also  applies  to  long-period  S.  This  con¬ 
sistency  perhaps  indicates  that  the  determination  of  the  coefficients  a,  al¬ 
though  empirical,  does  in  fact  measure  some  real  physical  effect.  The  slopes 
of  the  discriminant  lines  for  surface  waves  differ  from  the  unit  slope  of  the 
classical  discriminant,  but  this  difference  was  anticipated  since  the  P- 

wave  magnitude  P^  is  measured  at  a  higher  frequency  than  is  n^.  A  steeper 
slope  should  result  for  the  discriminant  if  m^  is  measured  on  the  de¬ 

caying  asymptote  of  the  P-wave  spectrum  rather  than  near  the  corner  frequency. 


Rivers,  D.  W.  ,  D.  H.  von  Seggem,  B.  L.  Elkins,  P.  J.  Klouda,  J.  A.  Burnetti, 
and  I.  Megyesi  (1979).  (S)  A  statistical  discrimination  experiment  for 

Eurasian  events  using  the  Priorities  I  and  II  networks  (1!) ,  SDAC  Report 
No.  fR-79-2,  Teledyne  Geotech,  Alexandria,  Virginia. 
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TABLE  XIII 


Discriminant  Lines  Fit  to  Figures  15-17 


P  =  1.41  P2-1.03 
P2  =  1.72  P  -1.45 
P  =  1.77  P3-0.87 


LP1 

= 

1.77 

P3-0.14 

LP2 

= 

1.77 

P3+0.02 

“3 

= 

1.77 

P3+1.22 

LS1 

= 

1.77 

P  3+0 . 2 1 

ls2 

= 

1.77 

P3+0.50 

LS3 

= 

1.77 

P  +1.22 

LR1 

= 

2.31 

P3=2.79 

LR2 

= 

2.31 

P3-2.94 

lr3 

= 

2.27 

P3-2.55 

LQ1 

= 

2.30 

P3-3.21 

LQ2 

= 

2.30 

P3-3.34 

lq3 

= 

2.30 

P3-2.39 

o„(- 

-2) 

=  -3fc(-2)+5.93 

fl0(-3)  =  -3fc(-3)+5. 97 
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It  can  be  seen  from  Table  XIII  that  the  slopes  for  all  six  discriminants  of  the 
form  M  :m^  were  not  quite  equal.  The  differences  are  probably  only  a  stati¬ 
stical  fluctuation,  and  a  slope  of  2.30  could  probably  have  been  imposed  upon 
each  of  the  six  discriminants  without  significantly  changing  the  separation 
between  the  two  populations. 

No  discriminants  involving  L  were  computed  since,  as  is  shown  in  Figures 
16  o-q,  insufficient  explosion  data  existed  to  permit  such  discriminants  to  be 
found.  This  deficiency  is  due  partly  to  nondetection  but  mainly  to  the  omis¬ 
sion  of  the  L  window  from  the  archived  data  set,  especially  at  stations  such 
g 

as  KAAO  and  MAIO  which  are  at  regional  distances  to  many  events.  The  failure 

to  construct  a  discriminant  utilizing  L  ,  at  certain  stations  if  not  averaged 

g 

over  an  entire  network,  and  to  compare  its  effectiveness  with  that  of  the 
teleseismic  discriminants  must  be  regarded  as  one  of  the  major  disappointments 
of  this  experiment. 

The  discrimination  value  of  the  spectral-fitting  algorithm  which  was  il¬ 
lustrated  in  Figure  2  lies  in  the  determination  of  low-frequency  spectral  level 
as  a  function  of  corner  frequency.  For  a  given  level,  explosions  should  have 
a  higher  corner  frequency  than  do  earthquakes  because  their  source  time  dimen¬ 
sions  are  smaller  (Hanks  and  Thatcher,  1972).  Figure  17  at  least  suggests  that 

this  effect  can  be  measured  from  the  spectral  fitting  algorithm,  using  either 

-2-3 

the  f  or  the  f  decay  models,  but  it  hardly  demonstrates  that  it  contains 
much  discrimination  power.  The  lines  drawn  between  the  two  populations  in 
Figures  17a  and  17b  have  slope  -3,  which  is  the  theoretical  value  for  events 
having  a  given  stress  drop.  The  intercepts  of  these  lines  are  irrelevant, 
since  no  upper  bounds  are  applicable  to  noise  measurements  of  this  pair  of 
parameters. 


Hanks,  T.,  and  W.  Thatcher  (1972).  A  graphical  representation  of  seismic 
source  parameters,  J.  Geophys.  Res.,  77 ,  4393. 
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Figure  17a.  f  spectral  fit.  ,  , 

Low-frequency  spectral  level  of  presumed  earthquakes  and 

explosions  as  a  function  of  corner  frequency.  Earthquakes 
are  denoted  by  crosses,  single  explosions  by  circles,  and 
multiple  explosions  by  triangles. 
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Figure  17b.  f  spectral  fit. 

Low-frequency  spectral  level  of  presumed  earthquakes  and 
explosions  as  a  function  of  corner  frequency.  Earthquakes 
are  denoted  by  crosses,  single  explosions  by  circles,  and 
multiple  explosions  by  triangles. 


The  failure  of  the  two  populations  to  separate  clearly  in  Figure  17  is 
probably  due  to  the  distortion  of  the  source  spectra  by  signal  propagation 
effects.  As  was  mentioned  in  the  description  of  the  spectral-fitting  algorithm, 
the  values  chosen  for  and  f£  are  critically  dependent  on  the  value  chosen 
for  t*,  so  at  least  part  of  the  weakness  of  this  discriminant  may  be  attribut¬ 
able  to  an  erroneous  attenuation  correction.  It  may  be  that  this  discriminant 
is  too  strongly  influenced  by  errors  in  the  values  which  were  calculated  for  t* 
for  each  epicenter-to-receiver  path  to  permit  meaningful  network  averages  to  be 
formed.  The  power  of  the  discriminant  might  thus  be  enhanced  if  the  propagation 
paths  were  made  constant  by  applying  the  versus  f£  discriminant  only  to  data 
from  events  occurring  at  the  same  depth  within  a  single  tectonic  region  as 
recorded  at  a  single  station.  Another  way  of  improving  the  discriminant's  power 
would  be  to  use  network- averaged  values  of  and  f^  as  was  done  in  Figure  17, 
but  to  use  better  values  of  t*  in  the  calculation  of  the  spectral  fits.  It  will 
be  recalled  that  a  value  of  t*  was  calculated  for  each  individual  spectral  fit, 
and  so  these  values  were  subject  to  a  significant  statistical  fluctuation.  This 
fluctuation  was  particularly  large  for  those  cases  in  which  the  data  from  only 
a  narrow  spectral  band  was  of  a  sufficient  signal-to-noise  ratio  to  be  used.  A 
more  reliable  attenuation  correction  might  thus  be  made  by  using  not  the  value 
of  t*  for  a  given  individual  case  but  rather  one  which  is  determined  by  averag¬ 
ing  the  individual  values  found  for  all  events  occurring  within  that  same  tec¬ 
tonic  region  as  measured  at  that  same  station  (and  using  the  same  theoretical 
model,  either  f-2  or  f-3) .  The  scatter  of  the  individual  values  about  these 
averages  is  shown  in  Appendix  III.  The  scatter  in  the  values  which  were  used 
for  t*  for  the  same  path  for  different  events  is  clearly  too  large,  so  more 
reliable  (or  at  least  more  self-consistent)  values  would  have  been  found  for 
and  f  had  the  same  value  been  used  for  t*  every  time.  Another  suggestion 
for  improvement  of  this  discriminant  was  made  earlier  in  this  report,  namely 
that  spectra  be  used  only  from  those  stations  for  which  30°<A<90°,  in  order  to 
avoid  uncertainties  introduced  into  the  calculation  of  the  magnitude  of  the 
low-frequency  level  by  the  correction  for  geometrical  spreading  of  the  wavefront 
It  should  also  be  noted  that  this  discriminant  is  probably  strongly  affected  by 
the  assumption  of  zero  depth  for  all  events.  In  light  of  these  errors  intro¬ 
duced  into  this  discriminant  by  propagation  effects,  it  is  hardly  surprising 
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that  Figure  17  shows  such  large  scatter.  For  the  purposes  of  this  experiment, 
we  merely  note  that  this  scatter  exists,  and  then  we  calculate  the  versus  f 
discrimination  variable  using  a  slope  of  -3. 

The  final  discrimination  variables  are  the  (logarithmic)  P-wave  comp¬ 
lexities.  The  distributions  of  the  complexities  for  the  earthquake  and  the 
explosion  populations  are  shown  in  Figure  18.  The  separation  between  the  two 
populations  is  seen  to  increase  with  increasing  length  of  the  time  window  during 
which  the  complexity  was  measured.  The  two  populations  have  a  large  overlap, 
and  certain  events  are  conspicuous  outliers.  Some  of  these  outliers  are  due 
to  insufficient  data,  being  based  upon  only  one  or  two  observations  (cf.  Appendix 
II),  but  others  reflect  real  source  or  path  effects.  In  particular,  deep 
earthquakes  are  characterized  by  small  complexities  since  their  waveforms  pass 
through  the  crust  once  rather  than  twice.  Two  of  the  multiple  explosions, 
events  265  and  271,  show  large  complexities  on  account  of  the  large  amount  of 
energy  arriving  from  the  second  explosion  during  the  coda  window  of  the  first. 
Although  the  other  multiple  explosion,  event  53,  should  also  be  expected  to 
have  an  anomalously  large  complexity  for  the  same  reason,  it  does  not  do  so. 

A  possible  explanation  for  the  explosion-like  complexity  of  event  53  is  that 
some  of  the  energy  from  the  second  explosion  arrives  within  the  signal  window 
of  the  first,  making  the  "signal"  of  the  hypothetical  single  event  look  large 
relative  to  its  "coda".  The  outlying  events  in  Figure  18  may  be  particularly 
prone  to  misclassif ication,  and  special  note  will  be  made  of  them  in  the  analysis 
of  the  discrimination  results. 

The  twenty  variables  which  were  used  as  the  basis  for  the  multivariate 
discrimination  are  listed  in  Table  XIV,  and  their  values  for  each  event  are 
given  in  Appendix  IV.  Also  in  Table  XIV  are  listed  the  number  of  events  to 
which  each  discriminant  was  inapplicable,  either  because  there  was  no  measure¬ 
ment  at  any  station  of  one  of  the  magnitude  variables  used  in  the  discriminant 
or  because  the  magnitudes  of  both  variables  in  the  discriminant  were  upper 
bounds  estimated  from  noise  measurements  only.  Next  are  listed  the  number  of 
events  for  which  one  of  the  magnitude  variables  was  an  upper  bound  and  for 
which  the  discriminant  variable  was  deleted  from  the  data  set  because  the  true 
value  could  lie  within  either  of  the  two  populations.  These  rejected  events 
are  those  which  are  shown  in  Figures  15  and  16  as  points  with  arrows  attached 
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Figure  18a.  Short  time  window. 

Complexity  for  presumed  earthquakes  and  explosions.  Certain 
outliers  are  identified. 


EXPLOSIONS 


TABLE  XIV 


Variables  Employed  in  Discrimination 


A)  Number 

of 

events 

to  which  the  discriminant  is  inapplicable  due 

to  insufficient  data. 

B)  Number 

of 

events 

for 

which  an  upper  bound  to  the  variable  was 

computed  but  then  rejected 

because 

;  the  true 

value  could  lie  on  either  side  of  the  discriminant  line. 

C)  Number 

of 

events 

for 

which  the  discrimination  variable  was  app 

roximated  by  an 

upper 

bound. 

D)  Number 

of 

events 

for 

which  the  true  value  of  the  discrimination  variable  was 

calculated . 

E)  Number 

of 

events 

for 

which  the  discrimination  variable  is  non- 

zero.  (Sura  of 

(C)  and  (D) 

.) 

Variable 

Number 

Discription  A  B 

C 

D 

E 

1 

P  -1.A1  P2  5  13 

2 

112 

11A 

2 

P2-1.72  P3  35 

7 

117 

12A 

3 

P  -1.77  P3  6  11 

6 

109 

115 

A 

ft0(-2)+3  f c (-2)  16 

— 

116 

116 

5 

Q„(-3)+3  f c (-3)  16 

— 

116 

116 

6 

corap^  19 

— 

113 

113 

7 

corap 2  22 

— 

110 

110 

8 

corap 3  AO 

— 

92 

92 

9 

LPX-1.77  P3  12  67 

9 

AA 

53 

10 

LP2-1.77  P3  11  75 

13 

33 

A6 

11 

LP3-1.77  P3  12  71 

17 

32 

A9 

12 

LSj-1.77  P3  10  60 

8 

5A 

62 

13 

LS2-1.77  F3  12  63 

8 

A9 

57 

1A 

LS3-1.77  P3  12  68 

10 

A2 

52 

15 

LR—2.31  P3  7  33 

13 

79 

92 

16 

LR2-2.31  P3  5  29 

9 

89 

98 

17 

LR3-2.27  P3  7  38 

7 

80 

87 

18 

LQ1-2.30  P3  5  35 

12 

80 

92 

19 

LQ2-2.30  P3  6  33 

8 

85 

93 

20 

LQ3-2.30  P3  9  38 

9 

76 

85 
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to  them  which  point  towards  the  (dashed)  discriminant  line.  The  next  column 
of  Table  XIV  shows  the  number  of  events  for  which  the  upper  bound  estimate  was 
retained  in  the  data  base  as  an  approximation  to  the  true  value  of  the  dis¬ 
crimination  variable.  These  events,  to  which  the  principle  of  negative  dis¬ 
crimination  is  applicable,  are  those  which  are  shown  in  Figures  15  and  16  as 
points  with  arrows  attached  which  point  away  from  the  discriminant  line.  Next 
are  shown  the  number  of  events  for  which  the  discrimination  variable  was 
actually  measured,  and  finally  there  are  listed  the  total  number  of  events  to 
which  the  discriminant  is  applicable.  Table  XIV  shows  that  for  every  discrim¬ 
inant  except  those  involving  long-period  body  waves  more  than  85%  of  the  mag¬ 
nitudes  retained  in  the  data  base  were  "true"  values,  rather  than  upper  bounds 


MULTIVARIATE  DISCRIMINATION 


Discrimination  can  be  performed  using  as  many  of  the  nineteen  variables 
as  desired;  instead  of  a  line  separating  the  earthquake  and  explosion  pop¬ 
ulations  as  with  the  classical  Mg  versus  m^  discriminant,  the  two  populations 
are  separated  by  an  N-l  dimensional  hyperplane  in  the  space  spanned  by  the  N 
chosen  variables.  A  practical  limit  is  imposed  upon  the  number  of  variables, 
however,  by  the  fact  that  as  more  variables  are  used  in  a  single  discriminant, 
there  are  fewer  events  in  the  data  base  for  which  all  the  chosen  variables  are 
non— zero.  The  N  coefficients  determining  the  hyperplane  are  referred  to  as 
the  "discriminant  function". 

The  discrimination  was  performed  using  the  stepwise  discriminant  analysis 
program  BMD07M  developed  by  Jennrich  (1977).  The  first  step  of  the  discrimin¬ 
ation  consists  of  the  a  priori  classification  of  a  suite  of  events  as  earth¬ 
quakes  or  explosions,  as  was  done  in  the  previous  section  of  this  report.  These 
events  are  used  as  a  "training  set"  to  calculate  the  discriminant  function  which 
best  separates  the  two  populations,  and  this  function  is  then  applied  to  the 
remaining  events.  As  a  check  against  misclassif ied  events'  being  used  in  the 
training  set,  the  discriminant  function  is  used  to  calculate  a  posteriori 
probabilities  of  correct  classification  for  all  events  of  the  training  set 
before  the  discriminant  is  applied  to  unknowns.  Unknown  events  are  classified 
as  explosions  or  earthquakes  depending  upon  on  which  side  of  the  hyperplane  they 
lie,  and  probabilities  of  misclassif ication  are  given  by  their  proximity  to 
the  other  population. 

Output  of  the  multivariate  discrimination  program  consists  not  only  of  the 
discriminant  function  and  the  a_  posteriori  probabilities  but  also  of  a  table 
which  ranks  the  N  chosen  variables  in  the  order  of  their  importance  in  deter¬ 
mining  the  equation  of  the  discriminant  hyperplane.  This  ranking  is  calculated 
by  a  stepwise  process  that  computes  which  of  the  N-m  remaining  variables  will 
yield  the  most  improvement  in  discrimination  of  the  training  set  when  added 
to  the  m  variables  already  used  in  the  previous  steps.  The  stepwise  process 
re-evaluates  the  weight  assigned  to  each  variable  whenever  a  new  one  is  added, 

Jennrich,  R.  I.  (1977).  Stepwise  discriminant  analysis,  in  Statistical  Methods 
for  Digital  Computers,  edited  by  K.  Enslein,  A.  Ralston,  and  H.  S.  Wiff, 

John  Wiley  and  Sons,  New  York. 
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and  it  may  then  delete  one  of  the  previously  used  variables  if  it  is  deemed 
to  yield  superfluous  or  contradictory  information  to  the  discrimination.  This 
ranking  of  variables  is  described  in  detail  by  von  Seggern  and  Rivers  (1979). 

All  132  events,  with  the  exception  of  the  three  multiple  explosions,  are 
to  be  used  in  the  training  set,  so  the  discrimination  experiment  will  essential¬ 
ly  be  a  test  of  whether  each  event's  a  posteriori  classification  agrees  with 
the  classification  which  was  assigned  to  it  _a  priori.  Although  the  a.  posteriori 
classification  in  effect  regards  the  training  set  events  as  unknowns,  the  only 
true  unknowns  (in  the  sense  of  not  having  been  used  to  calculate  the  discrimin¬ 
ant  functions)  are  the  three  shot  arrays.  The  question  then  arises:  how 
effective  will  the  discriminant  functions  be  when  they  are  applied  not  to  the 
129  events  of  the  training  set  but  rather  to  the  next  group  of,  say,  fifty 
events  which  are  true  unknowns?  The  a.  posteriori  discrimination  provides  an 
approximate  answer  to  that  question,  in  that  the  classification  probability  of 
each  event  computed  on  the  basis  of  the  129-event  training  set  is  a  good  approxi¬ 
mation  to  that  probability  which  would  result  from  the  use  of  only  the  other 
128  events  in  the  training  set.  Classification  of  the  next  fifty  unknowns 
should  thus  work  approximately  as  well  as  that  of  the  training  set.  In  practice, 
seismic  discrimination  ought  to  be  a  learning  process,  and  the  fifty  newly- 
classified  events  should  be  added  to  the  training  set  so  that  the  next  group 
of  events  can  be  classified  using  a  179-event  training  set.  The  influence 
exerted  on  the  discrimination  by  the  size  of  the  training  set  will  be  measured 
in  a  subsequent  section  of  this  report. 


DISCRIMINATION  RESULTS 


Using  the  number  of  events  to  which  each  discrimination  variable  is 
applicable,  as  listed  in  the  final  column  of  Table  XIV,  a  systematic  approach 
may  be  devised  to  the  problem  of  collecting  various  subsets  of  the  twenty 
variables  to  form  multivariate  discriminants.  Neglecting  for  the  moment  the 
three  multiple  explosions,  to  which  all  of  the  twenty  variables  are  applicable, 
Table  XIV  shows  that  the  least  frequently  observed  variable,  the  spectral 
ratio  of  the  middle-frequency  long-period  P-wave  magnitude  LP  to  the  high- 
frequency  short-period  P-wave  magnitude  P  ,  is  applicable  to  43  of  the  129 
events  to  which  an  a^  priori  classification  has  been  assigned.  A  twenty- 
variable  discriminant  can  thus  be  applied  to  at  most  43  events  in  the  train¬ 
ing  set.  It  turns  out,  however,  that  at  least  one  of  the  nineteen  other 
variables  is  inapplicable  to  19  of  these  43  events,  so  the  twenty-variable 
discriminant  can  actually  be  applied  to  only  24  events.  The  next  least  frequ¬ 
ently  observed  variable  is  applicable  to  46  events,  but  at  least  one  of  the 
remaining  eighteen  variables  is  inapplicable  to  eighteen  of  them,  so  the 
discriminant  formed  from  the  nineteen  most  frequently  observed  variables  can 
be  applied  to  only  28  events.  Proceeding  in  this  manner,  at  each  step  delet¬ 
ing  from  the  subset  of  variables  used  in  the  multivariate  discriminant  the 
next  least  frequently  observed  variable,  it  is  possible  to  construct  the 
systematic  incrementation  of  variables  which  is  shown  in  Table  XV.  No  two- 
variable  discriminant  is  listed  in  the  table,  since  it  was  not  desired  to 

create  a  multivariate  discriminant  which  included  the  ft0(-2)  versus  f  (-2) 

c 

variable  but  not  the  0o(-3)  versus  f  ( — 2)  variable,  or  vice  versa.  The  first 
attempt  at  an  ja  posteriori  classification  of  the  132  events  consisted  of 
applying  each  of  the  nineteen  discriminants  listed  in  Table  XV  to  the  train¬ 
ing  set  and  seeing  how  many  events  each  discriminant  classified  "correctly" 
(i.e.,  in  agreement  with  the  a  priori  classification). 

Results  of  the  a.  posteriori  classification  of  the  training  set  (and  also 
of  the  three  multiple  explosions,  events  53,  265,  and  271)  are  shown  in 
Table  XVI.  Listed  alongside  each  event  in  that  table  is  its  a_  priori  classi¬ 
fication:  "S"  for  "shot",  "Q"  for  "quake",  or  "M"  for  "multiple".  Next  are 
listed  the  a.  posteriori  classifications  generated  by  each  of  the  nineteen 
multivariate  discriminants  shown  in  Table  XV.  Events  are  classified  as  being 
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FESULTS  OF  TH3  DISCRIMINATION  USING  THE  SYSTEMATIC  INCREMENTATION  OF  VARIABLES 
REFFS  TO  TABLE  XTII  FOE  A  DESCRIPTION  OF  THE  MULTIVARIATE  DISCRIMINANTS 
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explosions  or  earthquakes  depending  upon  whether  the  computed  probability  of 
being  an  explosion  is  greater  than  or  less  than  one-half.  Asterisks  appear 
beside  those  classifications  for  which  this  probability  lies  between  twenty 
and  eighty  percent.  Denoting  the  probability  of  being  an  explosion  by  P(EX), 
the  classification  system  used  in  Table  XVI  may  be  written  as: 

S  if  0.8  <P (EX)  <1.0 

S*  if  0.5  <P (EX)  <0.8  (32) 

Classification  Q*  if  0.2  <P(EX)  <0.5 

Q  if  0.0  <P (EX)  <0.2 

For  each  multivariate  discriminant  Table  XV  shows  the  number  of  instances 
in  which  the  a  priori  and  ja  posteriori  classifications  were  in  agreement  or 
in  disagreement.  It  is  shown  that,  as  one  would  expect,  the  applicability  of 
the  discriminants  increases  but  their  reliability  decreases  as  the  number  of 
variables  decreases  from  twenty  to  one.  In  fact,  the  one— variable  discriminant 
P2  -1.72  P^  misclassif ies  more  explosions  than  it  classifies  correctly,  al¬ 
though  it  does  correctly  classify  some  95%  of  the  earthquakes  to  which  it 
is  applicable.  This  asymmetry  is  brought  about  by  the  difference  in  the  sizes 
of  the  two  populations  in  the  training  set  and  by  the  discrimination  algorithm's 
treatment  of  both  types  of  misclassif ication  as  being  equally  undesirable. 

It  will  be  shown  later  that  a  better  classification  percentage  may  be  obtained 

for  the  explosions  (at  the  expense  of  that  of  the  earthquakes)  by  weighting 

the  two  populations  unequally. 

Thus  far  121  of  the  132  events  have  been  examined  using  as  many  as  nine¬ 
teen,  or  as  few  as  one,  different  multivariate  discriminants.  Considering 
now  the  121  events  to  be  unknowns,  it  must  be  determined  which  one  of  these 
discriminants  should  be  used  to  assign  the  best  possible  classification  to 
each  event.  In  general  the  "best"  such  discriminant  would  be  expected  to  be 
the  one  which  utilizes  the  most  information  about  each  unknown  event,  i.e., 

the  one  containing  the  most  variables.  Table  XV  shows,  however,  that  the 

discriminants  which  employ  the  most  variables  are  poorly  calculated  since  the 
number  of  events  in  the  training  set  which  was  used  to  calculate  them  does 
not  greatly  exceed  the  number  of  variables.  Such  discriminants  should  be  re¬ 
garded  as  too  unreliable  to  use  for  classifying  unknowns  since  their  exact 
form  is  too  strongly  dependent  upon  the  selection  of  events  comprising  the 
training  set.  A  subjective  evaluation  must  be  made  of  the  trade-off  between 
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the  amount  of  information  used  by  each  discriminant,  as  measured  by  the  number 
of  variables,  and  the  discriminant's  reliability,  as  measured  by  the  number 
of  events  of  both  populations  in  the  training  set.  For  the  purposes  of  this 
report,  it  was  decided  to  reject  as  being  unreliable  those  discriminants  lisLed 
in  Table  XV  which  contain  fifteen  or  more  variables,  since  they  were  cal¬ 
culated  using  38  or  fewer  data  points.  The  discriminants  containing  fourteen 
or  fewer  variables  were  calculated  using  51  or  more  data  points,  however,  and 
it  was  somewhat  arbitrarily  decided  that  they  should  be  regarded  as  reliable. 
The  rule  which  was  used  for  classification,  then,  was  to  classify  each  of  the 
121  events  using  the  discriminant  with  the  greatest  possible  number  of  vari¬ 
ables,  up  to  a  maximum  of  fourteen.  The  results  of  applying  the  criterion  for 
classification  are  shown  in  Table  XVII.  Refering  to  the  description  of  the 
discriminants  as  given  in  the  second  column  of  Table  XV,  it  can  be  seen  that 
no  discriminant  containing  fourteen  or  fewer  variables  uses  any  of  the  long- 
period  body-wave  variables,  numbers  9-14  in  Table  XIV.  The  discrimination 
results  which  are  presented  herein  are  therefore  unaffected  by  all  the  measure¬ 
ments  which  were  made  of  long-period  P-and  S-waves.  Figures  16  c-h  neverthe¬ 
less  show  that  these  deleted  variables  are  effective  discriminants,  and  their 
worth  should  not  be  ignored  in  future  studies. 

The  systematic  incrementation  of  variables  which  was  described  in  Table 
XV  and  applied  in  Table  XVI  is  only  one  of  a  multitude  of  methods  for  creating 
the  multivariate  discriminants.  For  example,  there  are  (jq)  =  184,756  dif¬ 
ferent  ten-variable  discriminants  alone.  It  is  entirely  possible,  if  not  in¬ 
deed  probable,  that  one  of  the  184,755  other  combinations  would  work  better 
than  tl 2  particular  one  which  resulted  from  the  incrementation  shown  in  Table 
XV.  In  particular,  many  of  the  other  ten-variable  discriminants  might  be 
applicable  to  certain  events  to  which  the  chosen  ten-variable  discriminant 
could  not  be  applied.  Thus  in  Table  XVII  many  events  were  classified  using 
fewer  variables  than  had  been  determined  for  them.  For  example,  because  dis¬ 
crimination  variable  number  6  was  inapplicable  to  event  26,  that  event  was 
classified  using  five  variables,  although  it  is  shown  in  Appendix  IV  that  ten 
of  the  fourteen  variables  used  for  discrimination  could  have  been  applied  to 
it.  Worse,  nine  events  (designated  "unk."  in  Table  XVII)  were  not  classified 
at  all  by  the  systematic  incrementation  of  variables  because  discrimination 
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variable  number  2,  which  was  used  in  all  the  multivariate  discriminants, 
was  not  applicable  to  them.  Special  multivariate  discriminants  were  there- 

created  and  applied  to  events  such  as  number  z.6  in  order  to  classify  them 
using  as  much  information  as  possible.  Table  XVIII  lists  these  special  discri¬ 
minants,  their  effectiveness  is  classifying  the  129  events  of  the  training  set, 
and  the  classifications  which  result  from  applying  them  to  the  specific  events 
for  which  they  were  created.  These  classifications  are  to  be  considered  more 
reliable  for  these  particular  events  than  the  classifications  given  for  them 
in  Table  XVII. 

The  most  reliable  classification  for  each  event  is  thus  given  either  in 
Table  XVIII  for  all  the  events  listed  there  or  in  Table  XVII  for  all  other 
events.  A  summary  of  these  most  reliable  classifications  is  presented  in 
Table  XIX.  Alongside  the  a_  priori  classification  which  was  assigned  to  each 
event  ("EX"  for  explosion,  "Q"  for  earthquake,  "M"  for  multiple  explosion) 
there  are  listed  the  number  of  variables  which  were  used  for  the  most  reliable 
a  posteriori  classification.  Next  are  listed  the  a  posteriori  classifications  if 
they  are  uncertain  or  if  they  arc  in  disagreement  with  the  a^  priori  assigna¬ 
tions.  It  can  be  seen  from  the  table  that  the  multivariate  discrimination 
employed  in  this  report  resulted  in  some  six  misclassif ications  (assuming,  of 
course,  that  the  a  priori  classifications  are  in  fact  correct).  The  percentage 
of  misclassif ications  should  not  be  regarded  as  a  definitive  measure  of  the 
reliability  of  the  discrimination  process,  however,  since  it  is  obviously 
dependent  upon  the  particular  suite  of  events  chosen  for  use  in  the  analysis 
and  since  no  mention  has  been  made  of  the  uncertainty  associated  with  the 
classification  of  each  event.  It  is  difficult  to  make  any  quantitative 
assessment  of  the  confidence  with  which  any  given  event  may  be  classified 
since  uncertainty  enters  into  the  discrimination  at  several  different  stages. 
Sources  of  the  uncertainty  are: 

1)  Measurement  of  the  27  parameters  listed  in  Table  VI.  These  measure¬ 
ments  are  affected  by  the  starting  point  and  length  of  the  data 
windows  in  the  time  domain,  amplitude  leakage  between  discrete 
frequencies  in  the  computation  of  the  spectra,  etc. 

2)  Calculation  of  the  magnitudes  from  the  27  measured  parameters.  Errors 
arise  in  these  calculations  due  to  errors  in  the  assumed  values  of 
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the  B-factors,  attenuation  coefficients,  distances,  etc. 

3)  Calculation  of  the  maximum-likelihood  estimates  of  the  27  magnitudes. 

As  a  rule  of  thumb,  the  most  reliable  estimates  are  those  which  are 

(a)  based  on  magnitudes  from  the  greatest  number  of  stations,  and 

(b)  based  on  the  most  signal  detections  and  the  fewest  noise  levels. 

A  rough  idea  of  the  relative  uncertainty  of  the  network  estimates  of 
each  magnitude  for  each  event  may  thus  be  gained  from  inspection  of 
Appendix  II.  That  this  rule  of  thumb  is  only  an  approximate  guide¬ 
line  may  be  seen,  however,  by  considering  two  events  for  which  all 
the  measurements  of  one  of  the  parameters  were  based  on  signal  detec¬ 
tions  rather  than  noise  levels.  Although  the  standard  deviation  of 
the  mean  for  the  magnitude  of  that  parameter  is  likely  to  be  smaller 
for  the  event  which  was  observed  at  more  stations  than  was  the  other, 
there  is  certainly  no  guarantee  that  this  is  the  case.  The  less 
frequently  observed  event  may  in  fact  be  better  determined,  and  the 
only  sure  comparison  of  the  relative  uncertainty  in  the  two  event 
magnitudes  is  given  by  the  standard  deviation  of  the  means  for  both. 

It  should  be  noted  that  if  noise  levels  as  well  as  signal  detections 
were  used  to  compute  a  given  maximum-likelihood  estimate,  then  the 
standard  deviation  of  the  mean  cannot  be  used  but  an  analogous  con¬ 
fidence  limit  can  be  given  by  the  Cramer— Rao  bound  (Ringdal,  1976). 

4)  Determination  of  the  discrimination  variables.  Certain  of  the  co¬ 
efficients  oi  in  equation  (31)  are  more  poorly  determined  than  others, 
and  so  some  of  the  twenty  (later  decreased  to  fourteen)  discrimination 
variables  result  in  more  misclassif ications  than  do  others.  Although 
an  estimate  can  be  made  of  the  relative  reliability  of  each  discri¬ 
mination  variable,  it  is  not  necessary  to  do  so,  since  the  reliability 
of  the  multivariate  discriminant  functions  (considered  next)  can  be 
determined  independently  of  that  of  the  individual  variables  of  which 
they  are  composed. 

5)  Formation  of  the  multivariate  discrimination  functions.  It  has  already 
been  explained  that  the  reliability  of  the  discriminants  is  a  trade¬ 
off  between  the  number  of  events  in  the  training  set  used  to  deter¬ 
mine  them  and  the  fraction  of  those  events  which  are  misclassif ied 
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by  them.  This  information  is  presented  in  Tables  XV  and  XVIII. 

It  can  thus  be  seen  that  the  classification  of  events  such  as  number 
194,  which  Table  XVII  shows  to  have  been  classified  on  the  basis  of 
a  one-variable  discriminant,  should  be  regarded  as  highly  uncertain, 
since  Table  XV  indicates  that  this  discriminant  has  a  high  misclassi- 
fication  rate. 

6)  Computation  of  the  classification  probabilities.  Even  if  the  locations 
in  the  multivariate  space  of  the  discriminant  hyperplanes  (considered 
above)  were  taken  to  be  perfectly  determined,  uncertainty  wouJd  arise 
in  the  classification  of  a  given  event  due  to  the  proximity  of  that 
event's  data  poir t  to  the  particular  hyperplane  used  for  discrimination. 
Data  points  which  lie  distant  from  the  hyperplane  lie  well  within  the 
boundaries  of  either  the  earthquake  or  explosion  populations  and  hence 
have  a  low  probability  of  having  been  misclassif ied,  but  data  points 
near  the  hyperplane  are  within  a  region  of  some  overlap  of  the  two 
populations,  and  hence  their  classification  is  uncertain.  A  measure 
of  this  uncertainty  is  given  by  the  designations  which  are  explained 
in  equation  (32)  and  applied  to  Tables  XVII  and  XVIII.  Table  XIX 
shows  that  there  were  twelve  events  which  were  correctly  classified 
using  computed  probabilites  of  being  an  explosion  P(EX)  which  lie 
between  0.2  and  0.8. 

Uncertainty  in  the  discrimination  process  is  obviously  cumulative,  so  ;he 
factors  listed  above  (or  at  least  numbers  3,  5,  and  6,  which  are  highly  vari¬ 
able  from  event  to  event)  should  be  considered  in  evaluating  the  uncertainty 
associated  with  the  classification  of  any  given  event. 

Even  though  the  number  of  misclassif ied  events  cannot  be  used  by  itself 
to  evaluate  the  effectiveness  of  the  discrimination  process,  it  is  instructive 
to  examine  on  a  case-by-case  basis  the  six  misclassif ied  events  in  order  to 
determine  the  weaknesses  of  the  method  and  to  determine  whether  any  particular 
type  of  event  exists  which  is  more  prone  to  misclassif ication  than  others. 

These  six  events  are: 

1)  Event  8.  This  event  was  classified  on  the  basis  of  six  discrimination 
variables,  two  of  which  were  P-wave  signal  complexity  windows.  The 


complexity  was  based  on  a  measurement  at  only  one  station,  KSRS, 
and  hence  may  be  unreliable.  It  should  be  noted  that  at  KSRS  the 
P-wave  spectral  amplitudes  P  and  were  flagged  as  falling  below 
the  threshold  SNR  of  1.5;  since  the  signal  at  that  station  was  so 
noisy,  it  is  very  likely  that  the  complexity  of  the  coda  was  also 
dominated  by  noise.  When  a  four-variable  discriminant  not  using 
complexity  was  applied  to  this  event,  it  was  classified  as  an  earth¬ 
quake.  This  event  was  very  poorly  observed,  since,  as  is  shown  in 
Appendix  II,  the  single  measurement  at  KSRS  was  the  only  P-wave 
observation  in  the  entire  data  base.  (However,  there  did  exist  data 
windows  for  short-period  P  at  certain  other  stations  which  could  not 
be  measured  due  to  signal  processing  problems.)  Appendix  IV  shows 
that,  as  a  result,  the  P  :P2  and  P^P.^  spectral  ratios  could  not  be 
used  for  discrimination  and  that  the  P2:P3  spectral  ratio  which  was 
used  was  actually  an  upper  bound.  In  Figure  16b  event  8  is  indicated 
by  one  of  the  crosses  which  lie  on  the  left  of  the  discriminant  line 
and  which  have  arrows  attached  to  them  pointing  toward  the  left.  If 
the  magnitude  of  P^  had  been  represented  by  an  actual  measurement 
rather  than  by  an  upper  bound,  this  data  point  would  have  moved  farther 
into  the  earthquake  population.  The  misclassification  of  event  8  is 
thus  due  mainly  to  an  insufficiently  (and  possibly  incorrectly)  de¬ 
termined  value  for  the  complexity,  and  the  discrimination  was  influ¬ 
enced  heavily  by  the  unreliable  complexity  value  since  there  were  too 
few  observations  of  short— period  P— waves  to  compute  accurately  the 
spectral  ratios  P^:P2,  an<^  ^2:^3’ 

2)  Event  69.  The  NEIS  depth  for  this  event  is  255  km.  It  was  classified 
on  the  basis  of  the  two-variable  discriminant  consisting  of  the  spec¬ 
tral  ratios  LQ  :P3  and  LQ^Py  Appendix  II  shows  that  the  P^  magnitude 
was  actually  an  upper  bound  based  on  no  signal  measurements  and  on 
only  one  noise  measurement  (at  ILPA) .  The  LQ  magnitudes  were  based 
on  eight  noise  measurements  and  one  signal  measurement.  That  one 
measurement  was  made  at  KAAO,  at  a  distance  of  A  =  2.8°.  The  uncer¬ 
tainties  which  are  introduced  into  the  network-averaged  magnitude 
calculations  by  close-in  surface-wave  data  such  as  this  have  already 
been  discussed,  and  it  would  probably  have  been  better  to  delete  this 
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single  observation  from  the  data  base.  This  deletion  would  have  left 
no  discrimination  variables  applicable  to  this  event,  and  it  would 
have  been  classified  as  an  unknown.  Even  if  the  LQ  magnitudes  had 
been  retained,  the  two  spectral  ratios  should  be  regarded  as  unre¬ 
liable  since  the  magnitude  was  computed  using  a  single  measurement, 
and  one  which  furthermore  was  noise  rather  than  signal. 

3)  Event  79.  This  explosion  was  misclassif ied  using  a  seven-variable 
discriminant  including  the  two  complexity  measurements  comp^  and 
comp2.  It  was  noted  on  Figure  18b  that  comp2  is  anomalously  large 
for  this  event,  and  Appendix  II  shows  that  the  value  used  for  comp2 
was  based  upon  a  single  observation.  When  the  complexity  variables 
are  deleted  from  the  discriminant,  this  event  is  still  misclassif ied 
as  an  earthquake,  but  less  strongly  so.  Indeed,  comparing  the  mag¬ 
nitudes  listed  in  Appendix  I  for  this  event  with  the  short— period 
discriminants  plotted  in  Figures  16a-c  and  Figures  17a-b  shows  that 
event  79  is  at  least  weakly  earthquake-like  with  regard  to  each  dis¬ 
criminant,  and  Table  XVI  confirms  this.  This  event  was  misclassif ied 
because  no  long-period  discrimination  parameters  could  be  used,  and, 
as  Table  XV  demonstrates,  explosions  tend  to  be  misclassif ied  in  the 
absence  of  M  im^  data.  In  the  next  section  of  this  report  we  shall 
show  how  the  short-period-only  discriminants  can  be  made  more  effective 
for  identifying  explosions. 

4)  Event  145.  The  NEIS  depth  for  this  event  is  540  km.  When  complexity 
is  deleted  from  the  eight-variable  discriminant  used  for  this  event, 
it  is  classified  as  an  earthquake. 

5)  Event  147.  The  NEIS  depth  for  this  event  is  479  km.  It  was  classi¬ 
fied  on  the  basis  of  the  fourteen-variable  discriminant,  and  it  was 
explosion-like  with  regard  to  each  of  the  variables.  In  practice, 
however,  this  event  would  have  been  classified  as  an  earthquake, 
since  visual  inspection  of  the  short-period  P— wave  seismogram  recorded 
at  CHTO  (Figure  19)  reveals  a  strong  dilatant  first  motion. 

6)  Event  271.  It  was  noted  in  Figure  18  that  this  shot  array  has  a 
large  complexity.  When  the  three  complexity  variables  are  deleted 
from  the  fourteen-variable  discriminant  which  was  applied  to  it,  it 
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Figure  19.  Short-period  P  seismogram  of  event  147  recorded  at  CHTO. 
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is  classified  as  an  explosion. 


' 'i  As  previously  noted,  event  275  was  deleted  from  the  classification 
because  no  seismograms  for  that  event  were  retained  in  the  data  base. 
It  thus  remains  unknown. 


it  is  thus  seen  that  the  two  principal  causes  of  misclassif ication  were 
insufficient  detections  for  both  earthquakes  and  explosions  and  depth  of  the 
hypocenter  for  earthquakes.  Earthquakes  at  greater  than  crustal  depths  are 
seen  to  resemble  explosions  in  many  respects,  and  this  resemblance  is  especi¬ 
ally  strong  with  regard  to  complexity.  The  systematic  misclassification  of 
deep  events,  particularly  by  means  of  the  P-wave  signal  complexity,  has  been 
demonstrated  in  several  previous  investigations  (Dahlman  and  Israelson,  1977). 


Special  attention  should  be  paid  to  the  classification  of  the  three  mul¬ 
tiple  explosions,  events  53,  265,  and  271.  Figures  16  and  17  show  that  the 
spectral  ratios  for  these  events  lie  in  every  case  within  the  explosion  popu¬ 
lations,  but  as  we  have  seen,  event  271  was  misclassif ied .  Figure  15  shows 
that  the  earthquake-like  aspect  of  each  of  these  events  is  the  complexity.  It 
is  thus  seen  that  complexity  may  lead  to  the  misclassification  of  not  only 
deep  earthquakes  but  also  multiple  explosions.  This  suggests  that  in  practice 
perhaps  each  discriminant  function  should  be  applied  to  unknown  events  twice, 
once  with  and  one  without  the  complexity  variables.  A  classification  criterion 
might  then  be  established  as  follows: 


with  complexity 

earthquake 

earthquake 

explosion 

explosion 


without  complexity  classification 
earthquake  earthquake 

explosion  possible  multiple  explosion 

earthquake  possible  deep  earthquake  (33) 

explosion  explosion 


Application  of  such  a  criterion  would  of  course  fail  to  eliminate  all  deep 
earthquakes  from  the  list  of  events  classified  as  explosions,  since  earthquakes 
such  as  event  147  appear  explosion- like  with  respect  to  M  and/or  P-wave 


Dahlman,  R. ,  and  H.  Israelson  (1977).  Monitoring  Underground  Nuclear  Ex¬ 
plosions  ,  Elsevier  North  Holland  Inc,. 
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spectral  ratios  as  well  as  complexity.  Similarly,  for  certain  configurations 
of  multiple  explosions,  the  spectral  ratios  of  long-to-short  period  variables 
as  well  as  complexity  could  be  made  to  resemble  earthquakes.  In  such  cases, 
however,  the  short-period  P-wave  spectral  ratios  would  still  appear  explosion¬ 
like  . 


DISCUSSION 


It  was  explained  in  a  previous  section  of  this  report  that  the  multivariate 
discrimination  algorithm  is  a  stepwise  process  and  that  at  each  step  in  the 
computation  of  an  n-variable  discriminant  a  search  is  made  over  the  n-m  unused 
variables  in  order  to  find  the  one  which  will  most  strongly  improve  the  dis¬ 
crimination  capability  when  added  to  the  m  variables  already  entered  into  the 
discriminant.  After  the  new  m  +  1  variable  discriminant  is  computed,  the  m  +  1 
variables  are  ranked  in  order  of  their  contribution  to  the  separation  of  the 
earthquake  and  explosion  populations,  and  then  a  new  search  is  made  over  the 
remaining  n-m-1  variables.  The  results  of  this  procedure  for  the  fourteen- 
variable  case  are  shown  in  Table  XX.  It  should  be  noted  that  the  results  would 
be  somewhat  different  in,  for  example,  the  ten-variable  case  because  the  train¬ 
ing  set  would  be  different  and  because  the  relative  importance  of  each  of  the 
ten  common  variables  would  be  affected  by  the  absence  of  the  other  four. 

Because  the  fourteen-variable  discriminant  was  judged  to  be  the  "best"  one  ex¬ 
amined,  however,  the  ranking  of  the  fourteen  variables  shown  in  Table  XX  will 
be  assumed  to  impart  the  best  information  available  about  the  relative  value 
of  each  of  the  individual  discrimination  variables.  The  table  shows  that  the 
order  in  which  the  variables  are  entered  into  the  discriminant  function  differs 
from  the  order  in  which  they  are  ranked  following  the  final  step.  This  dif¬ 
ference  is  to  be  expected,  since  after  the  final  step  the  importance  of  each 
variable  is  compared  with  that  of  all  the  variables  which  were  eventually  en¬ 
tered  into  the  discriminant,  rather  than  with  that  of  just  those  which  had 
been  entered  before  it.  The  final  ranking  thus  presents  the  best  measure  of 
the  relative  value  of  each  of  the  individual  variables,  but  the  order  in  which 
the  first  m  variables  are  entered  into  the  discriminant  function  serves  to 
define  the  best  possible  m-variable  discriminant.  (It  should  be  emphasized 
that  in  this  context  the  "best  possible"  discriminant  is  the  one  which  yields 
the  best  separation  between  the  earthquakes  and  explosions  of  the  51  event 
training  set  used  for  the  fourteen-variable  discriminant.  This  "best  possible" 
m-variable  discriminant  may  be  of  limited  applicability  and/or  effectiveness 
as  a  discriminant  for  the  other  81  events  in  the  data  base.)  By  this  criterion 
it  may  be  seen  that  the  best  single-variable  discriminant  is  the  spectral  ratio 
of  middle-frequency  Love  waves  to  high-frequency  P  waves.  When  all  fourteen 
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TABLE  XX 


Ranking  of  the  Variables  in  the  Fourteen-Variable  Discriminant 
Order  in  which  the  variables  are  entered: 


number 

variable  entered 

F-statistic 

1 

LQ2  -2.31  P3 

79.9171 

2 

P2  -1.72  P3 

7 . 6636 

3 

no(-2)  +  3  fc  (-2) 

2.6590 

4 

comp2 

1.6107 

5 

LQ  -2.30  P3 

0.4096 

6 

LR  -2.31  P3 

1.1904 

7 

P  -1.41  P2 

0.6782 

8 

P  -1.77  P3 

0.3763 

9 

comp^ 

0.3536 

10 

LR3  -2.27  P3 

0.1770 

11 

LQ3  -2.30  P3 

0.1155 

12 

comp3 

0.1086 

13 

C2„(-3)  +  3  fc  (-3) 

0.0136 

14 

LR2  -2.31  P3 

0.0009 

After  all  14  variables  are  entered: 


degrees  of  freedom 

1,49 

1,48 

1,47 

1,46 

1,45 

1,44 

1,43 

1,42 

1,41 

1,40 

1,39 

1,38 

1,37 

1,36 


Order 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


variable 
LQ2  -2.31  P3 
LQX  -2.30  P3 
LR1  -2.31  P3 
P  -1.41  P2 
fl0(-2)  +  3  fc  (-2) 
F  -1.77  P3 
comp^ 

LR3  -2.27  P3 
LQ3  -2.30  P3 
comp3 
comp2 

P2  -1.7"  P3 
0o(-3)  +  3  fc  (-3) 
LR„  -2.31  P„ 


F-statistic  (degrees  of  freedom  —  1, 
2.3203 
1.4501 
1.2004 
0.7366 
0.4739 
0.3159 
0.1798 
0.1535 
0.0930 
0.0915 
0.0783 
0.0544 
0.0114 
0.0009 
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variables  are  used,  the  discrimination  algorithm  assigns  the  top  two  rankings 
to  the  spectral  ratios  LQ.,^  and  LC^^  and  the  third  highest  ranking  to  the 
spectral  ratio  LR^P^  It  is  thus  seen  that,  as  is  well  known,  is  a 

valuable  discriminant,  but  it  happens  that  in  this  experiment  its  most  im¬ 
portant  form  is  LQ:P  rather  than  the  traditional  LR:P.  It  is  perhaps  not 
surprising  that  Love  waves  are  a  more  effective  discriminant  than  are  Rayleigh 
waves,  since  in  an  idealized  case  Love  waves  should  be  an  infallible  dis¬ 
criminant,  their  very  existence  being  diagnostic  of  an  earthquake.  The 
Rayleigh-wave  radiation  of  an  idealized  explosion,  by  contrast,  is  not  more 
than  an  order  of  magnitude  smaller  than  that  of  a  shallow  earthquake  with  the 
same  seismic  moment  (von  Seggern  and  Rivers,  1979).  In  practice,  however. 

Love  waves  are  emitted  by  explosions  due  to  mode  conversion  and  to  deviation 
of  the  source  mechanism  from  a  purely  radial  presure  pulse,  but  they  are  still 
emitted  in  small  enough  amplitudes  that  they  remain  a  better  discriminant  than 
Rayleigh  waves.  In  fact,  even  though  the  spectral  ratio  of  middle-frequency 
Love  waves  to  high-frequency  P-waves  was  the  best  single-variable  discriminant, 
Table  XX  shows  that  the  corresponding  ratio  for  Rayleigh  waves  is  the  last 
variable  entered  into  the  fourteen-variable  discriminant  function.  It  should 
not  be  concluded  from  this  ranking  that  measured  by  means  of  Rayleigh 

waves  is  a  poor  discriminant;  Figure  16j  shows  that  this  is  clearly  not  the 
case.  Rather,  what  this  ranking  indicates  is  that  after  (LQ^  -2.31  P3>  has 
been  entered  into  the  discriminant  and  assigned  heavy  weight,  most  of  the 
information  contained  in  the  (LR2  -2.31  P3)  variable  becomes  redundant,  and 
the  discrimination  algorithm  regards  this  variable  as  superfluous.  It  should 
be  remembered  that  the  surface-wave  magnitudes  used  in  the  discrimination  are 
averages  taken  over  a  network  of  many  stations;  for  data  from  a  single  station, 
Love  waves  and  Rayleigh  waves  are  likely  to  be  complementary  rather  than  re¬ 
dundant,  since  for  many  source  geometries  the  nodes  of  one  type  of  surface- 
wave  radiation  pattern  are  coincident  with  the  lobes  of  the  other  type  of 
pattern.  Attention  should  be  paid  to  the  ranking  in  Table  XX  of  the  three 
complexity  variables.  It  has  already  been  mentioned  that  discrimination  based 
on  complexity  can  be  misleading,  particularly  for  deep  earthquakes  and  for 
multiple  explosions. 
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No  assessment  was  made  of  the  importance  of  the  six  variables  which  were 
not  employed  in  the  fourteen-variable  discriminant  function,  namely  the  spectral 
ratios  to  of  the  six  magnitudes  of  the  long-period  body  waves.  Although 
Tables  XIV  and  XV  show  that  these  variables  were  observed  too  seldom  to  permit 
the  accurate  determination  of  any  multivariate  discriminant  function  contain¬ 
ing  them,  Figures  16c-h  show  that  each  of  these  six  variables  considered  in¬ 
dividually  shows  considerable  discrimination  capability.  Some  of  this  capa¬ 
bility,  of  course,  is  due  to  the  fact  that  long-period  body  waves  are  usually 
detected  only  for  large  events,  and  in  such  cases  many  different  discriminants 
are  usually  effective.  Nevertheless,  these  six  variables  apparently  would 
supply  valuable  information  to  any  multivariate  discriminant  function  containing 
them,  and  in  any  future  investigations  of  statistical  seismic  discrimination 
techniques  they  should  be  routinely  measured  for  all  events  in  order  to  build 
up  a  suitable  training  set  of  events  to  which  they  are  applicable. 

It  was  noted  in  the  previous  section  of  this  report  that  of  the  earth¬ 
quakes  misclassif ied  as  explosions,  several  are  deep  events.  These  misclassi- 
fications  were  undoubtedly  influenced  by  noise  and  missing  data,  as  in  the 
cases  of  the  other  misclassif ications,  but  the  possibility  remains  that  system¬ 
atic  errors  were  made  which  affected  deep  events  more  severely  than  shallow 
ones.  One  such  systematic  error  could  be  the  use  of  zero-depth  B-factors  in 
the  computation  of  the  body-wave  magnitudes  P^,  P^,  P^»  LP^,  T^,  LP^,  LS^, 

LS^  and  LS^.  Another  would  be  the  use  of  a  zero-depth  geometrical  spreading 
correction  in  the  computation  of  the  low-frequency  spectral  levels  fi0(-2) 
and  fi0(-3).  Although  these  errors  are  systematic,  their  effect  is  unpredictable, 
being  dependent  upon  the  epicenter-to-station  distances  for  the  particular 
network  of  stations  used  for  a  given  event.  The  hypothesis  of  systematic 
error  is  thus  not  necessarily  incompatible  with  the  fact  that  several  deep 
events  were  correctly  identified  as  earthquakes.  As  has  been  pointed  out, 
there  could  also  be  certain  real  differences  between  deep  earthquakes  and 
shallower  ones  such  as  comprised  most  of  the  training  set.  Certainly  the 
surface-wave  excitation  is  less  for  deep  events.  The  P-wave  spectra  (and 

hence  the  spectral  ratios  P1:P2>  Pi:P3’  and  P2:IV  raay  also  be  different  for 
deep  earthquakes  on  account  of  certain  features  of  their  source  mechanisms. 

For  example,  at  least  some  deep  earthquakes  may  involve  smaller  source 
dimensions  and  larger  stress  drops  than  do  those  at  shallower  depths,  causing 
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them  to  have  higher  corner  frequencies  and  more  high-frequency  spectral  con¬ 
tent.  Their  mechanisms  might  also  be  explosion-like  in  that  they  could  involve 
a  volume  change  due  to  a  mineralogic  phase  transition.  The  spectra  of  deep 
earthquakes  might  also  differ  from  those  of  the  other  earthquakes  m  the  train¬ 
ing  set  on  account  of  differences  in  the  propagation  paths  of  the  signals;  P 
waves  from  deep  events  travel  through  the  earth's  crust  and  upper  mantle  only 
once,  not  twice.  The  P-wave  signal  complexity  is  also  affected  by  the  differ¬ 
ences  in  propagation  paths,  and  it  has  already  been  shown  that  the  anomalously 
low  values  of  the  complexity  for  deep  events  can  be  sufficent  to  cause  mis- 
classif ication. 

In  an  actual  operational  procedure  of  detecting,  locating,  and  discrimi¬ 
nating  unknown  events,  however,  events  which  had  been  determined  to  be  deep 
would  be  deleted  from  the  list  of  unknowns  prior  to  the  application  of  the 
multivariate  discriminant  functions.  In  order  to  evaluate  the  effect  of  using 
hypocentral  depth  as  a  discriminant,  a  three-stage  discrimination  process  was 
performed  which  resembled  ones  described  by  Dahlman  and  Israelson  (1977).  The 
results  of  this  three-stage  process  are  shorn  in  Table  XXI.  Listed  m  this 
table  are,  in  each  instance  in  which  they  are  available,  the  hypocentral  depths 
in  km  for  every  event  as  they  are  given  in  the  bulletins  produced  by  the  NEIS 
and  by  the  Network  Event  Processor  (NEP)  at  the  Seismic  Data  Analysis  Center. 

No  values  are  given  for  those  events  whose  depths  were  restricted  by  the  NEP 
to  be  zero.  The  first  stage  of  the  discrimination  consists  of  classifying  as 
earthquakes  all  events  whose  depths  are  believed  to  be  greather  than  100  loin 
Table  XXI  shows  that  there  are  24  such  events.  (One  of  them,  event  50,  was 
restricted  by  the  NEP  to  a  depth  of  40  kra,  but  was  classified  as  an  earthquake 
on  the  basis  of  its  NEIS  depth  of  400  km.)  The  second  stage  of  the  discrimi¬ 
nation  was  the  application  to  the  other  109  events  of  those  multivariate  dis¬ 
criminant  functions  which  contained  at  least  one  of  the  six  M^.n^  variables, 
numbers  15  through  20  in  Table  XIV.  This  stage  was  somewhat  different  from 
that  of  Dahlman  and  Israelson  (1977)  in  that  these  multivariate  discriminants 
contain  not  only  M^m^  but  also  the  short-period  spectral  ratios,  complexity, 
and  low-frequency  level  versus  corner  frequency.  Table  XXI  shows  that  the 
application  of  these  discriminants,  taken  from  Tables  XVII  and  XVIII,  results 
in  the  classification  of  82  more  events,  two  of  them  (events  8  and  271)  in¬ 
correctly.  The  third  stage  consists  of  the  classification  of  the  remaining 
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TABLE  XXI  (Continued) 
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26  events  using  only  short-period  discriminants.  One  of  these  events  was 
mis classified.  It  should  be  noted  that  the  training  sets  used  in  calculating 
the  discriminants  applied  in  the  second  and  third  stages  included  the  24  events 
identified  in  stage  one  as  being  deep.  It  is  doubtful  that  using  only  shallow 
events  in  the  training  sets  would  improve  the  discrimination,  because  some  o£ 
the  unknown  events  remaining  after  stage  one  are  actual  y  deep,  and  shallow- 
event  discriminants  might  systematically  misclassify  them.  In  comparing  Table 
XXI  of  this  report  with  Table  10.4  of  Dahlman  and  Israelson  (1977),  it  should 
be  remembered  that  short-period  discriminants  were  applied  in  both  the  second 
and  third  stages  in  this  report  but  in  only  the  final  stage  by  Dahlman  and 
Israelson.  It  ought  to  be  noted  that  several  of  the  24  events  which  could  be 
identified  as  earthquakes  on  the  basis  of  their  depth  were  assigned  uncertain 
classifications,  denoted  by  asterisks,  in  Tables  XVII  and  XVIII. 

So  far  in  the  calculation  of  all  the  discriminant  functions  examined  in 
this  report  it  has  been  assumed  that  the  misclassif ication  of  earthquakes  is 
as  undesirable  as  that  of  explosions.  In  practice,  unidentifed  explosions  are 
a  more  serious  error  than  are  false  alarms,  so  one  might  wish  to  adjust  the 
discriminant  functions  to  decrease  the  number  of  misclassif ied  explosions,  even 
if  such  adjustment  increases,  by  a  greater  amount,  the  number  of  misclassif ied 
earthquakes.  Such  an  adjustment  can  be  made  by  assigning  more  weight  in  the 
discrimination  algorithm  to  explosions.  This  assignation  of  weights  may  be 
interpreted  geometrically  as  moving  the  discriminant  hyperplane  parallel  to 
itself  in  the  direction  of  the  earthquake  population  so  that  all  explosions, 
including  the  outliers,  lie  on  the  same  side  of  the  hyperplane’s  new  position. 
For  example,  the ' three-variable  discriminant  function  which  uses  only  the  short 
period  P-wave  spectral  ratios  is  given  by 

D(Vr  V2,  V3>  -  -1.863  Vx  -1.727  V.,  -1.366  V.,  +  const.  (34) 

where  (cf.  Table  XIV) 

V  =  P  p2 

11  (35) 

V  =  P,  -1.72  P-  V  ' 

V3  -  F,  -1.77  P3 

The  value  of  the  constant  in  equation  (34)  was  -6.602  in  the  three-variable 
discriminant  function  which  was  used  for  classifying  event  155  in  Table  XVIII. 
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As  is  shov. there,  this  discriminant  is  highly  unreliable,  misclassify ing 
six  explosions  and  four  earthquakes.  The  effect  upon  this  discriminant  of 
assigning  more  weight  to  explosions  is  shown  in  Table  XXII.  It  is  demon¬ 
strated  that  by  changing  the  constant  from  -6.602  to  -5.321,  the  total  number 
of  mis classifications  is  raised  from  ten  to  sixteen,  but  the  number  of  mis- 
classi.ied  explosions  is  reduced  to  two.  By  allowing  the  constant  to  take  on 
some  intermediate  value  (i.  e.,  by  assigning  an  intermediate  weight  to  explo¬ 
sions)  it  is  possible  to  satisfy  some  arbitrary  criterion  chosen  as  the  most 
acceptable  trade-off  between  the  number  of  undetected  explosions  and  the  number 
of  false  alarms.  In  an  actual  discrimination  routine  it  would  probably  be 
advantageous  to  choose  some  such  optimum  trade-off  and  then  weight  the  two 
populations  in  the  discriminant  functions  accordingly. 

The  multivariate  discriminants  which  have  thus  far  been  described  used 
a  training  set  of  between  51  and  121  events.  The  question  will  now  be  addressed: 
since  all  the  unknown  events  (except  the  multiple  explosions)  had  in  fact 
been  assigned  a  priori  classifications,  would  the  discrimination  have  worked 
as  well  on  a  suite  of  true  unknowns?  In  order  to  investigate  this  question, 
certain  of  the  multivariate  discriminant  functions  weire  re-computed  using  only 
a  fraction  of  events  previously  used  in  the  training  set.  These  new  discrimi¬ 
nants  ware  then  applied  to  all  the  events,  including  those  which  were  not  used 
in  the  training  set  and  which  may  therefore  be  regarded  as  true  unknowns.  If 
the  discrimination  works  almost  as  well  using  a  small  training  set,  then  the 
discriminant  functions  are  fairly  stable,  and  the  results  described  in  this 
report  may  be  applied  to  an  actual  discrimination  routine  with  more  confidence 
than  they  could  be  otherwise.  It  is  shown  in  Table  XXIII  that  the  four  variable 
discriminant  was  in  fact  rather  stable  against  perturbations  introduced  by 
changing  the  training  set.  The  performance  of  this  discriminant  was  not  greatly 
deteriorated  by  the  use  of  69  or  even  34  events  in  the  training  set  rather  than 
the  full  suite  of  103  events.  Table  XXIV  shows  that  the  classifications  based 
on  the  eight-variable  discriminant  were  scarcely  changed  when  the  training  set 
was  reduced  from  71  events  to  23,  but  this  is  misleading,  since  the  misclassi- 
fication  rate  tripled  when  48  events  were  used.  It  is  not  surprising  that 
an  eight-variab1 c  discriminant  function  is  poorly  determined  by  48  events,  and 
it  was  on  account  of  the  likelihood  of  a  similarly  poor  determination  that  no 
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TABLE  XXII 


Effect  on  a  Three-variable  Discriminant 
Function  of  Assigning  More  Weight  to  Explosions 

discrimination  variables  (cf.  Table  XIV)  :  1,  2,  3 
constant  term  in  the 

discriminant  function  classification 


shot 

quake 

-6.602 

shot 

14 

6 

quake 

4 

84 

shot 

quake 

-6.219 

shot 

15 

5 

quake 

4 

84 

shot 

quake 

-5.968 

shot 

16 

4 

quake 

6 

82 

shot 

quake 

-5.740 

shot 

17 

3 

quake 

8 

80 

shot 

quake 

-5.321 

shot 

18 

2 

quake 

14 

74 
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TABLE  XXIII 


Effect  on  the  Four-variable  Discriminant 


Function  of  Changing 

the  Size  of  the  Training  Set 

number  of  events  in 

the  training  set 

classification 

shot 

quake 

7  shots 

shot  14 

7 

27  quakes 

quake  6 

76 

shot 

quake 

14  shots 

shot  14 

7 

55  quakes 

quake  4 

78 

shot 

quake 

21  shots 

shot  16 

5 

82  quakes 

quake  4 

78 

complete  training  set 
delete  event  189  only 
delete  event  195  only 


probability  of  being  an  explosion: 


event  189  event  195 


0.500  0.012 
0.468  0.011 
0.498  0.013 
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TABLE  XXIV 


Effect  on  the  Eight-variable  Discriminant 
Function  of  Changing  the  Size  of  the  Training  Set 


number  of  events  in 
the  training  set 

5  shots 
18  quakes 

11  shots 
37  quakes 

16  shots 
55  quakes 


complete  training  set 
delete  event  189  only 
delete  event  195  only 


classification 


shot  quake 

shot  15  1 

quake  3  52 

shot  quake 

shot  11  5 

quake  4  51 


shot  quake 
shot  15  1 

quake  2  53 

probability  of  being  an  explosion: 

event  189  event  195 

0.696  0.001 

0.595  0.001 

0.692  0.001 
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use  was  made  of  the  discriminants  involving  fifteen  or  more  variables  but 
determined  by  38  or  fewer  events  (cf.  Table  XV).  Tables  XXIII  and  XXIV  thus 
show  that  at  least  the  four-and  eight-variable  discriminants  are  fairly  well 
determined  and  may  be  used  to  classify  events  not  contained  in  the  training  set 
with  about  the  same  degree  of  success  as  for  the  training  set  events. 

The  tables  also  show  the  results  of  another  test  of  discriminant  stability, 
namely  the  deletion  of  specific  single  events  from  the  training  set.  The  events 
which  were  selected  for  deletion  were  number  189,  marginally  classified  by 
both  the  four-and  eight-variable  discriminants  as  an  explosion,  and  number  195, 
strongly  classified  by  both  as  an  earthquake.  Table  XXIV  shows  that  had  either 
event  not  been  assigned  an  _a  priori  classification,  it  would  still  have  been 
correctly  classified  by  the  eight-variable  discriminant  computed  on  the  basis 
of  the  other  80  events  in  the  training  set.  Table  XXIII  shows,  however,  that 
the  deletion  of  event  189  from  the  training  set  of  the  four-variable  discrimi¬ 
nant  function  results  in  the  misclassif ication  of  that  event  as  an  earthquake. 
This  test  shows  that,  as  one  would  expect,  the  classification  of  data  points 
lying  close  to  the  discriminant  hyperplanes  can  change  when  the  positions  of 
the  hyperplanes  shift  slightly,  and  such  shifts  do  in  fact  occur  when  data 
points  close  to  the  hyperplanes  are  added  to  or  are  deleted  from  the  training 
set.  It  was  for  this  reason  that  the  classification  of  events  in  Table  XVII 
and  XVIII  was  noted  as  uncertain  if  the  probability  of  being  an  explosion  was 
between  0.2  and  0.8.  In  an  actual  discrimination  process,  then,  in  which  event 
189  was  a  true  unknown  and  hence  was  assigned  no  a_  priori  classification,  the 
four-variable  discriminant  would  classify  that  event  marginally  as  an  earth¬ 
quake.  Had  more  events  been  used  in  the  training  set,  however,  the  discrimi¬ 
nant  function  might  have  been  slightly  different  (such  slight  differences  are 
not  incompatible  with  the  previous  assertion  of  discriminant  stability),  and 
event  189  might  have  been  classified  margl.nally  as  an  explosion.  An  operative 
discrimination  routine  would  be  a  learning  process  in  which  newly-classified 
events  would  be  added  periodically  to  the  training  set  and  the  discriminant 
functions  would  then  be  re-calculated.  The  a_  posteriori  classification  result¬ 
ing  from  this  re-calculation  might  reveal  that  some  of  the  previous  classifi¬ 
cations  of  training  set  events  are  wrong  or  at  least  uncertain.  Events  with 
marginal  classifications  like  event  189  should  not  be  added  to  the  training 
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set  unless  their  classifications  could  be  verified  by  independent  information 
such  as  dilatational  first  motions,  photographic  evidence  of  cratering,  etc. 
After  the  discriminant  functions  are  re-calculated  they  should  then  be  re¬ 
applied  to  the  events  with  marginal  classifications,  and  the  resulting  new 
classifications  should  be  regarded  as  more  reliable. 

Although  the  tests  which  are  shown  in  Tables  XXIII  and  XXIV  demonstrate 
that  the  discriminant  functions  are  fairly  insensitive  to  the  total  number  of 
events  used  in  the  training  set,  there  remains  the  possibility  that  the  dis¬ 
criminants  are  sensitive  to  the  geographical  distribution  of  these  events. 

No  assessment  can  be  made  of  the  value  of  the  discrimination  techniques  until 
it  is  determined  to  what  extent,  if  any,  the  discrimination  between  the  earth¬ 
quake  and  explosion  populations  was  actually  discrimination  between  certain 
characteristics  of  the  source  regions  of  the  two  populations.  If  the  dis¬ 
crimination  is  due  to  epicenter-to-receiver  path  effects  rather  then  to  source 
mechanisms,  it  will  lose  its  value  when  applied  to  explosions  in  seismically 
active  regions  or  to  earthquakes  at  test  sites.  Some  assurance  that  the  dis¬ 
crimination  was  not  solely  a  path  effect  may  be  gained  by  noticing  that  the 
data  base  did  in  fact  include  at  least  two  examples  of  earthquakes  and  ex¬ 
plosions  occurring  close  to  each  other.  Table  I  shows  that  event  21,  an  ex¬ 
plosion,  occurred  at  a  distance  of  about  150  km  from  event  65,  an  earthquake. 

The  multivariate  discriminants  correctly  classified  both  of  these  Lake  Baikal 
events.  Table  I  also  shows  that  event  264,  an  earthquake,  occurred  close  to 
the  Semipalatinsk  explosions.  This  event  was  also  classified  correctly,  so 
it  seems  that  path  effects  did  not  dominate  the  multivariate  discrimination. 
Certain  individual  discriminants  may  have  been  affected  by  differences  between 
the  earthquake  and  explosion  source  regions,  however;  in  particular,  it  has 
already  been  noted  that  attenuation  differences  may  have  enhanced  the  discrimi¬ 
nation  power  of  the  short-period  P-wave  spectral  ratios.  The  differences  in  the 
attenuation  coefficient  t*  for  various  source  region-to-station  paths  are 
shown  in  Appendix  IV. 

It  has  been  emphasized  that  the  training  set  of  events  used  for  this  re¬ 
port  was  not  geographically  homogeneous.  The  question  is  now  addressed  of 
whether  there  does  indeed  exist  a  single  earthquake  population  and  a  single 
explosion  population  for  all  Eurasian  events,  as  has  been  assumed  in  this 
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report,  or  whether  there  are  in  fact  several  distinct  populations  charac¬ 
teristic  of  different  regions.  If  region-dependent  populations  exist,  then 
the  same  discriminant  functions  may  not  work  equally  well  for  all  of  them. 

There  may  in  fact  exist  anomalous  regions  for  which  the  earthquake  population 
has  a  large  intersection  with  the  total  "Eurasian"  explosion  population  or  for 
which  the  explosions  are  contained  within  the  "Eurasian"  earthquake  population. 

It  must  be  known  whether  the  discrimination  is  region-dependent  before  the 
results  of  this  report  may  be  applied  to  the  classification  of  events  occurring 
in  regions  not  adequately  sampled  by  the  129  events  which  were  used  in  the 
training  set.  In  order  to  test  for  the  existence  of  regional  sub-populations 
of  the  "Eurasian"  earthquake  population,  a  test  was  made  of  whether  earthquakes 
in  different  tectonic  regions  may  be  discriminated  from  each  other.  This 
procedure  will  first  be  described,  and  then  the  results  will  be  discussed. 

Throughout  this  report  it  has  been  assumed  that  the  process  of  discrimi¬ 
nation  consists  of  classifying  unknown  events  as  being  members  of  one  or  the 
other  of  exactly  two  groups.  Inspection  of  the  formulas  for  the  calculation 
of  the  discriminant  functions  (von  Seggern  and  Rivers,  1979),  however,  shows 
that  it  is  a  straightforward  process  to  generalize  the  discrimination  to  in¬ 
clude  an  arbitrarily  large  number  of  groups,  and  the  computer  program  (Jennrich, 
1977)  which  was  used  for  the  multivariate  analysis  has  the  capability  of  per¬ 
forming  multi-population  discrimination.  Discrimination  between  several 
groups  involves  the  calculation  of  a  different  multivariate  classification 
function  for  each  pair  of  groups,  a  process  which  is  geometrically  equivalent 
to  the  determination  of  a  set  of  hyperplanes  whose  intersection  partitions 
space  into  cells,  one  cell  per  group*  The  classification  function  may  then  be 
applied  to  unknown  events  in  order  to  determine  in  which  cell  they  lie,  i.e., 
to  which  population  they  belong.  This  process  ought  not  to  be  confused  with 
cluster  analysis.  Cluster  analysis  attempts  to  find  how  many  statistically 
significantly  different  populations  exist  within  a  given  training  set,  and  then 
unknown  events  are  assigned  to  one  or  another  of  those  thus  determined  populations 
("clusters").  Discriminant  analysis,  on  the  other  hand,  begins  with  an  a_ 
jjr ior i  specification  of  the  number  of  populations  and  of  which  events  in  the 
training  set  belong  to  each  one,  and  then  unknown  events  are  assigned  to  one 
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or  another  of  those  given  populations.  Cluster  analysis  thus  determines  the 
most  statistically  significant  partitioning  of  the  training  set  into  groups, 
and  the  investigator  must  attempt  to  infer  the  physical  basis  for  the  simi¬ 
larity  between  the  events  in  each  group  (explosions,  shallow  dip-slip  earth¬ 
quakes,  earthquakes  in  Tibet,  etc).  In  discriminant  analysis  the  investigator 
specifies  the  criteria  whereupon  the  events  are  to  be  grouped  (explosions, 
strike-slip  earthquakes  in  the  Ural  Mountains,  earthquakes  occuring  on  a 
Tuesday,  etc.)  and  then  the  analysis  reveals  whether  the  groups  are  statistically 
significant,  i.  e.,  non-overlapping. 

In  order  to  evaluate  the  influence  of  the  geographical  distribution  of 
events  upon  their  classification,  a  multi-population  discrimination  was  per¬ 
formed  which  attempted  to  determine  whether  earthquakes  occurring  in  different 
tectonic  regions  were  significantly  different  from  one  another.  The  results 
of  this  discrimination  are  shown  in  Table  XXV.  In  the  first  test  shown,  the 
classification  functions  based  on  the  three  short-period  P-wave  spectral  ratios 
were  applied  to  three  populations  of  events,  namely  earthquakes  in  Kamchatka, 
in  Tibet,  and  in  the  Pamirs-Hindu  Kush  region.  The  table  shows  how  many  events 
were  classified  a  posteriori  into  each  group  and  also  which  events  they  were. 

It  can  be  seen  that  with  respect  to  the  three  chosen  variables  the  Kamchatka 
and  Tibet  populations  are  rather  different  from  each  other,  and  the  Pamirs  popu¬ 
lation  is  quite  different  from  both  of  them.  The  second  test  shows  the  addi¬ 
tion  of  the  three  LR:P  spectral  ratios  further  separates  the  Kamchatka  earth¬ 
quakes  from  those  in  Tibet.  Although  these  tests  were,  of  necessity,  based  on 
a  sparse  data  set,  it  appears  that  there  may  indeed  exist  regional  sub-popu¬ 
lations  of  the  large  population  of  earthquakes  which  was  used  for  calculating 
the  multivariate  discriminant  functions.  Unfortunately,  the  data  is  too  scanty 
to  determine  whether  there  also  exist  regional  sub-populations  of  the  explosion 
population.  If  such  sub-populations  do  exist,  then  the  discriminant  functions 
may  differ  significantly  from  region  to  region,  in  which  case  the  results  of 
this  report  cannot  be  applied  confidently  to  the  discrimination  of  other  suites 
of  events.  It  is  only  when  the  discriminant  functions  are  understood  quanti¬ 
tatively  on  the  basis  of  the  physics  of  event  source  mechanisms  and  propagation 
path  effects,  rather  than  just  on  an  empirical  basis  as  in  this  report,  that 
they  may  be  applied  to  the  classification  of  events  occurring  in  regions  not 

sampled  by  the  training  set. 
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TABLE 


It  has  already  been  pointed  out  that  systematic  errors  may  have  occurred 
in  the  computation  of  the  magnitudes  of  deep  events.  Certain  other  systematic 
errors  may  have  been  made  which  affected  all  the  events  equally.  Specifically, 
no  station  corrections  were  applied  for  any  station,  and  magnitude  biases  may 
have  been  thereby  introduced  into  the  network  averages.  The  station  correc¬ 
tions  which  should  be  applied  for  a  given  station  are  a  function  both  of  the 
phase  (e.g.,  LP  or  LS)  and  of  the  frequency  (e.g.,  or  P3>.  These  correc¬ 

tions  should  be  determined  by  comparing  the  magnitudes  of  each  measured  para¬ 
meter  at  a  single  station  with  the  values  obtained  by  averaging  over  the  net¬ 
work.  Unfortunately,  so  many  measurements  used  in  computing  the  network  mag¬ 
nitudes  were  noise  levels  and  so  many  different  sub-networks  of  stations  were 
used  for  different  events  that  such  a  calculation  employing  this  data  base 
could  yield  only  questionable  results.  Station  corrections  would  best  be 
determined  by  using  a  data  base  of  large  events ,  which  would  be  detected  with 
a  high  signal- to-noise  ratio  at  all  stations  of  a  network  which  remains  con¬ 
stant  from  event  to  event.  In  any  case,  we  note  that  the  addition  of  station 
corrections  to  the  magnitudes  would  have  changed  the  values  calculated  by  the 
Ringdal  (1976)  algorithm,  perhaps  (at  least  in  some  cases)  by  a  significant 
amount. 

A  subtle  bias  may  have  been  introduced  into  the  calculation  of  the  magni¬ 
tudes  for  many  of  the  events  by  the  geographical  distribution  of  stations  in 
the  network.  Specifically,  the  use  of  all  five  stations  BFAK,  UCAK,  NJAK, 

CNAK,  and  TNAK  may  have  placed  undue  emphsis  on  signals  recorded  in  one  small 
sector  of  the  radiation  pattern  for  events  teleseismic  to  Alaska.  This 
problem  could  be  particularly  acute  if  these  stations  were  close  to  a  node  in 
the  radiation  pattern  for  one  or  more  of  P,  SH,  SV,  LR,  or  LQ.  This  potential 
bias  might  be  investigated  by  computing  network  magnitude  estimates  using 
only  one  of  these  five  stations  and  then  repeating  the  discrimination  analysis 
in  order  to  see  whether  the  discrimination  capability  of  the  network  was  there¬ 
by  significantly  diminished  or,  at  least  for  certain' events ,  somewhat  enhanced. 
Another  way  of  testing  whether  too  much  weight  was  given  to  the  Alaskan  sites 
would  be  to  assign  each  of  the  stations  a  weight  of  one-fifth  when  computing 
the  network  magnitude  estimates.  The  procedure  of  deleting  certain  stations 
from  the  network  and  repeating  the  discrimination  could  also  be  applied  to 
each  of  the  twenty-seven  stations  one  at  a  time  in  order  to  gain  some  assessment 
of  the  contribution  of  each  station  (perhaps  in  a  negative  sense)  to  the 
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network  discrimination  capability  .  This  assessment  could  then  be  used  in 
assigning  weights  to  the  stations,  and  these  in  turn  could  be  used  to  calculate 
new  values  of  the  network  magnitude  estimates.  As  was  pointed  out  in  the  discus¬ 
sion  of  station  corrections,  however,  the  over-abundance  of  noise  measurements 
relative  to  signals  and  the  inconstancy  from  event  to  event  of  the  sub-network 
of  stations  which  were  actually  used  would  tend  to  introduce  a  large  element  of 
uncertainty  into  such  an  assignment  of  weights. 


CONCLUSIONS 


The  following  conclusions  may  be  drawn  from  the  statistical  discrimination 
experiment  described  in  this  report: 

1)  Inclusion  of  SRO/ASRO  stations  in  the  data  base  increases  the  volume 
of  useful  long-period  data,  thereby  enhancing  the  discrimination  cap¬ 
ability.  Stations  KAAO  and  MAIO,  being  close  to  the  Area  of  Interest, 
are  particularly  useful  for  detecting  long-period  signals.  The  beam¬ 
forming  capabilities  of  the  four  large  arrays  make  them  valuable  for 
this  purpose  also. 

2)  Inclusion  of  SRO/ASRO  stations  helps  better  determine  the  parameters 
which  are  used  for  short-period  discrimination.  The  SNR  at  these  sta¬ 
tions  is  good,  but  as  is  shown  in  Figure  1,  weak  signals  are  system¬ 
atically  unreported. 

3)  Whenever  an  event  is  large  enough  for  surface  waves  to  be  detected, 

M  :m,  is  a  valuable  discriminant.  Its  value  is  even  greater  when 

Love  waves  rather  than  Rayleigh  waves  are  used  to  measure  M  . 

s 

4)  If  a  seismogram  contains  no  visible  signal  in  the  data  window  sur¬ 
rounding  a  predicted  arrival  time,  the  level  of  the  background  noise 
in  that  window  should  be  measured  for  use  in  a  maximum-likelihood 
estimate  of  the  magnitude  of  the  unseen  signals. 

5)  Even  if  a  certain  signal  is  undetected  at  every  station  in  the  net¬ 
work,  an  upper  bound  to  the  magnitude  of  that  signal  can  be  computed 
from  noise  measurements.  This  upper  bound  may  (or  may  not)  be  useful 
for  negative  discrimination. 

6)  Deep  earthquakes  are  more  prone  to  misclassif ication  than  are  shallow 
ones.  Hypocentral  locations  and  the  detection  of  depth  phases  should 
therefore  be  used  in  order  to  eliminate  as  many  deep  events  as  possible 
from  the  data  base  of  unknowns  before  the  statistical  discrimination 

is  performed. 

7)  Deep  earthquakes  and  multiple  explosions  are  especially  prone  to  mis¬ 
classif  ication  on  the  basis  of  P-wave  signal  complexity.  Application 
of  the  criterion  shown  in  equation  (33)  may  reduce  the  number  of  mis¬ 
classif  ications  due  to  complexity. 
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8)  As  is  shown  in  Figure  19,  the  observation  of  dilatant  first  motions 
may  be  used  to  identify  as  earthquakes  some  events  which  the  purely 
statistical  discrimination  process  would  misclassify  as  explosions. 

9)  In  order  to  avoid  problems  introduced  by  a  magnitude  bias  between  the 
explosion  and  earthquake  populations,  spectral  ratios  of  the  form 
(magnitude^-coefficient'magnitude,^)  should  be  used  as  discrimination 
variables,  rather  than  the  magnitudes  themselves. 

10)  Multivariate  discrimination  is  preferable  to  any  particular  univariate 
discriminant.  The  multivariate  discrimination  should  employ  as  many 
variables  as  possible  for  any  given  unknown  event,  as  long  as  the  dis¬ 
criminant  function  was  calculated  using  a  suitably  large  training  set 
of  events  whose  classification  is  known  a  priori. 

11)  Statistical  discrimination  should  be  a  learning  process.  When  a  group 
of  unknown  events  has  been  classified  by  the  multivariate  analysis, 

it  should  then  be  added  to  the  training  set  used  to  calculate  the 
discrimination  variables  and  the  multivariate  discriminant  functions. 
As  the  data  base  enlarges,  the  discrimination  becomes  more  reliable, 
and  it  may  then  be  concluded  from  ji  posteriori  classifications  that 
some  of  the  classifications  of  events  in  the  training  set  are  wrong. 
Events  which  could  be  classified  only  marginally  should  not  be  added 
to  the  training  set  unless  their  classifications  can  be  verified  by 
independent  evidence.  As  the  training  set  expands,  classifications 
previously  noted  as  uncertain  should  be  repeated. 
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RECOMMENDATIONS 


The  authors  make  the  following  recommendations  for  future  research  per¬ 
taining  to  seismic  discrimination: 

1)  A  determination  of  station  corrections  should  be  made  for  P-,  S-, 
and  surface  waves  in  every  frequency  band  in  order  to  eliminate 
station  bias  from  the  network  estimates  of  event  magnitudes. 

2)  Weights  should  be  assigned  to  each  station  in  the  network  for  use 
in  the  calculation  of  event  magnitudes.  These  weights  should  be 
determined  on  the  basis  of  the  scatter  in  the  magnitude  calculation, 
i.  e.,  the  repeatibility  of  measurement,  at  each  station.  Consider¬ 
ing  a  large  group  of  events  having  the  same  network-averaged  magni¬ 
tude,  the  station  bias  at  some  given  station  is  the  difference  be¬ 
tween  that  common  magnitude  and  the  average  of  the  individual  mag¬ 
nitudes  measured  at  that  station  for  the  events  in  the  group,  and  the 
weight  which  should  be  assigned  to  the  station  is  inversely  propor¬ 
tional  to  the  size  of  the  "error  bars"  around  that  average  magnitude. 

3)  Because  the  individual  discriminants  employed  in  this  study  have  been 
empirically  determined,  an  investigation  should  be  made  of  their 
physical  basis.  Such  an  investigation  would  attempt  to  determine 
circumstances  under  which  particular  discriminants  might  be  useless 
or  misleading,  as  was  the  case,  for  example,  with  complexity  for  deep 
earthquakes.  This  study  would  address  the  question  of  whether  the 
discriminant  lines  shown  in  Figures  15  and  16  are  in  fact  straight 
lines  or  whether  they  bend  at  low  (or  high)  magnitudes,  whether  their 
slopes  can  be  theoretically  predicted,  as  was  the  case  in  Figure  17, 
whether  trese  discriminants  are  sensitive  to  the  ambient  stress  level 
or  to  the  geological  medium  at  the  source,  etc. 

4)  A  study  should  be  made  of  the  influence  exerted  upon  the  discrimina¬ 
tion  by  the  geographical  distribution  of  events  used  in  the  training 
set.  It  has  already  been  noted  that  serious  questions  may  be  raised 
about  the  reliability  of  a  discriminant  function  which  was  calculated 
using  earthquakes  in  one  geographical  region  and  explosions  in  an¬ 
other.  Still  more  questions  may  be  raised  about  the  applicability 
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of  this  discriminant  function  to  unknown  events  which  occur  in  some 
third  region.  Only  when  a  quantitative  assessment  can  be  made  of 
the  various  source  region  and  propagation  path  effects  can  the  sta¬ 
tistical  discrimination  techniques  be  applied  with  confidence  to 
those  circumstances  in  which  they  would  most  be  needed  in  actual 
practice,  namely  explosions  detonated  at  previously  unknown  test 
sites . 

5)  An  effort  should  be  made  to  identify  and  analyze  anomalous  events 
which  are  especially  prone  to  misclassif ication. 

6)  A  discrimination  experiment  should  be  performed  which  uses  Pearce's 
(1977)  algorithm  based  on  the  detection  or  non-detection  of  pP  and 
sP.  This  experiment  would  determine  both  the  reliability  and  the 
magnitude  threshold  of  the  applicability  of  this  technique  relative 
to  that  of  the  statistical  techniques  investigated  in  this  report. 

7)  A  comparison  should  be  made  between  the  tele.eismic  discriminants 
used  in  this  report  and  regional  discriminants  based  upon  such 
phases  as  Pg  and  Lg.  Such  a  study  would  reveal  whether  the  use  of 
stations  located  within  the  Area  of  Interest  would  enhance  discrim¬ 
ination  as  well  as  detection.  This  study  should  use  events  in  the 
Southwestern  United  States  in  order  to  insure  uniformity  of  the 
regional  paths  for  signals  from  both  explosions  and  earthquakes. 


Pearce,  R.  G.  (1977).  Fault  plane  solutions  using  relative  amplitudes  of 
P  and  pP,  Geophys  J.  R.  Astr.  Soc.,  50,  381. 
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Network  estimates  of  parameters  measured  for  each  event 
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T»  EOH  EACH  REGION-TO-STATION  PATH 
assuming  PALL-OFF  AS  OHEG1-SOU ARED 


REGION 

1 

EVPNT 

ANNO 

ANTO  ATAK 

3FAK 

SO 

0.40 

0.40 

SB 

64 

0.20 

0.  SO 

143 

0.40 

0.  50 

0.  SO 

146 

0.20 

148 

0.3  0 

156 

0.60 

1 57 

0.50 

0.  SO 

1 58 

16S 

0.  10 

0.  SO 

169 

0.  40 

0.20 

0.40 

171 

0.40 

0.60 

176 

178 

0.30 

180 

0.  10 

182 

0 .  SO 

183 

0.20 

0.40 

186 

193 

194 

<T*> 

0.  333 

0. 2S0 

0.464 

SIGMA 

0.  122 

0.173 

0.  112 

RbClON 

EVENT 

36 

47 

2 

ANNO 

ANTO  ATAE 

BFAK 

0.  40 

0.  20 

0.40 

49 

0.60 

0.20 

55 

56 

0.  30 

57 

60 

0.30 

62 

144 

0.  10 

147 

0.  30 

0.  30 

154 

0.40 

164 

0.  40 

0.  60 

0.30 

<T*> 

0.475 

0.  367 

0.206 

SIGMA 

0.126 

0.209 

0.  107 

REGION 

3 

EVENT 

ANNO  ANTO 

ATAK  HPAK 

3 

0.40 

18 

0.4  0 

37 

0.  70 

6  3 

0.30 

0.  3  0 

269 

090  0. 30 

0.40 

<T*> 

0.  510  0.  300 

0.  450 

SIGMA 

0.  265  0.0 

0.  17  3 

REGION 

EVENT 

4 

ANNO 

ANTO  ATAK 

OFAK 

1b 

21 

0.  40 
0.20 

0.20 

o 

o 

6S 

270 

0.  SO 

0.20 

<T*> 

0.  367 

0.  200  0.200 

o.  ion 

SIGMA 

0.  153 

0.0  0.0 

0.0 

BOCO 

CH  TO 

CNAK 

0.50 

CT  AO 
0.  30 

0.  30 

0.  50 

0.  50 
0.  SO 

0.  60 

0.  50 

0.  40 

0.  30 

7.60 

0.50 

0.60 

0.  40 

0.  20 

0.  10 
0.20 
0.  20 

0.  40 

0.  90 

0.262 

0.  106 

0.S78 

0.  1  30 

7.  400 
0.  100 

30C0 

CHTO 

CNAK 

CT  AO 

0.30 
0.  20 

0.  30 

0.  50 

0.  40 
0.20 

0.  SO 

0.  30 
0.40 
0.40 

0.  40 

0.27S 

0.096 

0.  380 
0.094 

0.450 

0.071 

BOCO 

CHTO 

0.  20 

0.  3  0 

CNAK 
0.  40 

7.40 
7.40 
0.  30 

CTAO 

0.40 

0 . 2  SO 
0.071 

0.  37S 
0.  050 

0.400 

0.0 

noco 

CHTO 

CNAK 

CTAO 

0.  30 

0.  30 
0.  20 

0.300 

0.0 

0.  250 
0.071 

GONO 

HNNE 

IR7 

0.  40 

K  A  AO 

0.  30 

KONO 

KSRS 

0.  10 

0.70 

0.20 

0.  10 

0.  SO 

0.  10 

0.  10 
0.10 

0.20 

0.  20 
0.70 

0.  SO 

0.  10 

0.  SO 

0.20 

0.  40 

0.  42S 

0.  250 

0.  162 
0.074 

0.400 

0.219 

G'J  NO 

HNNE 

1P7 

K  AAO 

KONO 

KSRS 

0.  10 

0.30 

0.20 

0.40 
0.  40 

0.  20 

0.30 

0.  50 

oco 

OOO 

0.20 

0.  10 
0.20 
0.20 

0.  2  0 

0.40 

ON' 

occ 

OC 

0.  280 
0.  084 

0.200 

0. 063 

0.350 
0.  100 

GONO 

HNNE 

0.20 

0.20 

IR7 
0.  30 
0.  20 

0.  30 

K  A  AO 
0.70 

0.  10 

KONO 

KSRS 

0.  30 

0.  80 

0.200 

0.001 

0.  267 
0.  0S8 

0.  ISO 
0.071 

0.  SSO 
0.354 

GMNO 

HNNE 
0.  10 
0.  40 

I  R  7 

0.  20 

K  A  AO 
0.  30 
0.  10 

KONO 

KSRS 
0.  30 

0.20 

0.  250 
0.212 

0.  200 
0.  0 

0.700 

0.141 

0.  250 
0.071 
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T*  FOR  EACH  REGION-TO-STATION  PATH 
JSSIINTNG  FALL-OFF  AS  ONEGA— SQUARED 


REGION 

EVENT 

'lao 

NAIO 

50 

0.50 

0.30 

58 

64 

o 

o 

143 

146 

148 

156 

157 

158 

0.60 

165 

169 

0.  40 

171 

0.80 

176 

0.  30 

0.  10 

178 

180 

0.  80 

182 

0.  30 

183 

186 

193 

194 

<T*> 

0.525 

0.  200 

SIGH  A 

0.  198 

0.  141 

REGION 

2 

NAIO 

EVENT 

LAO 

36 

U  7 

0.60 

0.40 

49 

0.60 

55 

56 

57 

60 

62 

0.  50 

144 

147 

0.  30 

154 

0.20 

0.20 

164 

0.  40 

<T*> 

0.  433 

0.  300 

SIGMA 

0.  163 

0.  141 

REGION 

EVENT 

3 

LAO 

NAIO 

3 

18 

37 

63 

269 

0.  60 

0.20 
0.  JO 

<T*> 

SIGMA 

0.  600 
0.  0 

0.  250 
0.071 

H1JO 

NAO 

NJ  AN 

0.80 

0.80 

0.30 

ooo 

•  •  ■ 

OOO 

0.30 

0.  10 

ooooooo 

f-njnoJirn 

ooooooo 

0.20 

0.50 

0.30 

0.340 

0.  270 

0.431 

0.214 

NAJO 

NAO 

NJAK 

0.  60 

0.20 

0.50 

0.50 

0.40 

0.20 

0.  30 

0.10 
0.  30 
0.  10 

0.  50 

0.30 

0.400 

0.0 

0.250 
0.  187 

0.420 

0.110 

NAJO 

NAO 

N  J  AN 

0.60 

0.20 
0.  20 

0.  10 

0.30 

0.20 

0.50 

0.2°° 

0.001 

0.350 
0.  354 

0.333 
0.  163 

REGION 

EVENT 

16 

21 

65 

270 


<T«> 

SIGMA 


4 


LAO 


NAIO 

0.20 

0.  200 

0.0 


HA. 10 


NAO 
0.40 
0.  JO 


NJAK 
0.  JO 


0.J50  0.300 

0.071  0.0 


NS  AO 

RKON 

SHIO 

SNZO 

TATO 

TNAK 

0.40 

OCAK 

0.40 

ZOBO 

0.  40 
0.40 

0.  60 
0.4  0 

0.  30 
0.20 

0.  40 
0.50 
0.70 
0.40 

0.40 

0.  20 

0.40 

0.  40 
0.20 

0.60 

0.20 

0.70 

0.  50 
0.60 

0.90 

0.20 

0.50 

0.  200 
0.0 

0.  400 

0.  128 

0.  492 
0.207 

NWAO 

RKON 

SHIO 

SNZO 

TATO 

TNAK 

OCAK 

ZOBO 

0.  30 

0.40 

0.  10 

0.  20 

0.  60 

0.  20 

0.  40 

0.40 

0.  20 
0.30 
0.  30 

OOO 

fNCNV" 

OOO 

0.300 

0.0 

0.  250 
0.212 

0.267 
0.  082 

0.  340 
0.  167 

NWAO 

RKON 

SHIO 

SNZO 

TATO 

TNAK 

0.40 

OCAK 

ZO80 

0.  20 

0.  20 

0.  20 
0.40 

0.40 

0.40 

0.  200 
0.  0 

0.  200 
0.  0 

0.  333 
0.  115 

g.  400 
0.001 

NWAO 

RKON 

SHIO 

SNZO 

TATO 

TNAK 

UCAK 

ZOBO 

0.40 

0.20 

0.  10 

0.  20 

0.30 

0.  300 
0.14  1 

0.  100 
0.0 

0.  200 
0.  0 

0.  300 
0.  0 

t 

A 


* 


v 


* 
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T*  ?:95.  EACH  HEGTCW:TO-STmO»  PATH 


BEGION 

EVENT 

77 

145 

160 

<T*> 

SIGMA 

5 

ANNO 

ANTO 

AT  AN 

0.  30 

0.100 

0.0 

BEAK 

0.40 

0.30 

0.30 

0.333 

0.058 

BOCO 

CHTO 

0.  50 
0.20 

0.350 

0.212 

CNAK 

0.40 

0.30 

0.  350 
0.071 

CTAO 

GONO 

HNHE 

IP7 

0.30 

0.  10 

0.200 

0.141 

KAAO 

0.60 

0.10 

0.350 

0.354 

KOIO 

KSBS 

0.10 

0.50 

0.300 

0.283 

REGION 

EVENT 

6 

ANNO 

ANTO 

ATAK 

BFAK 

BOCO 

CHTO 

CNAK 

CTAO 

GUNO 

HNNE 

IB7 

KAAO 

KONO 

KSBS 

6 

0.  30 

26 

195 

<T»> 

0.30 

0.30 

0.40 

0.20 

0.300 
n  n 

0.300 

0.  0 

0.  400 

0.  0 

0.300 

0.  0 

0^200 

0.0 

SIGNA 

u  •  u 

REGION 

EVENT 

7 

7 

ANNO 

ANTO 

ATAK 

BEAK 

0.60 

DOCO 

CHTO 

CNAK 

0.20 

CTAO 

GONO 

HNNE 

IR7 

0.10 

KAAO 

KOHO 

KSRS 

0.40 

24 

27 

29 

30 

31 

32 

34 

35 

0.40 

0.  10 

0.60 

0.  40 

0.50 
0.  50 
0.  30 
0.  40 

0.  30 

0.  30 

0.  90 
0.70 

0.  50 
1.00 

0.20 

0.60 

0.40 

0.20 

0.10 

0.20 

0.30 

0.30 

0.50 

0.20 

39 

46 

0.40 

0.60 

0.  40 

6B 

0.  10 

75 

0.30 

0.20 

0.10 

0.10 

78 

151 

159 

166 

190 

0.80 

0.50 

0.60 

0.60 

0.50 

0.  40 
0.40 
0.60 

1.00 

191 

192 

CTO 

SIGNA 

0.600 

0.283 

0.514 

0.  121 

0.4  20 
0.0B4 

0.533 
0.  242 

0.280 

0.217 

0.245 
0.  129 

0.533 

0.416 

REGION 

EVENT 

8 

ANNO 

ANTO 

ATAK 

BFAK 

BOCO 

CHTO 

CNAK 

CTAO 

g  nno 

HNNE 

IR7 

0.20 

KAAO 

0.10 

KONO 

KSRS 

69 

66 

0.  10 

0.20 

7  3 
74 

149 

150 
153 
170 

0.  40 

0.60 

0.40 

0.60 

0.40 

0.20 

0.40 

0.  30 

0.  50 

0.  30 

0.  30 

0.60 

0.10 

0.70 

0.  40 

172 

0.  50 
0.20 

173 

175 

0.50 

0.40 

195 

188 

272 

<T*> 

SIGH* 

0.467 

0.115 

0.30 

0.40 

0.440 

0.114 

0.300 

0.141 

0.  10 

0.  325 
0.171 

0.  300 
0.001 

0.  320 
0.  217 

0.  133 
0.058 

0.550 

0.212 

III-3 


T»  FOB  EACH  R  EG ION- TO- ST ATI  OH  PATH 
ASSOHING  FALL-OFF  AS  OBEGA-SQO ABED 


REGION 

eveit 

5 

LAO 

HAIO 

HAJO 

NAO 

HJAK 

,47? 

160 

0.50 

0.40 

0.40 

0.40 

0.50 

<T*> 
SIGH  A 

0.450 

0.071 

0.400 

0.0 

0.450 

0.071 

REGION 

EVENT 

6 

LAO 

HAIO 

HAJO 

NAO 

NJAK 

6 

25 

195 

0.  40 

0.50 

0.50 

0.40 

<T*> 

SIGHA 

0.400 
0. 0 

0.500 

0.0 

0.500 

0.0 

0.400 

0.0 

REGION 

EVENT 

7lao 

HAIO 

HAJO 

HAO 

NJAK 

7 

24 

27 

29 

30 

31 

32 

34 

35 
39 
46 
68 
75 
78 

151 

159 

166 

190 

191 

192 


1.00 


0.30  0.40 


0.80 


0.  60 


0.60 

0.50 

0.20 


0.50 

0.50 


0.20 


0.50 


0.50 


0.70 

0.10 

0.70 


0.80 

0.30 

0.30 

0.50 


0.60 

0.80 

0.50 


0.60 

0.90 


<T»> 

SIGHA 

0.  800 
0.283 

0.486 
0.  195 

0.400 

0.0 

0.450 

0.251 

0.589 

0.215 

REGION 

EVENT 

8 

LAO 

HAIO 

HAJO 

HAO 

NJAK 

59 

66 

73 

74 

149 

150 
153 
170 

172 

173 
175 
185 
188 
272 

0.  10 

0.50 

0.  40 

0.  10 
0.  10 
0.  10 

0.30 

0.  10 
0.60 

0.60 

0.10 

0.  10 

0.  10 
0.  10 

0.40 
0.  20 

0.  10 

0.30 

0.40 

0.50 

0.30 

0.40 

<T*> 

SIGHA 

0.440 

0.207 

0.  100 
0.  000 

0. 260 
0.114 

0.400 

0.082 

HU  AO  RICOH  SHIO 
0.20 

0.200 

0.0 


NU  AO  EICON  SHIO 


HU  AO  RICOH  SHIO 


0.  30 


0.  300 

0.0 


NU  AO  R  ICO  H  SHIO 


0.  60 


0.  600 

0.0 


SHZO  TATO 


THAK 


o. : 
o.: 


OCAK 

0.40 

0.30 


0.200  0.350 

0.001  0.071 


SHZO  TATO  THAK  UCAK 

0.60  0.40 

0.600  0.400 

0.0  0.0 


SHZO 


TATO 


THAK 

0.40 

0.  50 


OCAK 


0.50 

0.60 


0.50 

0.40 


0.  60  0.40 

0.  50  0.70 

0.40 

0.20 


0.  500 
0.082 


0. 462 
0.  151 


SNZO  TATO  TNAK  UCAK 


0.40  0.30 


0.40  0.50 

0.30 
0.60 
0.50  0.30 


0.433  0.400 

0.  058  0.141 


ZOBO 


ZOBO 


ZOBO 


70BD 


III-4 


T*  POE  EACH  S EG ION-TO-ST  ATI ON  PATH 
ASSOHIHG  PALL-OPF  AS  OMEGA- SQOAMBD 


REGION 

EVENT 


<T*> 

SIGMA 


9 

AMMO 
0.  10 

0.30 

0.30 

0.30 


0.10 

0.  30 
0.  30 
0.  10 


0.225 
0.  109 


ANTO 


REGION  10 
EVENT  ANNO 
9 
9 

10 
23 
26 
28 
38 
91 
95 
98 
67 

69 

70 
72 
76 

155 
161 
179 

<T*> 

SIGMA 


EVENT 

61 

167 

<T*> 

SIGMA 


ANTO 


ATAK 

0.  30 

0.20 

0.  20 


0.  10 
0.  30 
0.30 
0.  10 


0.219 

0.090 


ATAK 


BPA  A 

0. 10 
0.30 

0.10 

0.10 

0.30 

0.30 


0.30 

0.10 

0.90 

0.30 


0.  230 
0.  116 


BEAK 


BOCO 


BOCO 


CHTO 

CNAK 

CTAO 

Gnno 

RENE 

IR7 

0.30 

0.  10 
0.  30 

0.20 

0.  30 

0.90 

0.40 

0.50 

0.  10 

0.  30 

0.  30 

0.30 

0.40 
0.  20 

0.40 

0.30 

0.40 

0.20 

0.60 

0.20 

0.20 

0.20 

0.60 

0.40 

0.  30 

0.  20 

I.2B0 

0.212 

0.414 

0.300 

0.400 

1.095 

0.09  9 

0.  146 

0.  100 

0.  100 

CHTO 

CM  AK 

CTAO 

Grino 

RNHE 

IB7 

KAAO 

0.10 

0.20 

0.10 

0.10 

o.Ao 


MONO 


0.20 

0.20 

■0.10 


0.  150 
0.053 


KAAO 

0.  10 
0.  10 


KOMO 


EVENT 

22 

279 

<T*> 

SIGMA 


KSHS 

0.  10 
0.  30 
0.30 

0.  30 


0.90 
0.  90 
0.20 

0.20 

0.  30 


0.  20 

0.270 
0.0  95 


K3»S 


0.40 

0.70 

0.30 

0.  10 

0.20 
0.  10 
0.20 
0.30 

0.10 

0.  10 

0.  20 

0.20 

0.  20 

0.  30 

0.  10 

0.  20 

0.400 

0.0 

0.500 

0.283 

0.200 

0.0 

0.  100 
0.0 

0.217 

0.075 

0.  100 
0.000 

0.200 

0.001 

1 1 

ANHO 

ANTO 

ATAK 

BPAK 

0.10 

BOCO 

CHTO 

CNAK 

0.30 

CTAO 

GUflO 

HNHE 

IR7 

KAAO 

0.20 

0.  10 

KONO 

K  SR  S 

0.  100 
0.0 

0.  300 
0.0 

0.  150 
0.071 

12 

ANHO 

0.  30 

ANTO 

ATAK 

8FAK 

BOCO 

CHTO 

CNAK 

CTAO 

Gn«o 

HNME 

I R  7 
0.  10 

KAAO 

0.  10 

KONO 

KSRS 
0.  20 

0.  300 
0.0 

0.  100 
0.0 

0.  100 
0.0 

0.200 

0.0 

III-5 


T*  POR  EACH  REGION-TO-STATION  PATH 
ASSUMING  FALL-OPP  AS  OMEGA-SOU  ARED 


REGION 

EVENT 


14 

17 

20 

H 

81 

152 

187 

189 

264 

265 

266 

267 

268 
271 

276 

277 

278 


<T*> 

SIGMA 


REGION 

PVENT 


4 

9 

10 

23 

26 

28 

38 

41 

45 

48 

67 

69 

70 
72 
76 

155 

1M 

179 


9 

LAO 

M  A  10 

MAJO 

NAO 

NJAK 

NU  AO 

RKON 

SHIO 

SNZO 

TATO 

TEAK 

OCAK 

0.30 

0.40 

0.10 

0.  10 

0.  10 

0.30 

0.20 

0. 20 

0.  30 

0.  10 

0.  10 

0.20 

0.20 

0.30 

0.30 

0.  10 

0.30 

0.  10 

0.30 

0.  10 

0.10 

0.  20 

0.10 

0.50 

0.30 

0.20 

0.30 

0.  20 

0.  30 

0.  10 

0.40 

0.20 

0.  30 

0.10 

0.10 

0.50 
0.  30 

0.20 

0.  1  0 

0.30 

0.20 

0.20 

0.  30 

0.10 

0.10 

0.  20 

0.30 

0.20 

0.  40 

0.  30 
0.  10 
0.30 

0.10 

0.10 

0.  10 
0.30 

0.  30 
0.10 

0.10 

0.30 

0.  10 

0.20 

0.20 

0.360 

0.084 

0.20 

0.  170 
0.082 

0.200 
0.  141 

0.  100 
0.000 

0.233 
0.  103 

0.200 

0.100 

0.  10 

0.  171 
0.076 

0.  300 
0.0 

0.  100 
0.001 

0.211 

0.060 

0.225 
0.  104 

10 

LAO 

MAIO 

MAJO 

NAO 

NJAK 

NU  AO 

RKCN 

SHIO 

SNZO 

TATO 

TNAK 

UCAK 

<T*> 

SIGMA 


0.  10 


0.  10 
0.  10 


0.20 


0.20 


0.  10 

0.60 

0.30 

0.20 

0.30 

0.80 

0.  10 

0.20 

0.20 

0.10 

0.114 

0.  367 

0.039 

0.208 

0.287 

0.20  0.20 
0.20 


o'. &oi  o^ooi 


REGION 

EVENT 

11  „ 
LAO 

MAIO 

MAJO 

NAO 

NJAK 

NU  AO 

RKON 

SHIO 

SNZO 

TATO 

TEAK 

61 

167 

0.60 

0.10 

<T*> 

SIGMA 

0.600 

0.0 

0.  100 
0.0 

REGION 

EVENT 

12  „ 
LAO 

MAIO 

MAJO 

NAO 

NJAK 

NU  AO 

RKON 

SHIO 

SNZO 

TATO 

TNAK 

22 

274 


0.20 


UCAK 


UCAK 


<T*> 

SIGMA 


0.10 

0.  100 

0.0 


0.200 

0.0 


T«  FOR  EACH  REGION- TO -STATION  PATH 
ASSUMING  FALL-OFF  AS  CNEGA-SOU AREO 


UPGION 

EVENT 

13 

ANNO 

AMTO 

ATAK  BFAK 

BOCO 

CHTO 

CNAK 

CTAO 

s 

00 

162 

163 

160 

177 

184 

0.20 

0.  20 
0.60 
0.40 

<T*> 

SIGNA 

0.200 

0.  0 

0.  400 

0.  200 

REGION 

EVENT 

1 

19 

33 

273 

14 

ANNO 

0.70 

0.  30 

0.  30 

ANTO 

ATAK  BFAK 

0.10  0.20 

0.  30  0.  50 

0.40 

BOCO 

CIITO 

0.50 

0.  30 

0.  10 

CNAK 

0.  20 

0.  20 

0.  30 

CTAO 

<T*> 

SIGNA 

0.433 

0.231 

3.200  0.  367 

0.141  0.153 

0.300 

0.200 

0.  233 

0.  058 

REGION 

EVENT 

13 

LAO 

N  A 10 

NAOO  NAO 

NJ  AN 

NW  AO 

RKON 

SHIO 

8 

80 

162 

163 

168 

177 

184 

0.50 

0.50 
0.  40 
0.  60 

0.20 

0.  60 

0.  10 

0.  30 

0.60 
0.  20 

0.20 

<T*> 

SIGNA 

0.  500 
0.082 

0.  400 
0.283 

0.  200 
0.  141 

0.333 

0.231 

REGION 

EVENT 

14 

LAO 

NATO 

NAJO  NAO 

NJAK 

NW  AO 

RKON 

SHIO 

1 

19 

33 

273 

0.50 
0..  40 
0.50 

0.  40 

0.  30 

0.  30 

0.20 

0.  10 
0.  20  0.40 

0.  10 

0.50 

0.30 

0.  10 
0.20 

<T»> 

SIGNA 

0.467 
0.  058 

0.350 

0.071 

0.233  0.  250 

0.058  0.212 

0.300 

0.200 

0.  150 
0.071 

Gil  NO 


GMNO 


SNZO 


SNZO 


I1NNE 

IR7 

F  JULO 

KONO 

KSRS 

0.  10 
0.20 

0.  10 

0.  10 

ooo 

ooo 

0.  125 
0.050 

0.400 

0.100 

HNNE 

IF7 

ITAAU 

0.10 

0.  10 

0.30 

FONO 

K  SR  S 

oo 

r~=T 

oo 

0.  10 

0.  50 

0.  10 

0.  350 
0.071 

0.  100 
0.0 

0.  167 

0.  115 

0.309 

0.293 

tato 

TNAK 

I1CAK 

ZOBO 

OO 

oo 

0.  500 

0.  141 

TATO 

TNAK 

IICAK 

7080 

0.60 

0.70 

0.  30 
0.30 

o  oo 

*-  t\r 

o  oo 

0.600 
0.  0 

0.433 

0.231 

0.200 
0.  100 

III-7 


T*  FOR  EACH  REGION-TO-ST ATION  PATH 
ASSURING  FALL-OFF  AS  ONEGA-CUBED 


REGION 

EVENT 

60 

SB 

64 

143 

146 

14B 

156 

157 

158 
165 
169 
171 
176 
178 
180 
182 
183 
186 

193 

194 

<T*> 

SIGNA 


1 

ANNO 
0.  30 

0.  10 

0.20 

0.30 

0.40 

0.  40 
0.  30 

0.  30 
0.  10 


0.  267 

0.112 


ANTO 


ATAK 


0.  30 


BFA  A 
0.30 

0.50 

0.30 

0.10 

0.60 
0.  30 


0.10  0.50 

0.20  0.40 

0.60 


0.  10 


0.175 

0.096 


0.40 
0.3  0 


0.  391 
0.  151 


REGION 

EVENT 

36 

47 

49 

55 

56 

57 
60 
62 

144 

147 

154 

164 

<T*> 

SIGNA 


2 

ANNO 

0.20 

0.  30 


0.  10 

0.  20 

0.200 

0.  082 


ANTO 


ATAK 

0.  10 


0.  20 

0.  60 

0.  267 
0.  208 


BFAK 

0.20 

0.10 

0.  30 

0.3  0 

0.10 
0.  30 
0.30 

0.229 

0.096 


REGION 

EVERT 

3 

IB 

37 

63 

269 

<T*> 

SIGNA 


3 

ANNO 

0.20 

0.  20 

0.80 

0.  400 
0.  346 


ANTO  ATAK 


0.30 


0.300 

0.0 


BFAK 

0.20 

0.70 

0.30 

0.30 

0.  375 

0.222 


BOCO 

CHTO 

Cl!  AK 

0.  30 

CTAO 

0.  10 

GUNO 

HNNE 

IB7 

0.  20 

KAAO 

0.30 

KONO 

KSRS 
0.  10 

0.  10 

0.  ',0 
0.20 

0.30 

0.40 

0.  10 
0.10 

0.  40 

0.  50 

0.40 

0.50 

0.  10 
0.  60 

0.  20 

0.40 

0.40 

0.40 

0.  40 

0.  10 

0.  10 

0.  10 
0.20 

0.  10 

0.20 

0.  10 

0.  20 

0.  10 

0.  80 

0.  30 

0.  10 

0.  10 

0.  30 

0.  10 

0.175 
0.  1  16 

0.433 

0.  166 

0.  267 
0.153 

0.225 

0.126 

0.  137 

0,  074 

0.293 

0.194 

BOCO 

CHTO 

CNAK 

CTAO 

GUNO 

HNHE 

IP  7 

KAAO 

KONO 

KSRS 

0.  20 
0.  20 

0.  20 

0.  30 

0.10 

0.  10 

0.10 

0.  30 
0.  20 

0.  10 

0.20 

0.  30 
0.  20 

0.60 

0.  10 
0.  20 
0.  30 

0.  30 

0.30 

0.  30 
0.10 
0.  30 

0.  10 
0.  10 
0.10 
0.20 

0.  10 

0.20 

0.225 

0.050 

0.  260 
0.  152 

0.300 

0.001 

0.200 

0.141 

0.180 
0.  1  10 

0.  133 
0.052 

0.200 

0.082 

BOCO 

CHTO 

0.  10 

0.  20 

CNAK 

0.30 

0.  30 
0.  20 
0.  20 

CTAO 

0.  20 

GUNO 

HNNE 

0.  10 
0.  10 

IR7 
0.  10 
0.  10 

0.10 

KAAO 

0.20 

0.10 

KONO 

KSRS 

0.  10 

0.  70 

0.  150 
0.071 

0.  250 
0.058 

0.200 
0.  0 

0.  100 
0.001 

0.  100 
0.001 

0.  150 
0.071 

0.400 

0.424 

REGION 

EVENT 

4 

ANNO 

ANTO 

ATAK 

BFAK 

BOCO  CHTO 

CNAK 

CTAO 

GUNO  HNNE 

0.  10 

IP7 

KAAO 

0.20 

16 

21 

66 

270 

0.  20 
0.  10 

0.40 

0.20 

0.  10 

0.20 

0.  30 

0.  30 

0.  10 

0.30 

0.  20 

0.  10 

<T*> 

SIGNA 

0.233 
0.  153 

0.  2  00 
0.  0 

0.100 

0.0 

0. 20  0 
0.0 

0.  300 
0.0 

0.  200 
0.141 

0.200 
0.  141 

0.200 
0.  0 

0.  160 
0.071 

KONO  KSRS 
0.  10 


0.10 

0.  100 
0.001 


III-8 


T*  FOB  EACH  RBGION-TO-*ST ATIOi  PATH 
ASSURING  FALL-OPP  AS  ONBGA-COBID 


REGION 

EVENT 

50 

58 

64 

143 

146 

148 

156 

157 

158 
165 
169 
171 
176 
178 
180 
182 
183 
186 

193 

194 


LAO 

NAIO 

NAJO 

NAO 

0.  30 

0.  10 

0.70 

0.  40 

0.20 

0.  60 
0.  30 
0.  80 
0.  10 

0.70 

0.  10 


0.  10 


0.  1  0 

0.  10 


0.20 


NJAK 

0.60 

0.50 

0.30 

0.40 

0.  10 
0.20 
0.40 
0.20 
0.20 
0.30 
0.  10 


0.40 
0.  30 


NW  AO 


peon  shio  snzo  tato 


0.  10 


<T*> 

SIGNA 

0.412 

0.264 

0.  100 
0.001 

0.  260 

0.  251 

0.308 

0.150 

0.  100 
0.0 

REGION 

EVENT 

2 

LAO 

NA  10 

NAJO 

NAO 

NJAK 

N  W  AO 

RKON 

SHIO 

SNZO 

TATO 

36 

47 

49 

55 

56 

57 

0.50 
0.  40 

0.  20 

0.  30 

0.20 

0.30 

0.50 

0.  10 

0.20 

0.  10 

60 

62 

144 

147 

154 

164 

0.  40 

0.  10 
0.  10 
0.30 

0.  10 

0.  30 

0.  10 

0.10 
0.  10 
0.  10 

0.  10 

0.50 

0.30 

<T*> 

SIGPA 

0.  300 
0.167 

0.  150 
0.  071 

0.300 

0.0 

0.  150 
0.  084 

0.340 
0.  167 

0.100 

0.0 

0.  150 

0.  071 

REGION 

EVENT 

3 

LAO 

N  A  10 

NAJO 

NAO 

NJAK 

NUAO 

RKCN 

SHIO 

SNZO 

TATO 

3 

18 

37 

63 

269 

0.6  0 

0.  10 
0.20 

0.20 

0.20 

0.60 

0.10 

0.10 
0.  10 

0 .30 

0.  20 

0.  10 

<T*> 

SIGNA 

0.600 
0.  0 

r.  150 
0.071 

0.200 

0.001 

0.  350 
0.354 

0.  167 
0.  1 15 

0.200 
0.  0 

0.  100 
0.0 

REGION 

EVENT 

4 

LAO 

N  A  10 

NAJO 

NAO 

NJAK 

NU  AO 

RKON 

SHIO 

SNZO 

TATO 

16 

21 

65 

270 

0.20 

0.40 

0.30 

0.  10 

0.  40 
0.10 

0.  10 

<T*> 

SIGNA 

0.200 

0.0 

0.  350 
0.071 

0.  100 
0.  0 

0.  250 
0.212 

0.  100 
0.0 

TNAK  UCAB 


0.20 

0.30 

0.40 

0.  40 
0.60 

0.  30 
0.20 

0.  10 

0.  10 
0.30 
0.50 
0.20 

0.20 

0.  10 
0.40 

0.  10 

0.50 
0.  20 
0.70 

0.  50 
0.30 

0.90 

0.10 

0.20 

0.300 

0.  165 

0.350 
0.  250 

TNAK 

□  CAK 

0.  10 

0.30 

0.  10 

0.40 

0.40 

0.10 
0.  20 
0.20 

0.  10 
0.20 
0.10 

0.  1 B  3 
0.117 

0.220 

0.130 

TNAK 

UCAK 

0.  20 

0.  10 
0.  20 

0.40 

0.30 

0.  167 
0.  058 

0.350 
0.  071 

TNAK 

OCAK 

0.  10 

0.30 

0.  100 
0.0 

0.  300 
0.0 

ZOBO 


ZOBO 


ZOBO 


ZOBO 


1 1 1-9 


T*  fP5_E«H  ;«GIOII-TO-S«tIOj|  «TH 


REGION 

EVENT 

77 

145 

160 

<T*> 
SIGH  Ik 


REGION 

EVENT 

6 

25 

195 


<T*> 

SIGHA 


REGION 

EVENT 

24 

27 

29 

30 

31 

32 

34 

35 
39 
46 
6  8 
75 
7B 

151 

154 

166 

190 

191 

192 


<T*> 

SIGHA 


5 

ANNO 

ANTO 

AT  AN 

0.  20 

BEAN 

0.40 

0.  20 

0.  20 

BOCO 

CHTO 

0.40 

0.  10 

CNAK 

0.  30 

0.  30 

CT  AO 

0.200 

0.0 

0.267 

0.115 

0.250 

0.212 

0.  300 

0.  001 

6 

ANNO 

ANTO 

ATAK 

BFAK 

BOCO 

CHTO 

CNAK 

CTAO 

0.20 

0.  10 

0.  20 

0.200 
0.  0 

0.  100 

0.  0 

0.  200 
0.0 

7 

ANNO 

ANTO 

atak 

BEAK 

0.40 

BOCO 

CHTO 

CNAK 
0.  20 

CTAO 

0.  10 


0.20  0.20 

0.40 

0.30 


0.  20 

0.  20 
0.  20 

0.  10 

0.  20 


0.  20 
0.  90 
0.  70 

0.  30 
1.00 
0.  40 


0.10 


0.40 


0.  60  0.50 

0.40 


0.  20 
0.  30 
0.  60 


0. 400 
0.283 


0.  357 
0.  098 


0.180 

0.048 


0.  442 
0.  294 


G  U  HO 

BNHE 

IH7 

0.  20 

0.  10 

KAAO 

0.50 

0.10 

KONO 

0.  150 
0.071 

0.300 

0.283 

G  OHO 

HNHE 

IH7 

KAAO 

KONO 

0.20 

0.  10 

0.  200 
C.O 

6.  ioo 

0.  0 

G  OHO 

HKHE 

IR7 
0.  10 

KAAO 

KONO 

0.  10 

0.  10 
0.10 

0.  30 


0.  10 

0.10 

0. 10  0. 10 
0.40 

0.20 

0.  10 

0.20 

0.10 

0.10 

0.10 


0.  140  0.  145 

0.009  0.093 


REGION 

EVENT 

59 

66 

73 

74 

149 

150 
153 
170 

172 

173 
175 
135 
18R 
272 


R 

A  NHO 


<T«> 

SIGHA 


ANTO  ATAK 

JFAK  BOCO 

CHTO 

CNAK 

CTAO  GU  HO 

HNHE 

IR7 

0.  10 

KAAO  KONO 

0.  10 

0.  10 

0.  10 

0.  10 

0.  10 

0.  20 

9.30 

0.  10 

0.20 

0.  30 

0.  30 
0.  JO 

0.20 

0.  30 

0.  30 

0.20 

0.  20 

0.  20 

0.20 

0.  10 

0.267 

0.058 

0.  10 

0.40 

0.  240 

0.  114 

0.200 

0.141 

0.  10 

0.  200 
0.082 

0.  150 

0.071 

0.  160 
0.  089 

0.  100 

0.  001 

KSBS 

0.  10 
0.40 

0.250 

0.212 


KSRS 


KSRS 

0.  20 


0.  20 


0.  90 


0.433 

0.404 


KSRS 


0.50 


0.20 


0.350 

0.2  12 


T*  FOP  EACH  REGION-TO-STATION  PATH 
ASSUMING  EALL-OEE  AS  OHEGA-COBED 


REGION 

EVENT 

5 

LAO 

1 A  10 

PI  A  JO 

NAO 

NJ  AA 

NW  AO 

77 

145 

UO 

0,40 
0.  20 

0.20 

0.  10 
0.20 

<T»> 

SIGHA 

0.  300 
0.  141 

0.  200 

0.  0 

0.  150 

0.  071 

REGION 

EVENT 

6 

LAO 

H  A 10 

MAJO 

NAO 

NJAA 

NW  AD 

6 

25 

1  45 

0.  30 

0.50 

0.  30 

0.  30 

<T*> 

SIGHA 

0.  300 
0.0 

0.500 

0.0 

0.300 

0.0 

0.3  00 
C.O 

REGION 

7 

EVENT 

LAO 

HATO  HAilO 

NAO 

NJAK 

NW  AO 

7 

0.  10 

o 

o 

24 

1.00 

27 

0.  30  0.  30 

0.30 

29 

30 

0.20 

0.20 

31 

0.30 

32 

0.70 

34 

35 

0.  30 

0.40 

39 

0.20 

0.80 

46 

0.  30 

6B 

0.  10 

0.  30 

75 

7fl 

151 

0.20 

0.  60 

0.  30 

159 

0.10 

166 

0.  40 

O 

*M 

O 

0.60 

190 

0.  bO 

191 

192 

<T*> 

0.700 

0.  206  0.300 

0.  300 

0.411 

SIGHA 

0.424 

0.195  0.0 

0.210 

0. 1  90 

REGION 

8 

EVENT 

LAO 

HA  I  O 

HAJO  NAO 

N.7AK 

NW  AO 

59 

0.  10 

66 

71 

0.  30 

0.  10 

0.  30 

74 

0.  10 

149 

o 

o 

0.  10 

150 

153 

170 

0.  10 

0.  10 

0.  20 

0.  20 

172 

0.  50 

0.  10 

0.  30 

173 

0.10 

175 

0.  40 

O 

o 

0.  10 

186 

0.  10 

IRA 

0.  10 

0.  20 

272 

0.  20 

<T*> 

0.320 

0.  100 

0.180 

0.200 

SIGHA 

0.  1  4  fl 

0.  000 

0.084 

0.  OR 2 

RKON 

SHIO 

SNZO 

TATO 

THAN 

UCAA 

ZOBO 

0.  20 

0.  10 

0.2U 

0.  10 

0.  30 

1. 200 

0.  100 

0.  250 

1.0 

0.001 

0.071 

RKCN 

SHIO 

SNZO 

TATO 

TNAK 

0CAK 

ZOBO 

0.40 

0.30 

0.400 

0.  300 

0.  0 

0.0 

RKGN 

SHIO 

SNZO 

TATO 

THAN 

IICAK 

ZOBO 

0.  20 

0.20 

0.  30 

0.20 

0.40 

0.  30 

0.  30 

0.40 

0.30 

0.40 

0.40 

0.30 

0.10 

l.  200 

0.325 

0.287 

.0 

0.096 

0.099 

RKON 

SHIO 

SNZO 

TATO 

THAN 

IICAK 

ZOBO 

0.  30 

0.20 

0.40 

0.20 

0.20 

0.  10 
0.30 

0.  30 

0.10 

.400 

0.  267 

0.  180 

'.  0 

0.058 

0.084 

III-ll 


T*  rOHjlCHfWIOH-TO-STKIOJJ  MTB 


REGION 
EVENT 
14 
17 
20 
53 
79 
R  1 

152 

187 

189 

264 

265 

266 

267 

268 

271 

276 

277 

278 

<T»> 

SIGMA 


REGION 

EVENT 

4 

9 

10 

23 

26 

28 

38 

41 

45 

48 

67 

69 

70 
72 
76 

155 

161 

179 

<T*> 

SIGMA 


REGION 

EVENT 

61 

167 

<T*> 

SIGMA 


REGION 

event 

22 

274 

<T»> 

SIGMA 


9 

AN  MO 
0.  10 

0.  10 
0.  10 

0.  10 


0.  10 

0.  20 
0.  10 
0.  10 


0.112 

0.035 


10 

ANMO 


11 

ANMO 


12 

ANMO 
0.  10 


0.  100 

0.  0 


INTO 

ITAK 

BFA  K 

0.  10 

0.  30 

BOCO 

ONTO 

0.30 

0.30 

CIAE 

0.  10 
0.  30 

0.  30 

0.  1C 

0.  10 

0.  20 

0.  10 

0.  30 

0.30 

0. 10 

0.  10 
0.20 

0.  10 
0.  10 

0.  10 

0.30 

0  10 

0.  10 

0.30 

0.1  0 

0.  20 

0.  10 

0.40 

0.20 

0.  30 

0. 20 

0.186 

0.090 

0.  190 
0.  110 

0.260 

0.089 

0.162 

0.092 

ANTO 

ATAK 

3FA  K 

DOCO 

CHTO 

CNAK 

0.  30 

0.70 

0.  10 

0.  10 

0.  300 
0.0 

0.  400 
0.424 

0.100 

0.0 

ANTO 

ATAK 

BFAK 

0.10 

BOCO 

CI1TO 

CNAK 

0.20 

0.  ICO 

0.  0 

0.200 
0.  0 

ANTO 

ATAK 

BFAK 

BOCO 

CHTO 

CNAK 

CTAO 

0.20 

0.20 

GOMO 

HN  HE 

0.  10 

0.30 

IR7 

0.50 

0. 10 

KAAO 

0.10 

0.10 

0.10 

0.10 

KONO 

0.20 

0.20 

0.40 

0.10 

0.10 

0.  10 

0.50 

0.  10 
0.50 
0.20 

0.10 

0.20 

0.10 

0.271 

0.160 

0.200 

0.100 

0.333 

0.208 

0.110 

0.032 

CTAO 

GO  HO 

hnhe 

IR7 

KAAO 

KONO 

0.  10 

0.10 

0.  10 

0.  10 

0.  10 

0.  10 
0.20 

0.10 
0.  10 

0.  10 

0.  10 

0.10 

0.100 

0.0 

0.117 
0.04  1 

0.100 

0.000 

CTAO 

GOMO 

HNHE 

IR7 

KAAO 

0.20 

0.10 

KONO 

0.  150 
0.071 

CTAO 

GOHO 

HNHE 

IR  7 
0.10 

KAAO 

0.  10 

KONO 

0.100 

0.0 

0. 100 
0.0 

KSBS 

0.  10 
0.20 
0.  30 

0.  20 


0.  40 
0.40 
0.  20 

0.20 
0.  20 


0.  10 

0.230 

0.106 


K  SR  5 


0.  10 


0.  10 

0.100 

0.001 


KSRS 


K  SR  S 
0.  10 


0.  100 
0.0 
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T*  FOR  EACH  REGION-TO-STATION  PATH 
ASSOHIN  G  FALL-OFF  AS  OHEGA-CtlBED 


EVEHT 

LAO 

HAIO 

HAJO 

NAO 

NJAK 

NW  AO 

RKON 

14 

0.  10 

0.  10 

0.  10 

0.30 

0.  10 

17 

20 

0.  30 

0.  10 

0.  20 
0.  10 

53 

0.20 

0.  10 

0.30 

0.  10 

0.10 

0.  10 

79 

0.30 

0.  10 

a  i 

0.  20 

0.20 

152 

0.10 

187 

189 

0.20 

0.  10 

0.  10 

0.  20 

264 

0.  30 

0.  10 

265 

0.30 

0.  10 

266 

0.20 

0.  10 

0.1  0 

0.10 

267 

0.30 

0.  30 

268 

0.  10 

0.10 

0.  20 

271 

0.30 

0.  10 

0.  20 

0.  10 

276 

0.  10 

277 

0.  10 

278 

0.  20 

0.  10 

<T*> 

0.  240 

0.  150 

0.200 

0.  100 

0.  167 

0.133 

0.  129 

SIGH  A 

0.070 

0.085 

0.141 

0.00  0 

0.082 

0.058 

0.049 

REGION 

10 

EVEHT 

LAO 

HAIO 

HAJO 

NAO 

NJAK 

nwao 

RKON 

4 

0.  10 

9 

10 

23 

26 

0.  10 

28 

0.  10 

38 

41 

0.  10 

45 

48 

67 

0.  10 

0.60 

0.30 

69 

70 

0.10 

72 

0.  10 

0.  80 

76 

155 

161 

0.  10 

0.10 

179 

0.10 

0.  10 

<T*> 

0.  100 

0. 2b7 

0.325 

SIGHA 

0.000 

0.  289 

0.330 

REGION 

1  1 

EVENT 

LAO 

HAIO 

HAJO 

NAO 

NJAK 

NWAO 

RKON 

61 

0.  30 

0.  10 

167 

<T*> 

0.  300 

0.100 

SIGHA 

0.0 

0.0 

REGION 

12 

FV  ENT 

LAO 

HAIO 

HAJO 

NAO 

NJAK 

NWAO 

RKON 

TATO 

TNAK 

OCAK 

0.10 
0.  20 

0.30 

0.10 

0.  10 

0.10 

0.10 

0.  10 

0.30 

0.  10 

0.  10 

0.  10 

0.20 

0.20 
0.  20 

0.30 

0.10 

0.  1  0 

0.10 

1.  100 

0.133 

0.  187 

1.001 

0.  050 

0.099 

TATO 

TNAK 

UCAK 

0.10  0.10 
0.10 


0.  100  0.  100 
0.001  0.001 


S  NZO  TATO  TNAK  OCAK  ZOBO 


S  NZO  TATO 
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mmm 


T*  r OR  EACH  REGION-TO-STATION  PITH 
ASSUMING  FALL-OFF  AS  OHEG A-C0BED 


REGION 

EVERT 

a 

80 

162 

163 

168 

177 

iau 

13 

UNO 

INTO 

ATAK  BEAK 

HOCO 

CHTO  CM  IK 

CTIO 

GONO 

HUBS  IR7 

KUO 

KORO 

KSRS 

0.  10 

0.  10 

0.  50 

0.  30 

0.  10 
0.  10 

0.  10 
0.  10 

0.10 

0.10 

0.50 

<T*> 
SIGH  A 

0.  100 

0.  0 

0.300 

0.  200 

0.  100 
0.001 

0.233 

0.  231 

REGION 

14 

EVERT 

!HHO 

!NTO  !TAK 

BEAK 

BOCO  CRTO 

CNAK 

CTAO 

GtJ  HO  HNHE 

IB7 

KAAO 

KORO  KSRS 

1 

0.70 

0.  10 

0.10 

0.30 

0.  10 

0.10 

19 

0.  10 

0.  20 

0.20 

0.  10 

0.10 

0.  40 

33 

0.  30 

0.40 

0.20 

0. 20 

273 

0.  10 

0.20 

0.  10 

0.  10 

0.20 

0.  10 

<T*> 

0.300 

0.200 

0.233 

0.200 

0.  133 

0.  200 

0.  100 

0.  133 

0.250 

SIGH! 

0.346 

0.141 

0.  153 

0.100 

0.058 

0.  001 

0.0 

0.058 

0.212 

REGION 

13 

EVENT 

LAO 

HAIO 

FIAJO  NAO 

NJAK 

RW  AO 

RKON 

SHIO 

SRZO 

TATO 

TRAK  OCAK 

ZOBO 

a 

80 

162 

0.  30 

0.10 

0.  10 

0.40 

0.  40 

163 

0.  10 

0.  30 

168 

0.  30 

177 

0.30 

184 

0.  30 

0.60 

0.  30 

0.20 

<T*> 

0.300 

0.  350 

0.  200 

0.233 

0.  350 

SIGN! 

0.001 

0.354 

0.  141 

0.  153 

0.071 

REGION 

14 

EVENT 

LAO 

HAIO 

HAJO 

NAO 

NJAK 

NW  AO  RKON 

SHIO 

SNZO  TATO 

TRAK 

OCAK 

ZOBO 

1 

0.  30 

0.40 

0.20 

0.  10 

0.50 

0.  10 

19 

0.  30 

0.  10 

0.  10 

0.  60 

33 

0.50 

0.  10 

0.20 

0.  10 

0.  10 

0.  10 

27  3 

0.  10 

0.20 

0.20 

0.20 

0.  20 

0.20 

<T*> 

0.367 

0.250 

0.  167 

0.  150 

0.  167 

0.  100 

0.600 

0.267 

0.133 

SIGHA 

0.  115 

0.  212 

0.058 

0.071 

0.058 

0.001 

0.  0 

0.208 

0.058 

4 

d 
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APPENDIX  IV 


Discrimination  variables  calculated  for  each  event 


VALUES  OF  THE  20  DISCI;  If*.  I  NAT  I  OK  VARIABLES  IN  TABLE  TIV 


ri *- ou r-. «-  p-iin 3  (N 

rjrj  XOO'lT  *“CNO  ^  ' 

cr.  ^irr“.comiriVX1'^,X3C'OOC  ^c^f-crn 

r^rriplT'CCCOOCOOC  OOINOCC1 

mi  i  i  1  1 


r* cnr  cl  o ir> ^ olt  u^ir  u  irir u-  r~ oct'ct  & 
co^fx  fN^O'0'vc=t  vixr*rs'iT(N'-»rxr'fNiC 

r-  •7'tnr*'rrTf*~Ov-r-'fNJcr(N»-0'^(N=Tf^(Nir 

^  O^ClTir,  CC  C*“  «“  v-fNfsjJNC  CCCCC 

Tit  I  . 


lt fMrr r-  u  u  > 

rj(Ncr  *-c  *-  M^roo  pH.(£j#“  (Nr.’irir^a  cm 
oaf^-t-.v£)vrvrC7»-rJr*‘*-(^'»^f'J0'rv‘CLa' 

-  r-'  .-'  J  J.  o' coco  c  o'  o'  -V'  rn  <N  rn  r!  r3 

III  I  I  I  I  I  I  I  I  I  I  I 


rolfrO'TIN  rMrlpj  l/-  IT  »*■  I-  f  I- 
*-(Nr".(NO"£  »—  O'  3  fNP-lQ-;  rnp^r^CtD^ 
u.irc  irr** c cNr-xcr-  ar- nn  -j 

Of-cnaoc  occr-ocof“ 
til  ill  i  l  i  I  i  i 


(N  fT  (N  CD  rvj  C7>  ri  v£  (N  r>i  00  IN  (N  CD  IT  IT.  fN  r .  rsj  r  j 
irfNirura'f-,r*c?ofNzjr-r‘3.-rcr 
Vf  r^X-^rtrt-XCNO.OXX^OX 
*“  ^ircoo'oocooooniNiNnnfN 

I  I  I  I  I  I  I  1  l„l  I  I  I  II  1 


r-  -r  *-C'ri  *-  cT'tro'  irinin 

f-^-rrvc^D  &  ct*“(n  ^r,_ 

cr  i/-cor-vccc,crccon3hOCCoo'0'CD 

^  of-V^cccoo«-cocooci'iriiN 
lit  ■  1  1  1 


-J  ^  fvj,-  r*»  r-  rN  <N  <N  CO  (N  CDCDlT  irirsHNfN  CN 

ccr-izjirfMfr(T«eNfr3r^fOiN=J  *~r*  r**-  =tcr 
- f  ror^^CiTO  srr*»(N  0^*00  IN 

rr-riT^ooo^ococ-nr-.nnr^ 

III  I  I  I  I  I  I  I  I  I  I  I  I 


^  ^  jw  q  p-  vD  ^ 1 3  X 

^  ^  pj  p  pg  pp  p*  l/* i  3  CD  'C  M 

r»  tta.th®r-Mhcocoooccc303ir® 

^  oCClT^OOCOCCOCOCOrOCO 

I .  . 


^  r- CD  *- r- O  IT  O'  COCO  lT.lT  or*> 

s£r*iLT  (NfNfNv— cr>  a>r*  fs*r"  O'Jp 

q  ocoo(NOu''iroa'coocc«r)r'Cir(Nc 
*""  dro^.coorco-occoocc-o 

I  ill  l  I  l  I 


sI3f_f-^f_  co*-xxr. 

Ttccno'*-  o^-rncNcn 

VC  *— vOfMXO  OGOC.GOOC  OO O*-  o 

psj  , . . 

COOITCCOCOCOCOOOOOOOO 

II  III 


ir.CT'f^ir'  lT>  &  &  i D  3  3  ^ 

cn^o  o m  r*  m ct'Cnco  ^  -t 

^  rur-O'OCOMMvCCO'OCOlTf-OOOO 
!•••••■•  •••••  ••••*** 

r-^OCOOOr''-OMOOOOOCCCO 

I  II  . 


irir»rsi^®  ^miT.  r^r^r^r'r'r- 

LpZ}r"D5®r^rslfN  X*-  Tj  CTiTO' 

PPPP3  ^lPM»-*‘OOOOCOMPnMM3  3 

occirinoococcocooccrco 

II  III 


v£>  IT  vC  COX  ° 

<n  r*  *-  c- 

CD  O3OCO**O0'OOOCCOO3OOnC 
ddoOOOOCOOOOOOOININOINO 

I  II  III 


n«NC<^  xod*-  o-cn  ronccr'r-»“ 

OLD  IT  CM  v£>  o  o  cnr-  *—  nO^ClTiM 

U'lNf-Monnr-oocrr-unorocn 

droXCOOOOCOININOOOOOOO 

I  I  I  I  I  I 


nr-hhconQCvonnnnrn  cotN^nr 
mfNco^,D»-,T)^LD'r'»-fN  (Nr-rixr^o 
p-.  cr^cT'irr-»-fNr'zj  iDr-r-vra'ocrirvD'^r*vr 

dodaiPOOOOCrOOO— --"■ 

Ml  I  I  I  I  I  I  I  I 


occooaxocooooooooooo 

o'r-OOOC-OOCOOOOCOOOOO 

I  I  I 


c  ooocco 

^  COvC^ont- 

oocococoooc'-ovooirirn  oo 

doc  o  co  oo  c  o  o  o  o  cr  *-  *-*-»-  *-o 


«-  m  cn  n »-  o  m  *-  cm  cn  oo  <n  cn  cc  <n  cm  ip  irur  ir> 

vuooc7'^DcjrNfN»-r^Lrcr'OO^x^JO'DCN 
(M  CTsojr**^£>cr'®«T'Oit  cn*- rrrsifN®ors)p-ip-io 

M  I  i  I  I  •  I  ••••••  . . *  * 

OfN*-LDiT,OC*-OOOOOOrnornrnr-irr 
Ml  II  I  I  I  II  I  I  I  I  I  I 


PP  rc  z}  OCCOO  cra'OCsirMCM 

cn  r-®  t-inrr  ®cn  *-m®r*ccN 

ocmo®  *- o  o  ocr*  ®  ^3  o' (no  r- ®  ox  «o  x 

c  *“  o  IT.  vt.  c  O  c*  O  C  *“  c  »-  O  *-  *“  *-  •“  *“  *“ 


r^oxcr-c  3  rsirrx®  cmxXcnit  (ncmcn 
^  ID st  Ct  O' 1 ^■NO'r'hr1'  XOOlTi 
*-  's443rnxcr'Lf'iri®<,,ncN®ro*-OLrccr^irfn 

*VJ  . . . 

rrriTirOOOCOorCrtNOO^nf* 

III  I  I  I  I  I  II  I  III 


-f  (Mr-vO(NrjlP(N  (N(N(N(N^  tr~’-crcr 
fj.iN(rnu-  u'C'-r-i'-'-'fiopo 

v0^^min<0'-"'^o'^a)0.-r»-'-' — ic*- 

OOCITU'OOO'-Ol'IfMM’-OOOCO 

II  III  III 


3  0'^INC'0  3V"£'0\£^)^^r'l''0^^3 
r-r'onfNr^iririfNcor^rsif-r^oirirpsi^*- 

o  ssovio  aainr~ircM/'3(N«-cD3eoin,f' 
^  ^(NiN^r-'o’oO'-'-O'-'-oiractna-a 

ill  1  1  1  1  1  1  1  1  I  1  1  1  1  1  1 


r'-.fNrr-.irr-'r-'^*  ^  *-»-rsooo 

(Nu-u- ^>-ir  a  o  a  Lri^r^.-ri^o 

c?tN'~ir  rin^o^'-ccrno'i-  i~ter»u- 

^T-'oVJoc'odo'-oconMN"r*r< 

III  Ml  II  1  I  I  I 


cOvi-’-nacNiN^cca:  rccaa  ni-rn^r 

^  r-roinO=T  OL7'‘7'OP*'fNir^'“CDII'  rTC*- 

^  3oh*-  =»ir  ro^nri^or'irrrr-i^op** 

dddddd'iO'-'-ot-coai-mi'' 
ill  1  1  1  I  I  I  1  1  I  1  1  I  I  I  1 


— in'"',  itinvor-  mtr  O'-iNrrair  jir-  in  oo 


[H  ^iMr- ctL^^nr* rro'o rnn  =jit  cr'cooc^ 


rv-i 
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values  OF  THE  20  DISC9I*1IHATICN  VAPTABLES  IK  TABLE  XIV 


IV- 3 


VALUES  CF  THE  20  DISC  H  IF  I NATIOK  VAHIARLFS  IN  TABLE  XT'* 


oa— rsixp-  z?  o 

*— r-ocr  lp  zj  cr 

<?•  ncfsocccnorPccococcoc 

ir.  . . .  . . . . 

•“  Oi-0»ChC(\rOCOOCOOOOCCO 
I  t  I 
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